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Preface 



Presently, the notion of function is not as finally crystallized 
and definitely established as it seemed at the end of the 19th cen- 
tury; one can say that at present this notion is still in evolution, 
and that the dispute concerning the vibrating string is still going on 
only, of course, in different scientific circumstances, involving other 
personalities and using other terms. 

Luzin N.N. (1935) [42] 

It is symbolic that in that same year of 1935, S.L. Sobolev, who 
was 26 years old that time, submitted to the editorial board of the 
journal “Matematicheskiy sbornik” his famous work [61j and pub- 
lished at the same time its brief version in “Doklady AN SSSR” M- 
This work laid foundations of a completely new outlook on the con- 
cept of function, unexpected even for N.N. Luzin — the concept of 
a generalized function (in the framework of the notion of distribu- 
tion introduced later). It is also symbolic that the work by Sobolev 
was devoted to the Cauchy problem for hyperbolic equations and, in 
particular, to the same vibrating string. 

In recent years Luzin’s assertion that the discussion concern- 
ing the notion of function is continuing was confirmed once again, 
and the stimulus for the development of this fundamental concept 
of mathematics is, as it was before, the equations of mathematical 
physics (see, in particular, Addition written by Yu.V. Egorov and 
|10L 1111 116L 1171 1181 1321 1491 167] j. This special role of the equa- 
tions of mathematical physics (in other words, partial differential 
equations directly connected with natural phenomena) is explained 
by the fact that they express the mathematical essence of the funda- 
mental laws of the natural sciences and consequently are a source and 
stimulus for the development of fundamental mathematical concepts 
and theories. 
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viii PREFACE 

The crucial role in appearance of the theory of generalized func- 
tions (in the sense the theory of distributions) was played by J. Ha- 
damard, K.O. Friedrichs, S. Bochner, and especially to L. Schwartz, 
who published, in 1944-1948, a series of remarkable papers concern- 
ing the theory of distributions, and in 1950-1951 a two- volume book 
[54], which immediately became classical. Being a masterpiece and 
oriented to a wide circle of specialists, this book attracted the at- 
tention of many people to the theory of distributions. The huge 
contribution to its development was made by such prominent math- 
ematicians as I.M. Gel’fand, L. Hormander and many others. As a 
result, the theory of distributions has changed all modern analysis 
and first of all the theory of partial differential equations. Therefore, 
the foundations of the theory of distributions became necessary for 
general education of physicists and mathematicians. As for students 
specializing in the equations of mathematical physics, they cannot 
even begin any serious work without knowing the foundations of the 
theory of distributions. 

Thus, it is not surprising that a number of excellent monographs 
and textbooks (see, for example, (12l [22l [23l [25l [3TI l40l [44l [54l 
1571 [5911621 [681 169] 1 are devoted to the equations of mathematical 
physics and distributions. However, most of them are intended for 
rather well prepared readers. As for this small book, I hope that it 
will be clear even to undergraduates majoring in physics and math- 
ematics and will serve to them as starting point for a deeper study 
of the above-mentioned books and papers. 

In a nutshell the book gives an interconnected presentation of 
a some basic ideas, concepts, results of the theory of generalized 
functions (first of all, in the framework of the theory of distributions) 
and equations of mathematical physics. 

Chapter [l] acquaints the reader with some initial elements of the 
language of distributions in the context of the classical equations 
of mathematical physics (the Laplace equation, the heat equation, 
the string equation). Here some basic facts from the theory of the 
Lebesgue integral are presented, the Riesz spaces of integrable func- 
tions are introduced. In the section devoted to the heat equation, 
the student of mathematics can get familiar with the method of 
dimensionality and similarity, which is not usually included in the 
university program for mathematicians, but which is rather useful on 
the initial stage of study of the problems of mathematical physics. 
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IX 



Chapter [2] is devoted to the fundamentals of the theory of distri- 
butions due to L. Schwartz. Section [16] is the most important. The 
approach to some topics can also be interesting for the experts. 

Chapter 3 acquaints the reader with some modern tools and 
methods for the study of linear equations of mathematical physics. 
The basics of the theory of Sobolev spaces, the theory of pseudodif- 
ferential operators, the theory of elliptic problems (including some 
elementary results concerning the index of elliptic operators) as well 
as some other problems connected in some way with the Fourier 
transform (ordinary functions and distributions) are given here. 

Now I would like to say a few words concerning the style of the 
book. A part of the material is given according to the scheme: defini- 
tion — theorem — proof. This scheme is convenient for presenting 
results in clear and concentrated form. However, it seems reason- 
able to give a student the possibility not only to study a priori given 
definitions and proofs of theorems, but also to discover them while 
considering the problems involved. A series of sections serves this 
purpose. Moreover, a part of the material is given as exercises and 
problems. Thus, reading the book requires, in places, a certain ef- 
fort. However, the more difficult problems are supplied with hints 
or references. Problems are marked by the letter P (hint on Parking 
for the solution of small Problems). 

The importance of numerous notes is essentially connected with 
a playful remark by V.F. D’yachenko: “The most important facts 
should be written in notes, since only those are read” . The notes are 
typeset in a small font and located in the text immediately after the 
current paragraph. 

I am very grateful to Yu.V. Egorov, who kindly agreed to write 
Addendum to the book. I would also like to acknowledge my grati- 
tude to M.S. Agranovich, A. I. Komech, S.V. Konyagin, V.P. Palam- 
odov, M.A. Shubin, V.M. Tikhomirov, and M.I. Vishik for the useful 
discussions and critical remarks that helped improve the manuscript. 
I am also thankful to E.V. Pankratiev who translated this book and 
produced the CRC. 

While preparing this edition, some corrections were made and 
the detected misprints have been corrected. 



A.S. Demidov 



Notation 



N = {1,2,3,...}, Z = {0, ±1, ±2, . . . } and Z+ = {0, ±1, ±2, . . . } are 
the sets of natural numbers, integers and non-negative integers. 

X x Y is the Cartesian product of the sets X and Y, X n = X n ~ 1 x X. 
i = y/—l is the imaginary unit (“dotted i”). 
i = 2m the imaginary 2tt (“i with a circle”). 

R” and C n are n-dimensional Euclidean and complex spaces; R = R 1 , 
C = C 1 9 z = x + iy, where x = $tz £ R and y = < Az £ R. 
x<y, x<y, x>y, x>y are the order relations on R. 
a 1 means the a sufficiently large. 

{a: £ X | P} is the set of elements which belong to X and have a 
property P. 

]a, b] = {x £ R | a < x < 6}; [a, b], ]a, 6[ and [a, b[ are defined similarly. 
{a ra } is the sequence {a n }“ L 1 = {ai, 02, 03 , . . . }. 

/ : X 9 x 1— > f(x) £ Y is the mapping / : X — > Y, putting into 

correspondence to an element x £ X the element f(x) £ Y. 

1 a is the characteristic function of the set A, i.e. 1^ = 1 in A and 

1^4 = 0 outside of A. 

7 r 

arccot a = arctan a. 

2 



x — > a means that the numerical variable x converges (tends) to a. 

=>■ means “it is necessary follows”. 

•<=>■ means “if and only if” (“iff’), i.e an equivalence. 

A <e S! means A is compactly embedded in Q (see Definition 3.2). 
C7 m (fi), CT"(fi), C m (Cl), PC m (D), PC b m (D), C£*(fi), CJ*(ft) see Defi- 



9.1 



nition |3iT] (for 0 < m < 00 ). 

L p (fl), P°°(S2), Zf oc (D) see Definitions L |L 

P^Q), P#(fl), P(f2), P'(f2) see Definitions 

£(fi), see P[l6l3 

5(R n ), S'(R”) see Definitions [mol \V7A8 
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CHAPTER 1 



Introduction to problems of 
mathematical physics 



1. Temperature at a point? No! In volumes contracting to 

the point 

Temperature. We know this word from our childhood. The tem- 
perature can be measured by a thermometer. . . This first impression 
concerning the temperature is, in a sense, nearer to the essence than 
the representation of the temperature as a function of a point in 
space and time. Why? Because the concept of the temperature as 
a function of a point arose as an abstraction in connection with the 
conception of continuous medium. Actually, a physical parameter 
of the medium under consideration (for instance, its temperature) 
is first measured by a device in a “large” domain containing the 
fixed point £, then using a device with better resolution in a smaller 
domain (containing the same point) an so on. As a result, we ob- 
tain a (finite) sequence of numbers {aq, . . . , «m} the values of the 
physical parameter in the sequence of embedded domains contain- 
ing the point £. We idealize the medium considered, by assuming 
that the construction of the numerical sequence given above is for 
an infinite system of domains containing the point £ and embedded 
in each other. Then we obtain an infinite numerical sequence {a m }. 
If we admit (this is the essence of the conception of the continuous 
medium 1 ^, that such a sequence exists and has a limit (which does 
not depend on the choice of the system of embedded in each other 
domains), then this limit is considered as the value of the physical 
parameter (for instance, temperature) of the considered medium at 
the point £. 

11 In some problems of mathematical physics, first of all in nonlin- 
ear ones, it is reasonable (see, for example, HUTTUMI3IS1 HU) to 
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consider a more general conception of the continuous medium in which a 
physical parameter (say, temperature, density, velocity. . . ) is character- 
ized not by the values measured by one or another set of “devices”, in 
other words, not by a functional of these “devices”, but a “convergent” 
sequence of such functionals which define, similarly to nonstandard anal- 
ysis [5] 1 13 . 73] . a thin structure of a neighbourhood of one or another 
point of continuous medium. 

Thus, the concept of continuous medium occupying a domain 2 ) 
id, assumes that the numerical characteristic / of a physical param- 
eter considered in this domain (i.e. in Q) is a function in the usual 
sense: a mapping from the domain f 1 into the numerical line (i.e. into 
R or into C). Moreover, the function / has the following property: 

(f^ii) = a m, m = 1 , . . . , M. (1.1) 

Here, a m are the numbers introduced above, and the left-hand side 
of which is defined 3 ) by the formula = J f(x)(p^(x)dx, 

represents the “average” value of the function /, measured in the 
neighbourhood of the point £ € by using a “device”, which will 
be denoted by (•, p^). The “device” has the resolving power, that 
is determined by its “device function” (or we may also say “test 
function”) tffc ; Q, — > R. This function is normed: f p^(x)dx = 1. 

2) Always below, if the contrary is not said explicitly, the domain Q 
is an open connected set in R", where n > 1, with a sufficiently smooth 
(n — l)-dimensional boundary <9f2. 

3 ) Integration of a function g over a fixed (in this context) domain 
will be often written without indication of the domain of integration, and 
sometimes simply in the form J g. 

Let us note that more physical are “devices”, in which has 
the form of a “cap” in the neighbourhood of the point £, i.e. p^(x) = 
p{x — ^) for x £ fi, where the function ip : R n 9 x = (x \, . . . , x n ) \ — > 
ip(x) £ R has the following properties: 



‘P > 0 , 



/ 



<p = l, 



<p = 0 outside the ball {t6R n |cc| = 




+ • • • + xl < p}. 



(1.2) 



Here, p < 1 is such that {x £ H | \x — £| < p} C H. Often one 
can assume that the “device” measures the quantity / uniformly 
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in the domain to £ Q. In this case, tp = l w /|a;|, where is the 
characteristic function of the domain to (i.e. 1^ = 1 in to and 1^ = 0 
outside to), and |w| is the volume of the domain to (i.e. |w| = f l u ). In 
particular, if = R” and to = {a; £ R"| |ar| < a}, then <p(x) = S a (x), 
where 

c / n S«~ n /\Bn\ for |x| < a , 

I (J jor |x| > a, 

and \Br,\ is the volume of the unit ball B n in R n . 



(1.3) 



1.1. P. It is well known that |i?2| = 7 r and | -B3 1 = 47 t/3 . Try to 
calculate \B n \ for n > 3. We shall need it below. 

Hint. Obviously, \B n \ = cr n /n, where a n is the area of the 
surface of the unit (n — l)-dimensional sphere in R", since \B n \ = 
Jo rn ~ 1(J ndr. If the calculation of a n for n > 3 seems to the reader 
difficult or noninteresting, he can read the following short and unex- 
pectedly beautiful solution. 

Solution. We have 



/ oxj \ n 

(/ e "‘ I<s ) =/ 

Wno / UJ ri 



= / e~^dx = 



e r r n 1 a n dr = (a n /2) ■ T(n/2), 



where T(-) is the Euler function defined by the formula 



(1.4) 



OO 

T(A) = J t x ~ 1 e~ t dt , where SKA > 0. (1.5) 



For n = 2 the right-hand side of (1.4 1 is equal to 7r. Therefore, 



e 4 dt= Wn. 



( 1 . 6 ) 



Thus, cr ra = 27T ri / 2 r -1 (n/2). By taking n = 3, we obtain 2T(3/2) = 
\fi r. By virtue of the remarkable formula T(A + 1) = A • T(A), (ob- 
tained from ( 1.5 1 by integration by parts and implying r(n+l) = n!), 
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this implies that T(l/2) = Now we get that 



02 n — 



2 ? t h 

(n — 1)! ’ 



27T n 

CT2n+1 “ (74-1/2) -(74-3/2) 



3/2 • 1/2 ' 
(1.7) 
□ 



2. The notion of (5-sequence and 5-function 

In the preceding section the idea was indicated that the defini- 
tion of a function / : Q — > R (or of a function / : Q — » C) as a 
mapping from a domain C R n into R (or into C) is equivalent to 
determination of its “average values” : 

= J f(x)ip(x)dx, (2.1) 

n 

where d> is a sufficiently “rich” set of functions on Q. A sufficiently 
general result concerning this fact is given in Section [10] Here, we 
prove a simple but useful lemma. Preliminary, we introduce for 
e e]0, 1] the function 5 e : R ra — > R by the formula 

S E (x) = ^p{x/e)/£ n , v? > 0, J (p = 1, <p — 0 outside B n . (2.2) 
Let us note that for l/g > 1 we have 



J 5 e (x) dx = J S e ( x — £) dx = 1, £ £ fl. 

R" n 

2.1. Lemma. Let f £ C(Q), i.e. f is continuous in fl C R n . 
Then 

/(£) = !im [ f(x)S e (x - £) dx, (eff, (2.3) 

n 

i.e. the function f can be recovered by the family of ‘‘average values ” 



f(x) ■ 5 e (x 



l ;)dx 



jen, e>o 



PROOF. For any rj > 0, there exists e > 0 such that |/(a;) — 
/(£)| < 7? if |x — £| < e. Therefore, 



2. THE NOTION OF ^-SEQUENCE AND ^-FUNCTION 



5 




f f(x) 5 e(x-£)dx 


-m 


= 


[ (f(x) - f(Q)fie(x - £) dx 


\_J 

n 






J 

n 



< J \f(x) - f(€)\6 e (x-t;)dx 

l*-£l<e 




2.2. Definition. Let $ be a subspace of the space C'(fl) and 
f £ fh A sequence {S e (x — £),x £ fi} e6 R ]£ ^o of functions x i — > 
5 e (x — £) such that equality (2.3) holds for any / £ C(f2) (for any 
/ £ <&) is called 5-sequence (on the spaced d>) concentrated near the 
point £. The last words are usually skipped. 



b The notion of 5-sequence on the space <f> allow to obtain a series of 
rather important results. Some of them are mentioned at the beginning 
of Section [4] 

In Section [4] some examples of 5-sequences on one or another 
subspace C C(Jl) are given. Important examples of such sequences 
are given in Section [3] 
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2.3. Definition. A linear functional 2 -* 5 e defined on the space 
C( Cl) by the formula 

: C7(fi) 9 / ^ /(0 e K ( orC), (2.4) 

is called the 5-function , or the Dirac function concentrated at the 
point £. 



2 * A (linear) functional on a (linear) space of functions is defined as 
a (linear) mapping from this functional space into a number set. 



Often one writes (5-function (2.4) in the form Ex- 
action on a function / £ C(O) writes (see formula @0 



- £) and its 
in the form 



{f(x),S(x-0) = f(0 ° r (H x ~ £)>/ 0 )) = /(£)■ ( 2 - 5 ) 



The following notation is also used: (/, 6(f) = /(£) or (5^, f) = 
/(£). The Dirac function can be interpreted as a measuring instru- 
ment at a point (a “thermometer” measuring the “temperature” at 
a point). If £ = 0, then we write <5 or 5(x). 



3. Some spaces of smooth functions. Partition of unity 

The spaces of smooth functions being introduced in this section 
play very important role in the analysis. In particular, they give 
examples of the space $ in the “averaging” formula (2.1 1 . 

3.1. Definition. Let Cl be an open set in R", Cl the closure of 
Cl in R", and m £ Z+, i.e., m is a non-negative integer. Then 1 -* 

11 If to = 0, then the index m in the designation of the spaces defined 
below is usually omitted. 



13.11 1. C m {Cl) (respectively, C™(Cl)) is the space! C m (Cl) of m- 
times continuously differentiable (respectively, with bounded 
derivatives) functions tp : Cl — > C, i.e., such that the func- 
tion d a ip is continuous (and respectively, bounded) in Cl for 
| a | < m. Here and below 

c?l“l u>(x) 

d a ip(x) = Q x a! Qx an ’ = ai + ' ’ oi j £ Z + = {0, 1, 2, . . . }. 

The vector a = (or, . . . , a n ) is called a multiindex. 
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13.11 2. C m (f2) = C m (M")| n , i.e. 2 ) C m (fl) is the restriction of the 
space C' m (R n ) to ft. This means that <p £ C m (Cl) -£=> 
there exists a function ip £ C m (R") such that p(x) = ip{x) 
for x £ fl. 

2 ' ) The space C m (fl), in general, does not coincide with 
the space of functions m-times continuously differentiable up to 
the boundary. However, they coincide, if the boundary of the 
domain is sufficiently smooth. 

13.11 3. PC m (fl) (respectively, PC™( fi)) is the space of functions 
m-times piecewise continuously differentiable ( and respec- 
tively, bounded) in Q: this means that p £ PC m (fA) (re- 
spectively, p £ PC™(fl)) if and only if the following two 
conditions are satisfied. First, p £ C m (fl \ K 0 ) (respec- 
tively, ip £ C™( \ K 0 )) for a compact 3 ) K 0 C fi. Second, 
for any compact I\ C Q there exists a finite number of do- 
mains Clj C Cl, j = 1, . . . , N, each of them is an intersection 
of a finite number of domains with smooth boundaries, such 
that K C Ujlifi, and p\^ £ C m (tu) for any connected 

( ( N \ f N 

component oj of the set n U n jJ\ 

3 ) A set K C R n is called compact, if K is bounded and 
closed. 

13.11 4. The support of a function ip £ C(fl), denoted by supp^, 
is the complement in f l of the set {x £ Q \ p(x) ^ 0}. 
In other words, Slippy is the smallest set closed in f 1 such 
that the function p vanishes outside this set. 

13.11 5. C™(0,) = {p £ C m (fl ) | suppip is a compact}. 

13.11 6. = {<£ £ C™(fl) j suppip C fl}. 

[3Ti 7. C7°°(n) = n c m (n),. . . , c$°(n) = n cg*(n). 

m m 

08. If p £ C™{ Cl) (or ip e Cq°(CI)) and supp<p C lo, where w 
is a subdomain of fi, then the function ip is identified with 
its restriction to u. In this case we write: ip £ C™(uj) (or 
p £ 

3.2. Definition. We say that a set A is compactly embedded in 
fi, if A is a compact and A C tt. In this case we write A <£ Cl. 
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Obviously, C™( fl) = {ip G | supp<p <e fl}, and 

cy(f2) c cj*(f2) c c m (fi) c c m (o) c pc m (n), 

where the first inclusion and the third one should be replaced by =, 

if n = I™. 



3.3. Example. 

C 0 m (l n ) 9 : x i — ♦ <p(x) 



(1 — |x| 2 ) m+1 for |x| < 1, 
0 for |x| > 1. 



3.4. Example. 

C^(R n ) 9 <p : x i — ► <p{x) 



exp(l/(|x| 2 - 1)) for |x| < 1, 
0 for |x| > 1. 



3.5. Example (A special case of ( |2.2[ )). Let e > 0. we set 

<Pe(x) = <p(x/e), (3.1) 

where ip is the function from Example |3.4| Then the function 

S e : x i — * <p e (x)/ J ip e (x)dx , x G R n , (3.2) 



belongs to C'g°(K ra ) and, moreover, 

5 e (x) >0, Vx G R", S e (x) = 0 for |x| > e, 



S e (x) dx = 1. (3.3) 



3.6. Example. Let x G 1, ip : x i — » y>(x) = f _ 0 0 g{j)dr, where 
(see Fig.) g(—x) = —g{x) and g(x) = S e (x + 1 + e) for x < 0 (S e 
satisfies (3.3 1). 
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We have tp £ C^°(M), 0 < ip < 1, ip(x) = 1 for |x| < 1. 

3.7. Example. Let (xi, . . . ,x n ) £ M” be Euclidian coordinates 
of a point x £ R n . Taking ip from Example 3.6 we set 

ijjvix) = ^2 <p(xi+2ki) <p(x n +2h n ), kj £ Z, |fc| = fcH \-k n . 



\k\=v 



/ oo 



Then the family {<^„}£Z 0 of the functions <p v (x) = if v (x)/ ( ^v{x) 

\ i >=0 

form a partition of unity in Cl = K n , i.e. , p v £ C'o°(f2), and 

(1) for any compact I\ C 12, only a finite number of functions 
ip u is non-zero in K; 

(2) 0 < p „(x) < 1 and <Pv(x) = 1 Vx £ 12. 

V 

3.8. Proposition. For any domain w g fl, there exists a func- 
tion ip £ Cq°(C2) such that 0 < ip < 1 and tp(x) = 1 for x £ u>. 



PROOF. Let e > 0 is such that 3e is less than the distance from 
u) to <912 = Cl\Cl. We denote the e-neiglrbourhood of u by w e . Then 
the function 



= J 6 e (x - y)dy, 



x £ 12 



(<5 e from (3.2 I ) has the required properties. 



□ 



3.9. Definition. Let {12„} be a family of subdomains 12„ <g 12 
of a domain 12 = Uf2„. Suppose that any compact K (e 12 has 
a nonempty intersection with only a finite number of domains 12„. 
Then we say that the family {12„} forms a locally finite cover of Cl. 



3.10. Theorem (on partition of unity). Let {12„} be a locally 
finite cover of a domain Cl. Then there exists a partition of unity 
subordinate to a locally finite cover, i.e., there exists a family of 
functions p v £ Cfi°(Cl„) that satisfies conditions (l)-(2) above. 
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The reader can himself readily obtain the proof; see, for example, 
The partition of unity is a very common and convenient tool, by 
using which some problems for the whole domain fi can be reduced 
to problems for subdomains covering f 1 (see, in particular, in this 



connection Sections 11 20 and 22). 



4. Examples of 5-sequences 

The examples in this section are given in the form of exercises. 
Exercise P 14.ll will be used below in the deduction of the Poisson 
formula for the solution of the Laplace equation (see Section [5j, 
P 14.2I will be used for the Poisson formula for the solution of the 
heat equation (see Section [6j, P 4.3 will be used in the proof of the 



theorem on the inversion of the Fourier transform (see Section 17 1, 
and with the help of P |4.3| the Weierstrass theorem on approximation 
of continuous functions by polynomials can be easily proved (see 



Section 19 1 . 



4.I.P. Show that the sequence }j, >-q-o of the functions S y (x) = 

\ x 2 + v i ! where x £ R, is a 8-sequence on the space C&QR) (see Defi- 
nition 3.1\1 ) and is not a 8 -sequence on C(R). 



4.2.P. Show that the sequence {5 t } t _, + o of the functions 8t(x) = 






e x / 4t where x 



is not a 8 -sequence on C'(K), but is a 



8-sequence on the space 4) C C(R) of the functions that satisfy the 
condition \/ip £ 4> 3a > 0 such that \ip{x) exp(— ax 2 )\ — > 0 for |cc| — > 
oo. 



4.3.P. Show that the sequence {8 u }i/„^ 0 of the functions 8 u (x) = 
'r-, where x £ R, is a 8 -sequence on the space 4) C C 1 (K) of the 



functions ip such that 

/ |y>(x)| dx < oo, 



|y/(x)| dx < oo, 



/ 2 \ fc 

4.4.P. By taking the polynomials 8k{x) = f 1 — , where 

x £ K. and k is a positive integer, show that the sequence 
is a 8 -sequence on the space Co(M), but is not a 8 -sequence on the 
space C'h(K) (compare P f.l) 
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4.5. Remark. For exercises P |4.1| -P |4.4| it is helpful to draw 
sketches of graphs of the appropriate functions. Exercises P |4.1| - 
P |4.2| are simple enough, exercises P |4.3| P |4.4| are more difficult, 
because the corresponding functions are alternating. In Section |13| 
Lemma |13.10| is proven that allows us to solve readily PlOI-PMl 
While solving P |4.2| -P |4.4| one should use the well-known equalities: 

OO OO 

[ e~ v dy = y/n, [ ^ X dx = 7 r, lim (1 — alvY = e~ a . 

J J X 

— OO — OO 

5. On the Laplace equation 

Three pearls of mathematical physics. Rephrasing the title of 
the well-known book by A.Ya. Khinchin [33], one can say so about 
three classical equations in partial derivatives: the Laplace equation, 
the heat equation and the string equation. One of this pearls has 
been found by Laplace, when he analyzed 1 - 1 Newton’s gravitation 
law. 

b See in this connection Section 1 of the book by S.K. Godunov m- 



5.1. Definition. A function u £ C 2 (fl) is called harmonic on 
an open set f l C M™, if it satisfies in fl the (homogeneous 2 )) Laplace 
equation 

r\ 2 o2 

Am = 0, where A : C 2 (Ll) 9 u i — * Am = — ^ + • • • + y £ C(Q), 

(jX 2 CJ x 

( 5 . 1 ) 

and xi, ... ,x n are the Euclidian coordinates of the point x £ fl C 
R". Operator 3 ^ (H3 denoted by the Greek letter A - “delta” — is 
called the Laplace operator or Laplacian. 

The equation Am = / with non-zero right-hand side is sometimes 
called the Poisson equation. 

3) By an operator we mean a mapping / : A' — > Y, where A' and Y 
are functional spaces. 
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5.2. P. Let a function u £ C 2 (Ll), where x £ ft C R", depend 
only on p = \x\, i.e., u(x) = v(p). Show that A u also depends only 



on p and Au = d g^ + 



n—L dv(p) 

P dp ' 



Harmonic functions of two independent variables are closely con- 
nected with analytic functions of one complex variable, i.e., with the 
functions 



w(z) = u(x, y) + iv{x, y), z = x + iy £ C, 
which satisfy the so-called Cauchy -Riemann equations 



- v y = 0 , 



= 0 . 



(5.2) 



Here, the subscript denotes the derivative with respect to the corre- 
sponding variable, i.e., u x = du/dx, . . . , u yv = d 2 u/dy 2 , . . .). From 
(5.2) it follows that 



T Uyy {U X 'dy)x T (V X T Uy)y 0 



u vv 



= 0 . 



Thus, the real and imaginary parts of any analytic function w{z) = 
u(x,y) + iv{x,y) are harmonic functions. 

Consider the following problem for harmonic functions in the 
half-plane = {(a :,y) £ K 2 | y > 0}. Let / £ PCf,(K), i.e., (see 
Definition 3. 1|3| ) / is a bounded and piecewise continuous function 
of the variable x £ R. We seek a function u £ C 2 
Laplace equation 



u vv ~ 0 



m 



l 2 + 



+ ) satisfying the 
(5.3) 



and the following boundary conditions 

lim u(x,y) = f{x), 



v — 



(5.4) 



where a; is a point of continuity of the function /. 

Problem (5.3 1— ( 5.4 1 is called the Dirichlet problem for the Laplace 
equation in the half-plane and the function u is called its solution. 
The Dirichlet problem has many physical interpretations, one of 

which is given in Remark EES 

According to Exercises P5.12 and P5.15 below, problem (5.3) 



(5.4) has at most one bounded solution. We are going to find it. 
Note that the imaginary part of the analytic function 



ln(a; + iy) = In \x + iy \ + i arg(a; + iy), (x, y) £ 
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coincides with arccot(x/y) g]0, 7r[. Hence, this function is harmonic 
in Moreover, 



x 

hm arccot — 

2/ — >+0 y 



7 r, if x < 0, 
0, if x > 0. 



These properties of the function arccot {x/y) allow us to use it in 
order to construct functions harmonic in with piecewise constant 
boundary values. In particular, the function 



Pe{x,y ) 



1 I" x — e x + e 

arccot arccot 

2ne [ y y 



is harmonic in and satisfies the boundary condition 



lim P e (x,y) = 6 e (x) for |x| ^ e, 
y —>+ o 



where the function 5 e (x) is defined in ( 1.3 1 . On the other hand, if x 



is a point of continuity for /, then, by virtue of Lemma 2.1 



OO 

f{x) = lim j S e (f - x)f(f)df. 

— OO 



This allows us to suppose that the function 



OO 

M 2 B {x,y) i — ■> lim J P e (f - x,y)f(f)d£, (5.5) 

— OO 



assumes (at the points of continuity of /) the values f(x) as 
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y — > +0 and that this function is harmonic in R 2 , since 

^dx^ + dy^) - x ’ y)f(€k)(€k+i - 6 )^ = 



N 



^2f{£k){£,k+i - 6 )) 



*= i 



d 2 ^2 

Ox 2 + 0y2 



= 0 . 



The formal transition to the limit in (5.5 1 leads us to the Poisson 
integral (the Poisson formula) 



i(x,y) = [ f(QP(x - £,y)d£, where P{x 7 y) = - V , (5.6) 
J 7 rr + y- 



since 



^ \ 1 d ( x\ 1 y 

lim P e (x, y) = w- arccot - = -- 5 - — 

£ — >o 7 t ax \ y J 7 r x* + y- 



Let us note that condition (5.4 1 is satisfied by virtue of P 4.1 More- 



over, note that the function u is bounded. Indeed, 

oo oo 

\u(x,y)\= J \m\P{x-Z,y)dt<C J P{x - £,y) dy = C. 



We are going to show that the function u is harmonic in R 2 . Differ- 
entiating (|5.6|, we obtain 



& +k u{x,y) = f &+ k 
dxJdy k J dxWy 



k P(x-£,y)d£ Vj > 0, Vfc > 0. 

(5.7) 



Differentiation under the integral sign is possible, because 

d^+ k P(x-ty) 



dxWy k 



< 



1 + ^| 2 for \^\<R, ^ < y < R, (5.8) 



where C depends only on j > 0, k > 0 and R > 1. From (5.7 1 it 
follows that 

OO 

A u(x,y)= J f(€)AP(x-£,y)d£. 
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i+> 



since 



However, A P(x — £, y) = A P(x, y) = 0 in 

. 1 d ( x\ ( x\ d d 

Fix, y) = -7— arccot — , A arccot — ] = 0 and A— = —A. 

tt dx \ yj \ 



V 



dx dx 



Thus, we have proved that the Poisson integral (5.6 I gives a solution 



of problem (5.3 1— (5.4 1 bounded in R 2 . 



5.3.P. Prove estimate ( 5.8 1 . 



5.4. Remark. The function P defined in (5.6 1 is called the Pois- 



son kernel. It can be interpreted as the solution of the problem 
A P = 0 in R+, P(x,0) = S(x), where S(x) is the (5-function 4 ). 



4 ' The formula ( 5.6 1 which gives the solution of the problem A u = 0 



in R+, u(x, 0) = f(x), can be very intuitively interpreted in the following 
way. The source “stimulating” the physical field u(x, y) is the function 
f(x) which is the “sum” over £ of point sources f(£)5(x — £). Since one 
point source 5(x — t;) generates the field P(x — £, y), the “sum” of such 
sources generates (by virtue of linearity of the problem) the field which is 
the “sum” (i.e., the integral) by £ of fields of the form f(£)P(x — £, y). In 
this case, physicists usually say that we have superposition (covering) of 
fields generated by point sources. This superposition principle is observed 
in many formulae which give solutions of linear problems of mathematical 



physics (see in this connection formulae ( 5.10 1, ( 6.15 1, ( 7.14 1,. . . ). In these 



cases, mathematicians usually use the term “convolution” (see Section 19 1 . 



5.5. Remark. Note that the function P is (in the sense of the 
definition given above) a solution unbounded in of problem 



(5.3) (5.4), if f{x) = 0 for x / 0 and /( 0) is equal to, for in- 
stance, one. On the other hand, for this (piecewise continuous) 
boundary function /, problem (5.3 1 (5.4 1 admits a bounded solu- 



tion u(x,y) = 0. Thus, there is no uniqueness of the solution of the 



Dirichlet problem (5.3) (5.4[) in the functional class C 2 (MU- I' 1 this 



connection also see P 5.6 



P 15712 and P 5.15 



5.6. P. Find an unbounded solution u £ C ,00 (RU of problem 



( 5.3 ) — ( 5.4 ) for /( s) = 0. 



5.7. P. Let k £ R, and u be the solution of problem (5.3 1 -(5.4) 



represented by formula (5.6 1 . Find lim (xq + ky,y) in the two cases: 
y — >+o 
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(1) / is continuous ; 

(2) f has a jump of the first kind at the point Xo ■ 

5.8. P. Prove 

5.9. Proposition. Let it and oj he two domains inM 2 . Suppose 
that u : LI 9 ( x , y) i — » u{x , y) G K. is a harmonic function and 

z(0 = if) + iy(£, rj), (£, if) e u C K 2 

is an analytic function of the complex variable f = £ + irj with the 
values in LI (i.e., {x,y) G LI). Then the function 

U{f,if) = u(x{f,r]),y{f,if), (f,r)) G w 

is harmonic in u>. 



5.10. P. Let p and ip be polar coordinates in the disk D = {p < 
R, ip G [0, 27t[} of radius R. Suppose that f G PC(dD), i.e., f 
is a function defined on the boundary dD of the disk D and f is 
continuous everywhere on dD except at a finite number of points, 
where it has discontinuities of the first kind. Consider the Dirichlet 
problem for the Laplace equation in the disk D: to find a function 
u G C 2 (D) such that 



Au = 0 in D, lim u(p, <p) = f(Rp>), (5.9) 

p—>R 

where s = Rip is a point of continuity of the function f G PCb(dD). 
Show that the formula 



2t tR 



u{p , p) = 



2irR 



m 



( R 2 — p 2 )ds 
R 2 + p 2 — 2Rpcos(<p — 9) 



s 

R 



(5.10) 



represents a bounded solution of problem ( 5.9 1 . Formula (5.10) was 
obtained by Poisson in 1823. 



Hint. Make the transformation w = Rfrf of the half-plane 



onto the disk D and use formula (5.6 1 . 



5.11.P. Interpret the kernel of the Poisson integral (5.101, i.e., 
the function 

1 (R 2 - p 2 ) 



2irR R 2 + p 2 — 2Rpcos(tp — 9) ’ 

similar to what has been done in Remark \5.f \ with respect to the 
function P . 
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5.12.P. Using Theorem \5.l3\ below, prove the uniqueness of the 
solution of problem (5.9 1 as well as the uniqueness of the bounded 
solution of problem (5.3 1 (5.4) in the assumption that the bounded 
function f is continuous. (Compare with Remark 5.5 ) 



5.13. Theorem (Maximum principle). Let LI be a bounded open 
set in R" with the boundary <912. Suppose that u £ (7(12) and u is 
harmonic in 12. Then u attains its maximum on the boundary of the 
domain 12, i.e., there exists a point x° = {x \, . . . , x°f) £ <912 such that 
u(x) < u(x°) \/x £ 12. 



Proof. Let m = sup u(x), M = supu(;r) = u(x°), x° £ 12. 
x€dQ 

Suppose on the contrary that m < M. Then x° £ 12. We set 



v(x) = u(x) + 



M — m , 
2d 2 



\x — x 



oi2 



where d is the diameter of the domain 12. The inequality \x — x°\ 2 < 
d 2 implies that 



j(x) < m + 



M — m M + m 
d 2 = 



< M, x £ 912. 



2d 2 " 2 

Note that v(x°) = u(x°) = M. Thus, v attains its maximum 
at a point lying inside 12. It is know that at such a point Av < 0. 
Meanwhile, 



. . M — m. ( v-^/ 0 , 2 \ M — m „ 

Av = A u H — A I 2_^{ x k — x k ) ) = — — • 2n > 0. 



2d 2 



2d 2 



The contradiction obtained proves the theorem. 



□ 



5.14. P. In the assumptions of Theorem 5.13\ show that u attains 
its minimum as well on <912. (This is the reason why the appropriate 
result (Theorem 5.13) is known as minimum principle.) 

5.15. P. Using Theorem \5.16 below , prove the uniqueness of the 
bounded solution of problem (5.9 1 as well as problem (5.3 1 (5.4 1 . 
Compare with Exercise P\5.1£\ 



5.16. Theorem (on discontinuous majorant). Let 12 be a 
bounded open set in R 2 with the boundary <912, and F a finite set 
of points Xk £ 12, k = 1 , ,N. Let u and v be two functions har- 
monic in 12 \ F and continuous in 12 \ F. Suppose that there exists 
a constant M such that |7t(a;)| < M, |t'(a;)| < M \/x £ 12 \ F. If 
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u(x ) < v(x) for any point x £ dLl \ F, then u(x) < v{x) for all 
points x £ Ll \ F. 



PROOF. First, note that the function In |x|, where x £ R 2 \ {0}, 
is harmonic. We set 



w £ (x) = u(x) 



N 

- v(x) - Y 



fe = 1 



2M d 

In (d/e) \x-x k [ 



Here, 0 < e < d, where d is the diameter of Ll, hence, ln(d/|x — x^l) > 
0. Consider the domain Ll € obtained by cutting off the disks of 
the radius e centered at the points x k £ F, k = 1, . . . , N, from Ll. 
Obviously, w e is harmonic in Ll £ , continuous in f l e , and w e (a;) < 0 for 
x £ dL l £ = \ Therefore, by virtue of the maximum principle, 

We(x) < 0 for x £ fl e . It remains to tend e to zero. □ 



5.17. P. Let OcK 2 he a simply connected open set hounded by 
a closed Jordan curve dCl. Suppose that a function f is specified on 
dfl and is continuous everywhere except a finite number of points at 
which it has discontinuities of the first kind. Using Proposition \ 5. 9| 
the Riemann theorem on existence of a conformal mapping from LI 
onto the unit disk (see, for instance, [38j ). prove that there exists a 
bounded solution u £ C 2 {Li) of the following Dirichlet problem: 

A u = 0 in LI, lim u(x) = f(s), (5.11) 



where s is a point of continuity of the function f £ PC(dLl). Us- 
ing Theorem 5. 1 6| prove the uniqueness of the bounded solution of 
problem (5.111 and continuous (make clear in what sense) depen- 
dence of the solution on the boundary function f. (Compare with 



Corollary 22.31 ) 



5.18. Theorem (on the mean value). Let u be a function har- 
monic in a disk D of radius R. Suppose that u £ C(D). Then the 
value of u at the centre of the disk D is equal to the mean u(x), 



x £ dD, i.e., (in notation of P 5.10) 



u \ a = 

l /»=0 



1 

2ttR 



2 nR 






(5.12) 



PROOF. The assertion obviously follows from (5.10). 



□ 
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5.19.P. Let u be a continuous function in a domain Ll C R 2 and 
u satisfy (5.121 for any disk D C Ll. Prove (by contradiction) that 
i/« / const, then u( x) < ||u|| \/x £ f l, where ||u|| is the maximum 
of |u| in Ll. 

By virtue of Theorem |5.18[ the result of Exercise |5.19| can be 
formulated in the form of the following assertion. 



5.20. Theorem (strong maximum principle). Let u be a har- 
monic function in a domain fl C R 2 . If u ^ const, then u(x ) < 
||it|| Vcc £ Ll, where ||u|| is the maximum of |u| in Cl. 

5.21. P. In the assumptions of P \5.li\ show that u is harmonic 
in the domain f 1. 



Hint. Let a £ Ll and D c ft be a disk centred at a. Suppose 
that v is a function bounded and harmonic in D such that v = u on 
dD. Using the result of P |5.19| show that the function w = u — v is 
a constant in D. 



5.22. Remark. It follows from P I5.21l and Theorem 15. 181 that a 
function continuous in a domain LI C R 2 is harmonic if and only if 
the mean value property (5.121 holds for any disk D C Ll. This fact 
as well as all others in this section is valid for any bounded domain 
LI C R n , n > 3, with a smooth [n — l)-dimensional boundary (with 
an appropriate change of the formulae). 

Let us note one more useful fact. 



5.23. Lemma (Giraud-Hopf-Oleinik). Let u be harmonic in LI 
and continuous in LI, where the domain LI <s R™ has a smooth (n— 1)- 
dimensional boundary T. Suppose that at a point x 0 £ T there exists 
a normal derivative du/dv, where v is the outward normal to T and 
u( x 0 ) > u(x ) Vx £ LI. Then du/dv\ x=x > 0. 

The proof can be found, for instance, in fl2l[23l [24l[28l [46l. 



6. On the heat equation 

It is known that in order to heat a body occupying a domain 
Hcl 3 from the temperature uq = const to the temperature u\ = 
const, we must transmit the energy equal to C ■ (u\ — no) ■ |0| to 
the body as the heat, where |H| is the volume of the domain Ll and 
C is a (positive) coefficient called the specific heat. Let u(x,t) be 
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the temperature at a point x = (xi,x 2 ,x 3 ) S O at an instant t. We 
deduce the differential equation which is satisfied by the function u. 
We assume that the physical model of the real process is such that 
functions considered in connection with this process (heat energy, 
temperature, heat flux) are sufficiently smooth. Then the variation 
of the heat energy in the parallelepiped 

II = {i £ R 3 | x% < x < x\ + hk, k = 1,2, 3} 

at time r (starting from the instant t°) can be represented in the 
form 



C ■ (u(a: 0 ,f 0 + t) — u(x°,t°)\ • |II| + o(r • |II|) 

= C ■ [u t (x°,t°) ■ T + o(t)] • |n| + o(r • |n|), (6.1) 



where |II| = h\ ■ h 2 • h 3 and o(A ) is small oofdsRasd- >0. 

This variation of the heat energy is connected with the presence 
of a heat flux through the boundary of the parallelepiped II. Ac- 
cording to the Fourier law, the heat flux per unit of time through 
an area element in direction of the normal to this element is propor- 
tional with a (negative) coefficient of proportionality to —k derivative 
of the temperature along this normal. The coefficient of proportion- 
ality k > 0 is called the coefficient of heat conductivity. Thus, the 
quantity of the energy entered into the parallelepiped II during the 
time t through the area element X\ = xj + hi is equal to 

f)u 

k{x\ + hi,x%,x%) • q^( x i + hi,xl,xl',t°) ■ t ■ h 2 ■ h 3 + o(t ■ |II|), 



and gone out during the same time through the area element aq = x\ 
is equal to 



k(x 



Xr,,X 



1; J '2> 



1 • —(x° • 
3) dxi [Xl " 



-2»*i;0 - r-h 2 -h 3 + o(t ■ |n|). 



Therefore, the variation of the heat energy in II caused by the heat 
flux along the axis Xi is equal to 



d 

dxi 



(k(x°) 



du 

dx\ 



(x°,t°))hi + o(h\) 



t ■ h 2 - h 3 + o(t ■ |II|). 



It is clear that the variation of the heat energy in II in all three 
directions is equal to the total variation of the heat energy in II, i.e., 
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(6.1). By dividing the equality obtained in this way by r • |II| 
tending r, h±, ft-2, and h 3 to zero, we obtain the heat equation 



and 



du <9 / du \ d / du \ d / du \ 

dt dx\ \ dxi ) 8x2 V ' dx2 ) + 8x3 \ 8x3 J 



( 6 . 2 ) 



If the coefficients C and k are constant, then equation (6.2 1 
rewritten in the form 



can be 



du _ fcPu cPu 8 2 u \ 

dt \ dx\ 8x\ 8x 3 ) ’ 



where a 



( 



> 0 . 



6.1. Remark. In the case when the distribution of the temper- 
ature does not depend on the time, i.e., Ut = 0, the temperature 
u satisfies the Laplace equation (if k = const). Thus, the Dirichlet 
problem for the Laplace equation (see Section [BJ can be interpreted 
as a problem on the distribution of stabilized (stationary) tempera- 
ture in the body, if the distribution of the temperature on the surface 
of the body is known. 



If we are interested in distribution of the temperature inside the 
body, where (during some time) influence of the boundary conditions 
is not very essential, then we idealize the situation and consider the 
following problem: 

du 

C— =div(fc • gradw), (i,!/,z)e8 3 , t> 0, 
u | t=0 =f(x,y,z), 

where / is the distribution of the temperature (in the body without 
boundary, i.e., in R 3 ) at the instant t = 0. This problem is sometimes 
called by the Cauchy problem for the heat equation. 

Suppose that / and k , hence, also u do not depend on y and 2. 
Then it is a solution of the problem 



du d / du\ 

c dt ~ dx v fc ‘ ’ 


(. x , t) € = {cc £ R, t > 0}, 


(6.3) 


u lt=o = 




(6.4) 
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6.2. Hint. The method which was used to solve problem (5.3 1 



(5.4) 1 ) suggests that the solution to problem (6.3 1 (6.4) can be rep- 



resented by the formula 



OO 

u(x,t)= j f(Qv(x-£,t)d£, 



where v is the solution of equation (6.3 1 , satisfying the condition 



t) = 6(x), where S is the 6-function. (6-5) 



b See note 4 in Section [ 5 ] 



Below (see Theorem |6.5| ) we show that this is true. 

Let us try to find the function v. It satisfies the following con- 
ditions: 



C 



dv 
d t 



±( k .*\ 
dx \ dx ) ’ 



Cv dx = Q , 



( 6 . 6 ) 



where Q is the total quantity of the heat that in our case is equal to 
C. Thus, we see that v is a function G of five independent variables 
x, t, C, k, and Q , i.e. , 



v = G(x, t , C , k, Q). 



(6.7) 



6.3. Remark. The method with the help of which we seek for 
the function v is originated from mechanics [[55 . It is known as the 
dimensionless parameters (variables) method. 



Note that the units of measurement of the quantities v, x, t , and 
Q in the SI system, for instance, are the following: [v] = K, [x] = 
m, [t] = sec, [Q\ = w. By virtue of ( |6.6| ), [G][u]/[t] = [fc][u]/[a:] 2 , 
[C\ [u] [.t] = [Q\ ; therefore the dimensions of the quantities C and k 
are expressed by the formulae: 

[C] = w/(m • K), [k] = w • m/(sec • K). 

Since C, k , and Q play the role of the parameters of the function 
v(x,t), it is preferable to express the units of measurement v and, 
say, x via [t], [C\, [ft] and [Q], We have 



N = Vlt) ■ « [v] = [Q]/VWV¥W\- 
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Let us take another system of units of measurement 

[t*] = vt[t], [<?*] = <j c [C]i [k*] = f T k [k), [Q*] = ctq[Q\, 

where at, ac, &k, <?Q are scaling coefficients, i.e., positive (dimen- 
sionless) numbers. We formulate the question: what are the values 
of the scaling coefficients a x and a v (for the variables x and v which 
are “derivatives” of the chosen “basic” physical variables t, C, k, and 
Q)7 We have 



x = 



MM = a x [x] = a.VMMM 

= cr x ^/ a c [t*}[k*]/ (a t a k [C*}) . 



Therefore, 

a x = y/a t ak/ac and similarly a v = aQ/^Ja t a k ac ■ (6.8) 

The numerical values t* , ,v* of the variables t, ... ,v in the new 
system of units are determined from the relation 



t*[t*]=t[t},...,v*[v*]=v[v]. 



Thus, 



t = 



<Jt 



C* = —, k* = Q* = 



<?c 

a~ta k 



JatOkac 

)— . 

a Q 



For instance, if [f] = sec and [f*] = hour, then at = 3600 and t* = 
t/3600. 

Let us now note that relation (6. 7 1 expresses a law that does not 
depend on the choice of the units. Therefore, 



v* = G(x*,t*,C*,k*,Q*) (6.9) 

with the same function G. Now, we choose the system of units such 
that t* = C* = k* = Q* = 1, i.e., we set a t = t, ac = C, a k = k, 
oq = Q. Then 



x* = xy / C/(kt), 



v* = vVtkC/Q. 
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Hence, by virtue of (6.9), we have 




where g(y) = G(y, 1,1, 1,1). 

( 6 . 10 ) 



6.4. Remark. We can come to formula (6.101 strongly mathe- 



matically. Namely, by making the change of variables 

t C k Q * x * v 

t* = — , C* = — , k* = — , Q = — , x* = — , v* = — , 

&t 

we require v* be equal to G(x * , t* ,C* ,k* ,Q*), i.e., we require that 

OO 

dv* d ( dv* \ f 
C 'W = 0^( k ’^)' / c "»Yi = (3*. 



Then (6.6 1 necessarily implies (6.8). Choosing, as before, the scaling 



coefficients a t , ac , a*,, and <tq, we again obtain (6.10). 

Nevertheless, it is helpful to use the dimension arguments. Firstly, 
they allow us to test the correctness of involving some parameters 
when formulating the problem: both sides of any equality used in the 
problem should have consistent dimensions. Secondly, the reasoning 
of dimension allows to find the necessary change of variables (not 
necessarily connected only with scaling coefficients). All these facts 
allow automatically (hence, easily) to get rid of “redundant” param- 
eters and so to simplify the analysis as well as the calculations 2 ). 
Moreover, the passage to dimensionless coordinates allows us to ap- 
ply the reasoning by similarity that sometimes essentially simplify 
the solution of rather difficult problems (see [55] ). 



2 * Consider the problem on the temperature field of an infinite plate of 
a thickness 25, with an initial temperature To = const, in the case, when 
there is the heat transfer on the surface of the plate (with the coefficient 
of the heat transfer a) with the medium whose temperature is equal to 
T\ = const. In other words, we consider the problem 

~dr ~ a <9£ 2 ’ T > 0’ l£l < S; =Ffc-7^r| 5=±s = a(T— Ti)| f=±s ; T| t=0 = To. 

The function T = /(r, £, a, S, k, a, Ti, To) depends a priori on eight pa- 
rameters. The tabulating the values of such a function, if each of the pa- 
rameters run over at least ten values, is unreasonable, because one should 
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analyze million pages. At the same time, the passage to the dimensionless 
parameters 

u = (T — Ti)/(Ti — To), x = £/S, t = ar/S 2 , <j = k/aS 



reduces this problem to the problem 



du 

dt 



d 2 u 

dx 2 



, t > 0, |a;| < 1; 



du 
u±o — 
ox 



= 0; u\ 



= 1 , 



( 6 . 11 ) 



whose solution u = u(t, x, o) can be represented (this is very important 
for applications) in the form of compact tables (one page for each value 
of a > 0). 



In order to find the function g and, hence, v we substitute ex- 
pression ( |6.10 I into the heat equation (6.2). We obtain 

Q^C/kt 3 [g(y)/2+yg\y)/2+g'\y)) = 0, i.e., (yg(y))' /2+g"(y) = 0. 

Thus, the function g satisfies the linear equation 

g(y) + yg{y)/ 2 = const. (6.12) 

If g is even, i.e., g(— y) = g(y), then g'(0) = 0; therefore, the func- 
tion g satisfies the homogeneous equation (6.12) whose solution, ob- 



viously, is represented by the formula g(y) = A ■ exp(— y 2 /4). The 
constant A is to be determined from the second condition in (6.6 1 : 



OO OO OO 

Q= J Cvdx = ACQ/VkCt j e~ Gx ^ ikt dx = 2AQ J e~^ d£, 



i.e., (accounting formula ( 1.6 1) A = l/(2y / 7r), hence, 



v(x,t) = ( Q/2VkC-Kt ) • exp(— Cx 2 /4:kt). 



(6.13) 



6.5. Theorem. Let f £ C(R), and for some o £ [1,2[, a > 0, 
and M > 0 the following inequality holds: 

|/0)| < M exp(a|*| CT ) VzeR. (6-14) 

Then the function u : R^_ = {0,t) £ R 2 \t > 0} — > R defined by the 
formula 

OO 

u(x,t) = J f(£)P(x - £,t)d£, P(x,t) = (l/2v / 7rf)exp(-a: 2 /4t), 



(6.15) 



— OO 
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is a solution of the heat equation 

in R^ = {0M)el 2 | i>0}. (6.16) 

This solution is infinitely differentiable and satisfies the initial con- 
dition (which is sometimes called the Cauchy condition) 

\hn o u ( x ,t) = f(x). (6.17) 

Moreover, VT > 0 3 C(T) > 0 such that 

\u(x,t)\ < C(T)exp((2a|a;|) IT ) Vs£l and VtG[0,T]. (6.18) 



PROOF. From construction of function ( 6. 13 1 , it follows that the 
function P(x,t) = (4-7rf) -1 / 2 exp(— x 2 /M) satisfies relations (6.6 1 in 
which C = k = Q = 1. Therefore, formula (6.171 follows from 



P 4.2 and (6.16) as well as the smoothness of the function u follows 



from the known theorem on differentiability of integrals with respect 
to the parameter (see, for instance, 72 )), because the appropriate 
integral converges uniformly, since Vi? > 1 Ve > 0 3 N > 1 such that 
f)j+k 

:(f(QP( x ~ £> t)) 



da ddt k 



df; < e 



(6.19) 



[1 /R,R]. For j + k > 0, the integrand in 



l£l >N 

for x G [— R, i?], t G 
(6.191 can be estimated for the specified x and t via Cnf(£)P(x — 
f,t). Therefore, in order to prove inequality (6.19) it is sufficient to 
establish estimate (6.18). Note that for a G [1, 2[, 

ier < 2 ff (kr + 1? - *r), ie - *r < - *o 2 + ^ - x\. 



Take e such that 1 — 4T • a 0 • e > 0, where a 0 = a ■ 2 CT . Then for 
t < T, (compare with [25] ) 

H x ,t) I < 7^ J e a|5r ’ exp(— (a: — f) 2 /At) 

< Mi e (2a|x|)<T [ e “( 2 l«-*ir . e -(C~0/zVt) 2 

J 2 Vt 



< M ie (2a|x|)<T J . e a 0 Cz(\x-i\/2Ci)2Vi. 

Setting r] = (£ — a;)(l — 4Ta 0 e) 1/,2 /(2v / t), we have 

\u(x,t)\ < C{T) e ( 2aM) ° J e-^+^^dr,. 



df_ 
2 yft 
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This implies estimate (6.18 1, if we note that 



OO OO OO 

J e -n 2 + av - 2 Vi dri = e c?t J e -( v -aVi) dri < g a 2 t J e -c 2 d £' 



□ 



6.6. Remark. In general, there exists a solution of problem 



(6.16l (6.17l different from (6.151. For instance, the function u(x,t) 



represented by the series 



u(x,t) = Y <P {r 

m—0 



\t) ■ x 2m /(2m)\, (. x,t ) G 






( 6 . 20 ) 



in which the function ip £ C° 



satisfies the conditions: 



supp^) C [0, 1], Vm £ Z + (t)\ < (7 m)!, where 1 < 7 < 2, 

( 6 . 21 ) 

is, obviously, a solution of problem ( 6. 16[ ) ( 6. 17 I for / = 0. (The 
condition 7 < 2 is needed for uniform (with respect to x and t, |x| < 



R < 00) convergence of series (6.201 ad its derivatives.) This simple 



but important fact was observed in 1935 by A.N. Tikhonov 



who, while constructing series (6.201, used the result by Carleman 



||9] on the existence of a non-zero function ip with properties (6.21 1. 
It is important to emphasize that the non-zero solution (6.20) of 



the heat equation constructed by Tikhonov (satisfying the condition 
u(x,0) = 0) grows for |x| — > 00 faster than expCx 2 VC > 0 (and 
slower than exp Cx a , where er = 2/(2 — 7) > 2). On the other hand, 
one can show, by using the maximum principle for the heat equation 
(see, for instance, [20l 1251 144 , 64] ) that the solution of problem 
(16. 16b (16. 17b is unique, if condition (16. 181 holds. The uniqueness 



theorem for a more large class of function was proved in 1924 by 
Holmgren m- 



From Remark 6.6 it follows 



6.7. Theorem. Let f £ C'(M), and f satisfy (6.14l. Then for- 
mula (6.15) represents a solution of problem (6.161 (6.17) and this 



solution is unique in class (6.18). 



28 1. INTRODUCTION TO PROBLEMS OF MATHEMATICAL PHYSICS 



7. The Ostrogradsky Gauss formula. The Green formulae 
and the Green function 



Let fl be a bounded domain in R n with a smooth [n — 1)- 
dimensional boundary <9f 1. Let / = (/i, . . . , /„) be a vector-function 
such that fk £ C(ft) and dfk/dxk € PC{ fl) Mk. It is known [ 63., 72 
that in this case the Ostrogradsky -Gauss formula ^ 




dfkjx) 

dx k 



dx 




■ a k dT 



(7.1) 



holds, where ak = a k {x) is the cosine of the angle between the 
outward normal v to T = <9fl at the point x £ T and /cth coordinate 
axis, and dr is the “area element” of T. For n = 1, formula (7.1 1 
becomes the Newton-Leibniz formula. 



^ Formula is a special case of the important Stokes theorem 

on integration of differential forms on manifolds with boundary (see, for 
instance, [63, 721 ). which can be represented by the Poincare formula: 
J n duj = f gn a ). The Poincare formula implies ( |7.1| ) for 

w = fk(x)dxi/\ . . . /\dx k -if\dx k+1 f\ . . . /\dxn , 

k 

because 

du> = ^2(df k (x)/dx k )dx, and ut| an = ^2 f k {x)a k dF. 



If fk{x) = A k { x)v(x), where v £ PC 2 (Q) fl then 

implies that 



(7.1) 




dx = — / A 



fc = i 



dv , 

k A dx+ 

OXk J 



Ak-va k dT . 



an 



fc = i 



(7.2) 



Setting A k = where u £ PC 2 (fl) fl G 1 (fl), we obtain the first 
Green formula 



J 



v ■ Audx = 



du 

v A~ dT 

ov 



on 



f du dv 

J “ dx k dx k 
n fc=1 



(7.3) 



where A is the Laplace operator (see Section [5]). Renaming it by v in 
( 7.3 ) and v by u and subtracting the formula obtained from ( 7.3 1 , we 
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obtain the so-called second Green formula for the Laplace operator 

(7.4) 



f f ( du 9v\ 

/ (v • A u — u ■ Av)dx = / I t; - — u ■ — J dT. 

n on 



Formula (7.4 1 implies (if we set v = 1) the remarkable corollary: 

du 



A u dx = 



dv 



dr. 



(7.5) 



n an 

In particular, it the function u G C rl (f2) is harmonic in Q, then 



J f^dr = 0. This is the so-called integral Gauss formula. 



on 



Rewrite formula (7.4) in the form 



J u(y)Av(y)dy = J v(y)Au(y) dy 
n n 



an 






dr. (7.6) 



Let us take a point x £ Q. Replace the function u in (7.6) by the 
function E a (x,-) £ PC 2 (G.) which depends on x as an parameter 
and which satisfies the equation 



n d 2 

A y E a (x,y) = ^2 -^~ 2 E ot{x,y) =S a (x-y), 
k = i ° Vk 



(7.7) 



where S a is defined in (1.3 1, and 1/a 1. Using the result of 



Exercise |7.1[ we tend a to zero. As a result, taking into account 
Lemma 12. 11 we obtain 



u{x) = J E(x - y)Au(y)dy 



an 



' dE{x-y) du(y ) 

u{y) E(x-y) 



dv 



dy. (7.8) 



7.1. P. Using the result of P\5.S\ and Theorem \ 5.13\ show that 
the general solution of equation (7.7 1 depending only on \x — y\ can 



he represented in the form E a {x — y) + const, where the function 
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E a £ C' 1 (K ra ) coincides for |x| > a with the function 
E(x) = 



(1/27 r) • In |x| 

-l/((n ^ 2 )<j r 



i/0, n = 2, 
n ~ 2 ) i/0, n > 3, 



(7.9) 



and for |x| < a, f/ie estimate |i? Q (x)| < |i?(x)| holds. Here, a„ 
denotes (see p\L IP i/ie area of the unit sphere in R n . 



Hint. By virtue of (7.5 1 and (1.3 1, 



(dE a /dv)d£ = 



AE n dx = 1. 



1 * 1=0 



\x\<c 



Let x £ fh We take 2 ) the function g(x,-) : Cl 9 y i — > g(x,y), 
which is the solution of the following Dirichlet problem for the ho- 
mogeneous Laplace equation with the special boundary condition 

& y g(x,y) = 0 in 11, g(x,y) = -E(x - y) for y G T. (7.10) 

Substituting the function g(x,-) into formula (|7.6[) for the function v 



and summing termwise the equality obtained with (7.8 1 , as a result, 
we have the following integral representation of the function u £ 
PC 2 (fd) (~l C 1 (f2): 



u(x) = J G(x,y)Au(y) 



dy- 



dG(x, y) 
dv 



u(y) dT, 



on 



where 



G(x, y) = E(x -y)+ g(x, y). 



(7.11) 



(7.12) 



In Section 22 the theorem is presented on existence of solutions for 
problems much more general than problem ( |7.10| ) . In the same section the 
theorem concerning the smoothness of the solutions is given. 



Function (7.121 is called the Green function of the Dirichlet prob- 



lem for the Laplace equation 

A u = f in n, 



u = tp on dfl. 



(7.13) 



This term is connected with the fact that, by virtue of (7.11|, the 
solution of problem (7.131, where / £ PC(fl), ip £ G(9H), can be 
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represented , using the function G , in the form 

dG(x, y) 



u(x)= f f(y)G(x,y)dy+ f ip(y)- 
n an 



dv 



dT. 



(7.14) 



Formula (7.14) is often called the Green formula. 



7.2. P. Let f l = R", where R™ = {x = (x',x n ) £ R” | x' £ 
R" _1 , x n > 0} and x' = (xi,. .. ,x n -i) £ R n_1 . Show that in this 
case G(x , y) = E(x , y) — E(x*, y), where x* = (a/, —x n ) is the flip of 
the point x over the hyperplane x n = 0. Verify (compare with (5.6)) 
that 

dG(x,y) | _ 2 x n 

dy Ij/»=o 



o-„ [(*i - j/i ) 2 H h (a; n _i - y n - 1 ) 2 + x n ] n/2 ' 



8. The Lebesgue integral 1 ^ 



The theory of the Lebesgue integral is, as a rule, included now in 
the educational programs for the students of low courses. Nevertheless, 
maybe some readers are not familiar with this subject. This section and 
the following one are addressed to these readers. At the first reading, one 
can, having a look at the definitions and formulations of the assertions of 
this sections, go further. For what follows, it is important to know at least 
two facts: 



(1) if a function / is piecewise continuous in P <s R n , then / is 
integrable in the sense of Riemann; 

(2) a function integrable in the sense of Riemann is integrable in 
the sense of Lebesgue, and its Riemann integral coincides with 
its Lebesgue integral (see P 8.15 below). 



If the need arises (when we consider the passage to the limit under the 
integral sign, the change of the order of integration and so on), it is ad- 
visable to come back to more attentive reading of Sections [5] [9] and the 
text-books cited. 



In Sections 01 we have outlined the idea of representability (in 
other words, of determination) of a function by its “average values”. 
This idea is connected with the notion of an integral. Let us recall 
that the definition of an integral, known to all the students since 
their first year of studies, has been given by Cauchy. It was the 
first analytic definition of the integral. Cauchy gave the first strict 
definition of continuity of a function. He proved that functions con- 
tinuous on a closed interval are integrable. In connection with the 
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development (Dirichlet, Riemann) of the concept of a function as a 
pointwise mapping into the real axis, the question arose concerning 
the class of functions for which the integral in the sense of Cauchy 
exists. This question has been answered by Riemann (see, for in- 
stance, [ 65 ]). This is the reason, why the integral introduced by 
Cauchy is called the Riemann integral. 

The space of functions integrable in the Riemann sense is rather 
large. However, it is not complete (see note 5 in Section [8]) with 
respect to the convergence defined by the Riemann integral similarly 
to the fact the set of rational numbers (in contrast to the set of real 
numbers) is not complete with respect to the convergence defined by 
the Euclidian distance on the line. Actually, let 

f n (x) = aC 1 / 2 for x &]l/n, 1] and f n {x) = 0 for a; G]0, 1/n]. 



Obviously, \f m ( x ) ~ fn{x)\dx — > 0 as m and n — * oo, i.e., the 
sequence {/&} is a fundamental sequence (see note 5 in Section [8]) 
with respect to the convergence defined by the riemann integral. It 
follows from the definition the Riemann integral that the function 
x i — > £ -1 / 2 (to which the sequence {f n (x)} converges pointwise) 
is not integrable in the Riemann sense. Moreover, one can easily 
show, using P |8.15| and Theorem [8T7J that there exists no function / 
integrable in the Riemann sense such that lim j \f(x) — f m (x) \ dx = 



0, i.e., the space of functions integrable in the Riemann sense, is not 
complete with respect to the convergence defined by the Riemann 
integral (see also Exercise P 8.23). 

This reason as well as some others stimulated (see, for instance, 
[ 65 ]) the development of the notion of integral. A particular role 
due to its significance is played by the Lebesgue integral. In 1901 
26-year Lebesgue introduced (see Definition 8.12 below) the space 
L(Q) of functions defined on an open set C M n and called now 
by integrable in the Lebesgue sense and the integral that is called 
now by his name (see Definition 8.121. This integral was defined by 
Lebesgue as the functional 



L(Sl) 9 / i 



/ £ K, 



which in the case f 1 =]a, b[ is noted by the standard symbol and 
possesses the following six properties: 
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b-\-h 



( 1 ) J f(x ) dx = J f(x — h) dx for any a, b, and h. 

(L Q,-\~h 

b c a 

(2) J f(x) dx + J /( x) dx + J f(x) dx = 0 for any a, b, and 

a b c 

C. 

b b b 

(3) J[f(x)+g(x)]dx = J f(x)dx + J g(x) dx for any a and b. 

a a a 

b 

(4) J f(x) dx > 0, if / > 0 and b > a. 

a 

1 

(5) J 1 • dx = 1. 

o 

(6) If, when n increases, the function f n {x) tends increasing 
to f(x), then the integral of f n (x) tends to the integral of 

f(x)- 



“Condition 6, - wrote Lebesgue p9 , - takes a special place. 
It have neither the character of simplicity as the first five nor the 
character of necessity” . Nevertheless, it is Condition 6 that became 
a corner-stone in Lebesgue’s presentation of his theory of integration. 

Below we give 2 ) the construction theory of the Lebesgue inte- 
gral and the space L(f2) and formulate results of the theory of the 
Lebesgue integral that we need in further considerations. 



Following mainly the book by G.E. Shilov and B.L. Gurevich ES; 
also see [Ml 1551 . 



8.1. Definition. A set A c D is said to be a set of zero measure, 
if Ve > 0 there exists a family of parallelepipeds 



n k = {x = (si, . . . , x n ) e R n I Xj £]aj k ,b jk [}, k £ N, 



such that 

OO 

(1) Ac U n fc ; 

k = 1 
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(2) m(I I fc) < e, where /r(IIfc) is the measure of the paral- 

k= 1 

OO 

lelepiped, i.e., g(n k ) = J] (b jk - a jk ). 

f=i 

8.2. Definition. We say that a property P(x) depending on 
the point x £ fl is valid almost everywhere, if the set of points x, at 
which P(x) does not hold, is of zero measure. We also say that P{x) 
is valid for almost all x £ fl. 

8.3. Definition. By a step function in fl we mean a function / : 
fl — > R that is a finite linear combination of characteristic functions 
of some parallelepipeds 11*,, k = 1, . . . , N, N £ N, i.e., 

N 

f(x) = y^cfc • ln fc (ff), c fc e R, x £ fl. (8.1) 

fc= l 
N 

In this case the sum c k • *), denoted by / /, is called the 

fc=i 

integral of the step function ( j8.1| ). 

8.4. Definition. A function / : Q, B x i — > /(x) £ C with 
complex values finite for almost all x £ fl is called measurable , if 
there exists sequences {g m } and {h m } step functions in such that 

lim [g m {x) + ihm{x)) = f(x) for almost all x £ fl. 

m — »oo 

8.5. P. Show that if f and g are measurable in ft and h : C 2 — > C 
is continuous , then the function fl 9 x i — > h(f(x),g(x)) £ C is also 
measurable. 

8.6. Definition. A set A c fl is called measurable, if 1^ is a 
measurable function. 

8.7. P. Show that any open set and any closed set are measurable. 
Show that the complement of a measurable set is measurable. Show 
that a countable union and a countable intersection of measurable 
sets are measurable. 

8.8. Definition. We say that / belongs to the class L + (more 
exactly, to L + (fl)), if in fl there exists an increasing sequence {h^^Lx 
of step functions such that f h m < C Vm for a constant C and, more- 
over, h m f /. The last condition means 3 ) that h\{x) < /^(z) < • • • < 
h m {x) < . . . and lim h rn (x) = f(x) for almost all x £ fl. 
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3) The notation h m J. / has the similar meaning. 

One can readily prove 

8.9. Proposition. If f £ L + (Q), then f is measurable in Cl. 



8.10. Definition. The (Lebesgue) integral of a function / £ 
L + ( Cl) is defined by the formula / / = lim / h m , where {h m }m=i 

m — kx) 

is an increasing sequence of functions determining / (see Defini- 
tion 



8 . 8 ). 



One can show that Definition 



8.10 



is correct, i.e. , / / depends 
only on / but not on the choice of the sequence h rn j /. 



8.11. Lemma (see, for instance, |56|h Let f n £ L + , f n f /, 
/ fn < C Vn. Then f £ L + and J f = lim f n . 

n — >-oo 

8.12. Definition. A function / : fl — > it. is said to be Lebesgue 
integrable (in f 1), if there exists two functions g and h in L + such that 
/ = g — h. In this case, the number f g — f h denoted by f n /( x) dx 
(or simply, J f ) is called the Lebesgue integral of the function /. The 
linear (real) space of functions integrable in 0 is denoted by L(Q) 
(or L) 4 ). If 1^4 £ L(fl), then the number g{A) = f 1^ is called the 
measure of the set A C LI. 



In this case two functions are identified, if their difference is equal 
to zero almost everywhere (see, in this connection, P |8.14| and P 8.18| ). 

8.13. Remark. One can readily check that / / depends only on 
/, i.e., Definition |8.12| is correct. 

8.14. P. Show that f / = 0, if f = 0 almost everywhere. 

N 

8.15. P. Let h^ix) = nik ■ lu k (x), where TO& = inf f{x) is a 

k = l 

N 

step function corresponding to the lower Darboux sum m fc/t(IIfc) 

k = 1 

of a function f which is Riemann integrable. Verify that Definitions 



8.8 8.10 and 8.12 immediately imply that f £ L + (thus, f £ L), 

N 

and, moreover, the Riemann integral equal to lim mkpijlk), 

fc=1 

coincides with the Lebesgue integral J f(x)dx = lim J h^(x)dx. 

N — >oo 
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8.16.P. Let f be measurable on [0,1] and bounded: m < f{x) < 
M. Consider the partition of the closed interval [ m , M ] by the points 
Vk ■ Vo < Ui < ■ ■ ■ < Un = M. Let a = ma x(y k - y k _ i). Consider 
the sum 

N 

s = ^2 y^{ x G [°> *] l 2 /fc- 1 ^ f( x ) ^ z/fc}- 

k= 1 

Prove that 3 lim S, and this limit is the Lebesgue integral f f(x) dx. 

a — >-0 



8.17. Theorem (Beppo Levi, 1906; see, for instance, [56] ) . Let 
f n £ L(ft) and f n (x) f f(x) Wx £ f 1. If there exists a constant C 
such that f f n < C Vn, then f £ L(Cl) and lim J f n = f f . 



8.18. P (Compare with P 8.14| ). Show that f 
where if f > 0 and f f = 0. 



0 almost every- 



8.19. P. Verify that if ip £ L, if £ L, then ma £ L, 
min V*) £ L. 



8.20. Theorem (Lebesgue, 1902; see, for instance, [56]). Let 
f n £ L(fl) and f n {x ) — » f{x) almost everywhere in Q. Suppose 
that there exists a function g £ L(f2), which is called the majorant, 
such that \f n (x)\ < g{x) Vn > 1, Va; £ f l. Then f £ L(Q) and 

//= lim I fn- 
n—> oo 



8.21. Lemma (Fatou, 1906; see, for instance, 56 i. Let f n £ L, 
f n > 0 and fn — * f almost everywhere. If f f n < C Vn, where 
C < oo, then f £ L and 0 < J f < C. 

8.22. Theorem (Fischer and F. Riesz, 1907). A space L with 
the norm ||</3|| = f \<p\, is a Banach space 5 \ 



5 ^ A norm in a linear space A is a function || • || : X 9 / i — > ||/|| € R, 
with the following properties: ||/|| > 0 for f A 0 € X, ||0|| = 0, ||A/|| = 
|A| • ||/|| for any number A, ||/ + g\\ < ||/|| + ||g||. The words “the space 
X is endowed with the norm” mean that a notion of the convergence is 
introduced in the space X. Namely, f n —■ ► / as n — > oo, if p(f n ,f) —> 0, 
where p(fn, f) = ||/n — /||. In this case one say that the space X is 
normed. The function introduced p has, as can be easily seen, the following 
properties: p{f,g ) = p(g,f), p(f,h) < p(f,g ) + p{g,h), and p(f,g) > 0, 
if f A g- If a function p with these properties is defined on the set 
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X x X, then this function is called the distance in X, and the pair (X,p) 
is called a metric space (in general, non-linear). It is clear that a normed 
space is a linear metric space. A metric space is called complete, if for 
any fundamental sequence {f n }n> 1 (this means that p(/n,/m) — > 0 as 
n, m — > oo ) there exists / £ X such that p(f n ,f) — > 0. A complete 
normed space is called a Banach space. 

PROOF. Let || ip n — tp\\ — > 0 as n, m — > oo. Then there exists 
an increasing sequence if indices {rife}fe>i such that \\ip n — tp m \\ < 

N - 1 

2~ k Vn > rife. Let us set fw(x) = |<^ nfc+1 (a;) - <Pn k (x)\. The 

fe= l 

sequence {/W}^ = 2 is increasing and / /jv < 1- Using the theorem of 

OO 

B. Levi, we obtain that the series ^ |^n fc+ i(^) — Vn k (%) | converges 

fc=i 

oo 

almost everywhere. Therefore, the series ^ |(/? rt) . +1 (a;) — </J nii .(:r)| 

fc= l 

also converges almost everywhere. In other words, for almost all x, 
the limit lim ip nk = <p{x) exists 6 ^. Let us show that ip £ L and 

k — »oo 

\\ip n — yj|| — > 0 as n — » oo. We have Ve > 0 3 N > 1 such that 
J\<Pn m (*r) — (pri k (x')\dx ^ 6 for n^n 3^ TV, rife ^ TV. Using the Fatou 
lemma, we pass to the limit as n m — > oo. We obtain ip — £ L , 

J |</?(:r) — yi nfc (a:)| dx < e; therefore, ip £ L, and ||^ — > 0 as 

k — * oo. Hence, ||<p — <p n || — > 0 as n — > oo, because ||y> — ip n || < 
\W~ Vn k W + \Wn k - Vnli# □ 

6) Thus, any fundamental sequence in L contains a subsequence that 
converges almost everywhere. We shall use this fact in Corollary |9.7| and 
in Lemma ll0.2l 

8.23.P. Construct an example of a bounded sequence (compare 
with the example at the beginning of Section^ which is fundamental 
with respect to the convergence defined by the Riemann integral but 
has no limit with respect to this convergence. 

Hint. For A £ [0, 1[, consider the sequence of characteristic func- 
tions of the sets C n that are introduced below, when constructed 
a Cantor set of measure 1 — A. Let {A ra } be a sequence of posi- 

OO 

five numbers such that J2 2" -1 A n = A £ [0,1]. The Cantor set 

n= 1 
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C C [0, 1] (see [21 [ [36 , 56] ) corresponding to the sequence {A„} is 

n 

constructed in the following way: C = fl C n , C n = [0, 1] \ ( (J /„), 

m = 1 

2 n ~ 1 

/„ = |J Im- Here, 1^ is the kth interval of the rath rank, i.e., an in- 

fc= 1 

terval of length A m whose centre coincides with the centre of the £:th 
(k = 1, . . . , 2 m_1 ) closed interval of the set C m _i. In other words, 
the Cantor set is constructed step by step in the following way. At the 
first step we “discard” from the closed interval Co = [0, 1] its “mid- 
dle part” Ji of length Ai, at the second step from the two remaining 
closed intervals of the set Cl = Co \ I we “discard” their “middle 
parts” each of which has the length A 2 . At mth step (m > 3) we 
“discard” from the remaining 2 m_1 closed intervals of the set C m _i 
their “middle parts”, each one the length A m . It can be easily seen 
(verify!) that the set C is measurable and its measure y,[C) is equal 
to/x=l — A>0. It can be rather easily shown that the set C is not 
countable and, moreover, there exists a one-to-one correspondence 
between C and R. At first sight, this seems astonishing, since the 
measure /i(C) = 0 for A = 1. Nevertheless, it is true! (Prove.) The 
Cantor set is often used as the base of constructing some “puzzling” 
examples. 

8.24. Theorem (Fubini, 1907; see, for instance, [58]). Let ft x 
be an open set in R fc and Cl y an open set in R m . Suppose that 
f : Q, 3 (x,y) 1 — > f(x,y) is an integrable functions in the direct 
product ft = ft x x fly. Then 

(1) for almost all y £ fl x (respectively, x £ ft y ) the function 
/(•, y) '■ fl x £> x 1 — > f{x,y) (respectively, f(-,y ) : ft y 9 
x 1 — » f{x,y)) is an element of the space L(fl x ) (respec- 
tively, L(fly)); 



(2) 


f{x,-)dx £ L(fly) (respectively, 


[ f(-,y) d y e L(fi x )),- 


J 

fi, 

( 3 ) 


f{x, y)dxdy = / 


/ f(x, y)dx 


a- « 

II 


1 

Ss 

sit 

1 


J 

a 


J 

Cly 


J 

Itc J 




Q,y 



dx. 
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8.25. Remark. The existence of two (iterated) integrals 






f(x, y)dx 



dy and 



f(x, y)dy 



dx 



implies, in general, neither their equality nor the integrability of the 
function / in x Tl y (see, for instance, [ 21 ]). However, the 

following lemma holds. 

8.26. Lemma. Let f be a function defined in fl = Ll x x fl y . Sup- 
pose that f is measurable and f > 0. Suppose also that there exists 



the iterated integral 



/ 


/ f{x,y)dx 


Q,y 





dy = A. Then f G L(Tl x xfl y ); 



therefore, property 3) of Theorem 8.24 holds. 



Proof. We set f m = min (/, H m ), where H m = ma x(/i 1; . . . , h m ), 
{hk} is a sequence of step functions such that hk —> f almost every- 



where (see Definition 8.4 1 . Note that f m = lim min(/tfc. H m ) almost 

k — ^oo 

everywhere and |/ m | < \H m \, since / > 0. Therefore, f m = lim gk m 

k — »oo 

almost everywhere, where gkm = max(min(/ifc, H m ), — \H m \). Fur- 
thermore, \gkm\ < \H m \ G L Vfc > 1. Hence, by the Lebesgue 
theorem, / G L. By virtue of the Fubini theorem, we have 



dy < A. 



Let us note that f n | /. Therefore, by the B. Levi theorem, / G 

L(n). ' □ 



[ fm= [ 


/ fm(x, y)dx 


J J 

Qy 


1 / 



8.27. Theorem (see, for instance, [ 56 ]). Let f G L(M), g G 

L(K), F{x) = f f(t)dt, G(x ) = f g(t)dt. Then 
o o 

b b 

[ F(x)g(x)dx+ ( f{x)G(x)dx = F{b)G(b) - F(a)G(a). 
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In this case the function F( x) has for almost all x £ R the deriv- 
ative F'(x ) = lim(.F(:r-|- <r) — F(x))/a, and F'{x ) = f(x) almost 

( 7 — >0 

everywhere. 

9. The spaces L p and L/ oc 

9.1. Definition (F. Riesz). Let 1 < p < oo. By the space 
L P (Q) (or simply L p ) of functions integrable in ptli power we call 
the complex space of measurable functions 1 ^ / defined in fl and such 
that \f\ p £ L(£l). If / € L 1 ( fl), then the integral of / is defined by 
the formula 



11 More exactly, of the classes of functions {/} : Q. — ► C, where 
g £ {/} 4==£- g = f almost everywhere. 

9.2. Lemma. Let p £ [l,oo[. Then the mapping 



i/p 



: L p 9 / , 



\f(x)\\ p dx 



(9.1) 



which will sometimes be noted by || • \\lp is a norm 2 \ 

^ That is the properties of a norm (see note 5 in Section |8| hold. 



PROOF. It is not obvious only the validity of the triangle in- 
equality, i.e., the inequality 

11/ + slip < Wf\\p + IMIp (9-2) 

which (in case of norm (9.1 1 ) is called the Minkowski inequality. It 
is trivial for p = 1. Let us prove it for p > 1, using the known 
Holder inequality 



11/ -sill < II/IIp ' llffllg, where 1/p + 1/q = 1, p > 1. (9.3) 

We have 

J i/ + g\ p < J i/+«?r 1 i/i + /(i/+<?r 1 i<?i) 




9. THE SPACES L p AND L‘ 



41 



However, 



C "i r r 

J 1/ + 9\ (p ~ iyq = J 



\f + 9\ l 



i-(i/p) 



□ 



Similarly to the proof of Theorem 8.22 one can prove 



9.3. Lemma. Let 1 < p < 00 . The space L p with norm (9.1 1 is 
a Banach space. 

9.4. Lemma. The complexification of the space of step functions 3 ^ 
is dense in L p , 1 < p < 00 . 

3) The complexification of a real linear space A' is the complex linear 
space of elements of the form / = g + ih, where g and h are elements of 
X. 



PROOF. It is sufficient to prove that V/ G L p , where / > 0, 
there exists a sequence {h k } of step functions such that 



11/ - h k \\ P -*• 0 as k — > 00 . 



(9.4) 



In case p = 1, we take the sequence {h^^Lx such that hk T / £ L + 
and f hk — > / /. Then we obtain ( 9.4 ( . If 1 < p < 00 , then we 
set E n = {x G ft | 1/n < /(x) < n}, where n > 1 and /„(x) = 
1-E„(a0 • /(x) (l_g„ is the characteristic function of E n ). We have 
/n T /; hence, (/ - f n ) p | 0. By the B. Levi theorem, ||/ - f n \\ p = 

(In \f(x) - f n (x)\ p dx) 1/2 -> 0 for n — > ex). Therefore, Ve > 0 3n > 1 
such that ||/ — /„|| p < e/2. Let us fix this n. Note that f 1 e k = 
f 1 e„ — f n p \f\ p < 00 . By virtue of the Holder inequality, J f n = 

f 1e„/ < (/ 1 | n ) 1/ ' ? /(/ p ) 1/p < 00 . Since f„ G L(Q) and f n (x) G 
[0,n] Vrefi, there exists a sequence {/ifc} of step functions defined 
in with values in [0, n] such that lim f \f n — hk\ = 0. Therefore, 

k — »oo 



||/« - h k \\p = 



I f n h k \ P 



1 !/P 



f I 1 /P 

J (| fn - h k \ P ~ 1 \fn - h k |) 



< n 



I/™ - M 



i/p 



0 for k — > 00 . 
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Choose K such that ||/„ — hk\\ p < e/2 for k > K. Then 

Il/-Mp< ll/-/n||p+ll/n-M <e Vfc > K. 



□ 



9.5. Theorem. Let f £ L 1 (f2) and / = 0 almost everywhere 
outside some K <s Cl. Let p > 0 be the distance between I\ and dCl. 
Then for any e €]0,p], the function 

Re(f) = fe-ClB xi — > f e (x) = j f(y)S e (x - y)dy , (9.5) 



where is defined in (3.2 1 , belongs to the space Cfi°(fl). Moreover 
Jim II/- / £ ||p = 0, l<p<oo. (9.6) 



4) Function ( |9.5[ ) is called the (Steklov) smoothing function of the 
function /. 



PROOF. Obviously, f e £ Cfi?(Cl). Let us prove (9.6). By virtue 
of Lemma |9.4| Viy > 0 there exists a function h = hi + i/ 12 , where 
h\ and /i 2 are step functions, such that \\f — h\\ p < 77 . We have: 
11/ - /clip < 11/ - h\\p + lift - -R<f(ft)||p + ||-R £ (/ - ft)llp- Let s show 
that ||i? e (< 7 )|| p < ||g||p. For p = 1 this is obvious: 



\g(y)\ ■ S e {x - y)dy 



n Lo 



n Ln 



dx 



S e (x — y)dx 



\g{y)\dy= / \g{y)\dy. 



9. THE SPACES L p AND L- 



43 



If p > 1, then by inequality (9.3): 
Il#e(s 0 ||£= / 1 9 e{x)\ p dx 



< 



< 



1 p 



(6 e (x - y)) (p 1)/p) (S e ( x - y) 1/p \g(y)\)dy 



dx 



(p-i )/p / \ i /p 

5 e {x - y)dy | f J 6 t (x - y)\g{y)\ p dy j 



V 

dx 



de{x - y)\g(y)\ p dy 



dx 



S e (x — y)dx 



\g{y)\ p dy= / \g(y)\ p dy. 



Thus, ||/ - f e \\ p < 2y + || h - R e {h ) || p . By virtue of (8.1), h = 



N 



^2 Ck • ln fc , where c& G C; hence, 



k = l 



\h — Re(h)\\p = 



I I Cfc ' “ ^(inj) dx 

a k=1 

< l Cfc ^ m f x / l ln * “ R e{ln k )dx <Ce, 



because (ln fc — i? e (ln fc )) = 0 outside the e-neighbourhood of the 
parallelepiped life. Taking e < y p /C , we obtain \\h — R t (h)\\ p < 
y. □ 

9.6. Corollary. is dense in L P (TL), 1 < p < oo. 

Proof. Let g £ L p (fl). Note that V 77 > 0 3K <s such that 
|| g — g ■ 1a || p < V- By Theorem 9.5 there exists e > 0 such that 
\\g ■ 1 K - Re(g ■ 1 a)||p < V- □ 

9.7. Corollary. Let f £ L 1 ( f2), / = 0 almost everywhere 
outside K <g Cl. Then there exists a sequence of functions f m £ 



44 1. INTRODUCTION TO PROBLEMS OF MATHEMATICAL PHYSICS 



Cq°(£ 1) such that / m — > / almost everywhere as m — » oo, if \f\ < M 
almost everywhere. 



Proof. By virtue of (9.5)-(|9l5|), ||/ — i? e (/)||i — > 0 as 



» 0 . 

Therefore, according to note 6 in Section [8] there exists a subse- 
quence {fm} of the sequence {R e (f)} e -, o such that / m — > / almost 
everywhere. The estimate \f m \ < M is obvious. □ 



9.8.P. 5 ) Prove that ||zi — R e (u)\\c —>0 as e — > 0, if u e Co(O). 



Here and below, ||/||c = sup \f(x)\ for / £ C(O). 



9.9. Definition. L°°(fl) is the space of essentially bounded 
functions in Q, i.e. , the space of measurable functions / : f2 — > C 
such that 



= inf sup \f(x)\ < oo, p(fl\u) = 0. (9.7) 



Condition (9. 7 1 means that the function / is bounded almost every- 
where, i.e., 3 M < oo such that \f{x)\ < M almost everywhere and 
ll/llno = inf M. 



One can readily prove 



9.10. Lemma. The space with the norm (9.7 ) is a Banach 

space, 



9.11. Remark. The symbol oo in the designation of the space 
and the norm (9.7 1 is justified by the fact that ||/||oo = lim 

P — KX 

eg R ra . This fact is proved, for instance, in m 



up) 



if 



9.12. Definition. Let A' be a normed space with a norm || • ||. 
Then X' denotes the space of continuous linear functional on X. 
The space X' is called dual to X. 



One can readily prove 

9.13. Proposition. The space X' equipped with the norm 

ll/ir=sup^^, /ex', 

I6X INI 

is a Banach space. Here, (f,x) is that value of f at x e X. 
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9.14. Theorem (F. Riesz, 1910). Let 1 < p < oo. Then ( L p )' = 
L q , where l/p+ 1 /q =1 (q = oo for p = 1). More exactly: 

1) V/ G L q (Ll) 3 F G {L p {Ll)) r , i.e., a linear continuous func- 
tional F on L P {Q) such that 

( F ,<fi) = J f(x)tp(x)dx \/tpGL p (n); (9.8) 

n 

2) VF G ( L p (fl ))' there exists a unique ^ element (function) 
f G L q {fl) such that ( |9.8| ) holds; 

6 - ) See note 4 in Section 0 

3) the correspondence I : ( L p )' 9 F i — > f G L q is an isomet- 
ric isomorphism of Banach spaces, i.e., the mapping I is 
linear bijective and ||/F|| g = ||T||p. 

Proof. Assertion 1) as well as the estimate ||F||p < ||/|| 9 are 
obvious for p = 1. For p > 1 one should use the Holder inequality. 
Assertion 2) as well as the estimate ||F||p > ||/|| g are proven, for 
instance, in m- Assertion 3) follows from 1) and 2). □ 

9.15. Definition. Let p G [l,oo]. Then Lf oc (SY) (or simply 

L p oc ) denotes the space of functions locally integrable in pth power 
f : H — » C, i.e., of the functions such that / • Ik G L p (Ll) \/K g H. 
We introduce in L v loc (Ll) the convergence: fj — » / in if and 

only if || lie ■ {fj ~ f) || P 0 as j -> oo \/K <g Q. 

Let us note the obvious fact: if 1 < r < s < oo, then 
PC C Lf% c C L s loc C L\ oc C L\ oc . 

10. Functions of L] oc as linear functional on Cq° 

The idea of representability (i.e., of determination) of a function 
by its “averagings”, outlined in Sections m can now be written in 
a rather general form as the following 

10.1. Theorem. Any function f G L\ 0 JP) can be uniquely 1 ’ 
reconstructed by the linear functional 

(f,<P) = / f{x)ip{x)dx G C (10-1) 
n 



(/,•>: C 0 °°(H) 9 ip 
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(i.e., by the set of the numbers j (/, ip) Moreover, the 

correspondence f < — > (/, ■), f £ L} oc , is an isomorphism. 

11 As an element of the space L] oc (£l) (see note 4 in Section [8|. 

PROOF. Suppose that two functions /i and /2 correspond to 
one functional. Then f (/i — f 2 ) l f > = 0 S C [ (°. This, by virtue of 
Lemma [10. 2| below, implies that /i = /2 almost everywhere. □ 

10.2. Lemma. Let f £ Lf oc (Q). If f n f(x)<p(x)dx = 0 Vip £ 
Co c (ft), then f = 0 almost everywhere. 

PROOF. Let w <s Cl. Note that |/| ■ l w = / • g , where g{x) = 
l w • exp[— * arg f(x)]. (If / is a real function, then g(x) = — sgn f(x) ■ 
l^(j;).) According to Corollary |9.7[ there exists a sequence of func- 
tions ip n £ such that almost everywhere in fl f ■ ip n — > / • g 

as n — > oo, and \ip n \ < 1. Since f |/| = fnf'9 an d by the Lebesgue 
theorem f n f ■ g = lim f n f ■ ip n , we have |/ = 0, because 

Jo / ' Pn = 0. Thus, / = 0 almost everywhere in to. Hence, by 
virtue of arbitrariness of ui <s H, / = 0 almost everywhere in Q. □ 

11. Simplest hyperbolic equations. Generalized Sobolev 

solutions 

In this section we illustrate one of the main achievements of the 
theory of distributions on the example of the simplest partial differ- 
ential equation ut + u x = 0, 1 ' which is sometimes called the transfer 
equation. It concerns a new meaning of solutions of differential equa- 
tions, more exactly, a new (extended) meaning of differential equa- 
tions. This meaning allows us to consider correctly some important 
problems of mathematical physics which have no solutions in the 
usual sense. This new approach to the equations of mathematical 
physics and their solutions, designed by S.L. Sobolev in 1935 (see, 
for instance, ( 621 ) under the title “generalized solutions”, allows, in 
particular, to prove the existence and uniqueness theorem for the 
generalized solution of the Cauchy problem: 

Lu = u t +u x = 0, (x,t) € R+ = {{x,t) £R~ \ t > 0}, (11.1) 

u\ t=0 =f(x), ifR (11.2) 
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for equation (|11.1 1 for any function / £ PC'(R) (and even / g Lj nr ; 
see Theorem |11.10 below). The theorem is also valid (see P 11.11) 
on the continuous dependence of the solution of this problem on 
/ e L} oc (R). 

11 Here, Ut and u x denote the partial derivatives of the function u(x, t) 
with respect to t and x. 



Let us clarify the essence of the problem. Equation (11.1) is 
equivalent to the system 

Ut + u x ■ dx/dt = 0, dx/dt = 1. 

Therefore, along the line x = t + a, where a is a real parameter, we 
have du(t + a,t)/dt = 0. In other words, u(t + a,t) = u(a, 0) Vi. 
Thus, the function /( x) = lim u(x,t ) must necessary be contin- 
uous; in this case u(x,t) = f(x — t ). 2 ^ If / is differentiable, then 
u(x,t) = f (x — y) is a solution of problem (11.1) (11.2). How- 
ever, this problem has no solution (differentiable or even contin- 
uous), if / is discontinuous, for instance, if f{x) = 9(x), where 
i — > 9(x) £ R is the Heaviside function, i.e. , 

6(x) = 1 for x > 0 and 9{x) = 0 for x < 0. (H-3) 



2) This formula implies that the graph of the function x i — > u(x,t) 
for any fixed t can be obtained by the transition (shift) of the graph of the 
function / to the right along the axis x on the distance t. It is the reason 
why equation (11.1 1 is sometimes called the transfer equation. 



However, consideration of problem (11.1 1 (11.2) with the initial 
function (11.3 1 is justified at least by the fact that this problem arises 
(as minimum, on the formal level) when studying the propagation of 
the plane sonic waves in a certain medium. The appropriate process 
is described by the so-called acoustic system of differential equations 

ut + -p x = 0, p t + p ■ c 2 u x = 0, p > 0, c > 0. (11.4) 



Here, p is the density, c is the characteristic of the compressible 
medium, and u = u(x, t) and p = p{x, t) are the velocity and the 
pressure at the instant t at the point x. Setting 



a = u + p/{p-c), (3 = u-p/{p-c), 
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we obtain the equivalent system at + ca x = 0 , fit ~ c/3 x = 0 of two 



transfer equations. Thus, problem ( 11.1 1 ( 11.2 I with the initial func- 



tion (11.31 arises when considering the propagation of sonic waves, 



say, for the initial velocity u(x,t) = 9(x) and zero initial pressure. 
ll.l.P. Show that any solution of the class C 1 of system (11.4) 



can be represented in the form 

u(x, t ) = [tp{ x — ct) + iffx + ct)\/ 2, 

p(x,t)=[(p(x — ct) — if(x + ct)]/2, where ipGC 1 , ?/> £ C 1 . 

(11.5) 

11. 2. P. Show that the following theorem is valid. 

11.3. Theorem. V/ £ C 1 (K) MF £ C(]R+) the Cauchy problem 
u t + u x = F(x,t) in R+, zt| t=a = f(x), x £ R 

has a unique solution u £ C 1 (fi 2 + ). 

As has been said, for f(x) = 9(x) problem (11.1 1— ( flL2| has 
no regular solution (i.e., a solution in the usual sense of this word); 
nevertheless the arguments which lead to the formula u(x, t) = f(x — 
t) as well as this formula suggest to call by the solution of problem 
(11.1 1— ( 11.2 1 the function f(x — t) for whatever function / £ PC(K) 
(and even / £ L\ oc ( R)), especially, as the following lemma holds. 



11.4. Lemma. Let f £ Lj oc (M.) and {f n } be a sequence of func- 
tions f n £ C^R) such tha& 



fn -> / in Lj a 



as n 



oo. 



Then the function u : 9 (x, t) i — * u{x, t) = f{x — t) belongs to 

Lj oc (M.+), and u = lim u n in I/) oc (M^_), where u(x,t ) = f(x — t). 4 > 



3 ' According to Lemma 11.5 below, such a sequence exists. 



4) Note that «„(i, t) is a solution of problem ( 11.1 ) — ( 11.2 1, if Mn| t _ 0 = 
fn(x). 



PROOF. Since f = f 1 + if 2 and f k = f k — f k , where f± = 
max(±/ fe , 0), it is sufficient to consider the case / > 0. Let us make 
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the change of the variables (x,t) i — » where y = x — t. Note 

that u(x,t ) = f(y) and 



b / d 



b / d—a 



j | J u{x,t)dx J dt < f | J f(y)dy J dt < oo 



a \c—d 



for any a, 6, c, d such that 0 < a < b, c < d. By virtue of Lemma [8. 26 [ 
this implies u £ L) oc (R+). Furthermore, for the same a, b, c, and d 



b / d 



J | J | U n (x,t) - u(x,t)\dx J dt < (b - a) j \ f n (y) - f(y)\dy -»■ 0 



a \c 

as n — > oo. 



c—b 



□ 



11.5. Lemma. V/ £ Lj oc (M) there exists a sequence of functions 
f n £ C°°(R) converging to f in Lj oc (R). 

Proof. Let {w}/j^=i be the partition of unity in R (see Sec- 
tion [ 3 }, Vm € Cg°(R) and ip^ ' = W We have Vv ' / e 

By Theorem 9.5 there exists a sequence of functions f% £ 

C 00 ^) such that for any fixed y: lim \\ifnf — f £ ||i = 0. Setting 

n — >-oo 
00 

/ n (x) = ^2 (x), x G R, we have f n G C°°(R). Note that 

m=i 

OO M 

Vc > 0 3M G N such that J2 wO®) = S < Pn( x ) for M < c. Hence, 

fl— 1 /i=l 



l/fa) - fn{x)\dx = 



M 



f{x) 






dx 



+ 



M 



- fn) 

n=l 



M J; 

dx -J2j - fn\ dx - 






Therefore, 



M 

lim f I/- /„|cte < V lim [ \tp^f - f£\dx = 0, 

i — >00 / z ' n — >-oo / 

u— 1 

— r. ^ — r 



because lim ||Vv/ ~ /£||i = 0. 



□ 
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The definition of the solution of problem ( ll.l| )-( 11.2 ), where 
/ € Lj oc , with the help of the formula u(x , t ) = f(x — t) is rather 
tempting, however, let us note that it has a serious defect: with 
the help of a concrete formula, one can define the solution of only 



a small class of problems. Lemma 11.4 suggests a definition free of 
this defect. 

11.6. Definition. Let / G Lj nr { M.). W e say t hat u G 
is a generalized solution of problem (11.1) (11.2), if there exists a 
sequence of solutions u n G C 1 (R^_) of equation (11.1) such that, as 
n — > oo, 

u n ->u in L} oc (R 2 + ) and u n \ t ^ Q -> / in L} oc (R). 

Approximative approach to the definition of a generalized so- 
lution can be applied to a large class of problems. So, it has been 
above constructed (but was not named), for instance, the generalized 
solution of the equation A E(x) = 6(x) (see ( |7.9| )) as well as the gen- 
eralized solution of the problem A P = 0 in Mil, P{ x,0) = S(x) (see 



Remark 5.4 1. However, the approximative definition, in spite of tech- 



nical convenience, also has an essential shortage: it does not show 
the real mathematical object, the “generalized” differential equation, 
whose immediate solution is the defined “generalized solution”. 

It is reasonable to search for the appropriate definition of the 
generalized solutions of differential equations (and appropriate “gen- 
eralized” differential equations), by analyzing the deduction of the 
equations of mathematical physics (in the framework of one or an- 
other conception of the continuous medium). The analysis fulfilled 
in Sections 013 Lemma [l0.2[ and the Ostrogradsky-Gauss formula 
(7.2) suggest the suitable definition (as will be seen from Proposi- 



tion 11. 



11.7. Definition. Let / G L^^M). A function u G L ; 1 oc (R^_) is 
called a generalized solution of problem (11.1 |-( [il.2| ), if it satisfies 
the following equation (so-called integral identity (in (/?)) 

(tpt + ipx)u(x,t)dxdt + / ip(x,0)f(x)dx = 0 \hp G C'g(R^). 






( 11 . 6 ) 

11.8. Proposition. Ifu G C 1 (1R^_), then ( 11.6 ) is equivalent to 



(11-1) (1L2). 
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PROOF. Let ip £ Cq(R+) and let Obea bounded domain in R^ 
with the boundary T = dfl. Formula (7.2 1 implies that 

J (u t + u x )ip dx dt + J (<p t + ip x )udx dt 

= / (tp • w)[cos(u, t) + cos(u, x)] dr. (11.7) 
on 

If slippy C D and (see Fig.) (supp ip D dfl) cRi = {(x,i) G R 2 | 
t = 0}, then formula (11.71 can be rewritten in the form 

J (ut + u x )ipdxdt + J (<pt + <p x )udxdt = — J (ipu)\ t _ Q dx. ( 11 . 8 ) 



R2 



Furthermore, by virtue of Lemma |10.2[ 



( 11.11 



(ut + u x )ipdxdt = 0 \/ip G C'g(R^) 



and 



( 11.21 



f(x)ip(x,0)dx = / u\ t _ Q ■ ip(x,0)dx V(p £ Cq 



This and (11.81 imply that ( 11.6 !<*=>( 11.1 1 (11.21. 



Proposition 11.8 shows that 
Definition [TT7T] is consistent with 
the definition of an ordinary (dif- 
ferentiable or, as one says, regu- 
lar) solution of problem (11.1 1 



(11.21. The following Theorem 
11.101 justifies the new features 






□ 




appearing in Definition EH and shows that the integral equality 
( 11.61 is the same “generalized” differential equation which has been 
spoken about. 

11.9. Remark. The proof of Proposition |11.8| comes from the 
deduction of the Euler-Lagrange equation and the transversality 
conditions in calculus of variations proposed by Lagrange (see, for 
instance, [56)). 

11.10. Theorem. V/ e 
(unique) generalized solution u £ L\ oc y^ +) 



problem (11.1 1 (11.2) has a 



,2 -)- 
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PROOF. First, we prove the existence. Since the function 
u : R+ 9 (x, t) i — ♦ u n (x, t) = f(x — t) 



is a regular solution of equation (11.1[) and satisfies the initial con- 



dition u n | ( _ 0 = /„(x), by virtue of proposition 11.8 
/ ( ipt + (p x ) ■ u n dxdt+ / f n (x)(p(x, 0) dx = 0 



we have 



G Cq(k+). 



(11.9) 

On the other hand, by Lemma 11. 4[ the sequence tends in 

L; 1 oc (R+) to the function u € L) oc (R 2 ) such that u(x,t) = f(x — t). 
It remains to verify that the function u satisfies (11.61. For this 
purpose we note: \hp £ Cq (R^) > 0 and b v > 0 such that 

Slippy C {(x, i) £ M 2 



fs u^,,0 / . 



Therefore, 

(u„(x, t) - u(x, t))(ip t + (p x )dxdt 



< [max | tp t + ip x |] 

(x,t) 



< M v ■ b v 




dt 



j \fn(x - t) - f(x ~ t)\dx 

~ CLip 

\f n {x) — f(x)\dx — * 0 for n — * oo. 




Vt + ‘fix = 9 (x,t), (x,t)GM. 2 + 



sufficient to show that the equation 

( 11 . 10 ) 

has a solution ip £ Cq(K^_) for any g £ Cg°(M^_). However, this 
follows from P 11.2 Indeed, let T > 0 be such that g(x,t ) = 0, if 
t> T. We set 
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t 

ip(x,t) = J g{x — t + t,t)cIt. 



Obviously, (see Fig.) ip G and p is a solution of ( 11.10 1. □ 



ll.ll.P. Prove that the generalized solution of problem (11.1 ) 



(11.21 depends continuously in L ; 1 oc (K^_) on the initial function f G 



L 



loc\ 



11.12. P. Analyzing the proof of Theorem \11.10\ prove that Def- 
inition \11. 7| is equivalent to Definition \11.0] 

11.13. P. Verify directly that the function u(x,t) = d(x — t) is a 
solution of problem ( 11.1 1 — ( 11.2 1 in the sense of Definition \ 1 1 . T\ if 
fix) = 6(x). 

In the exercises below, we assume Q = {(x,t) G K 2 | x > 0, 
t > 0}. 



11.14.P. Consider the problem 

u t +u x = 0 in Q , 




(11.11) 



U\ t=0 = fix), x > 0, (11.12) 



u=/ 



I x—0 



= h^t), t> 0. (11.13) 

This problem is called mixed , because it simultaneously includes the 
initial condition (11.12) and the boundary condition ( 11.13[ ). Show 
that problem ( 11.11 1 -( 11.13 ) has a (unique) solution u G C 1 (Q) if 
and only if f G C i (K + ), h G C^M-).), ans /( 0) = h(0), /'( 0) = 

-^(o). 



11.15.P. Show that the following problem 




u=f 



X 



u t u x — 0 in Q, 



u\ t=0 = fix), x>0 



(11.14) 

(11.15) 
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has a (unique) solution u £ C 1 (Q) if and only if f £ C 1 (]R+). Com- 
pare with problem ( 11.11 ) —( 11.13| ) . Compare the characteristics, i.e., 
the families of the straight lines dx/dt = 1 and dx/dt = —1 (shown in 
figures), along which the solutions of equations ( 11.11| ) and (11.14) 
are constant. 



11.16. P. Consider a mixed problem for the system of acoustic 
equations 

u t + (l/p)p x = 0, Pt + pc 2 u x = 0, (x,t) £ Q, (11.16) 

u\ t=0 = f(x), p\ t=0 = g(x), x>0 , (11.17) 

p| x=0 = M*)’ 1 > (11.18) 

where f, g, and h are functions from C 1 (1R+). 

(1) Draw the level lines of the functions u ± (1 / pc)p. 

(2) Show that problem (11.16)— (11.18) has a (unique) solution 
u £ C 1 (Q ), p £ C 1 (Q) if and only if 

h(0) = g(0) and /'(0) + (l/pc 2 )h'(0) = 0. (11.19) 

Show also that this solution ( u,p ) can be represented by formulae 
(11.51, where 



v(y) = f(y ) 



and 



(1 /pc)g(y), Ip(y) = f(y) - (1 /pc)g{y), if y > 0 

( 11 . 20 ) 



¥>(y) = (2 /pc)h(-y/c) + f(-y) - (1 /pc)g(-y), if y< 0. (11.21) 



11.17. Remark. Often instead of system (11.41 of the acoustic 
equations, the following second order equation is considered 

d 2 p/dt 2 - c 2 • d 2 p/dx 2 = 0. 



This equation obviously follows from system (11.4), ii p £ C 2 , u £ 
C 2 . This equation is called the string equation, since the graph of 
the function p can be interpreted as a form of small oscillations of a 
string. The string equation is a special case of the wave equation 

Pu — c 2 Ap = 0, p = p(x,t), x £ R", t > 0. (11.22) 



Here, A is the Laplace operator. For n = 2 the wave equation 
describes the oscillations of a membrane, for n = 3 it describes the 
oscillations of 3-dimensional medium. 
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11.18. Remark. The string equation is remarkable in many as- 
pects. It was the first equation in partial derivatives that appeared 
in mathematical investigations (B. Teylor, 1713). It was a source of 
fruitful discussion (see, for instance, [42, 1501 ) in which the notion of 
a function was developed (d’Alembert, Euler, D. Bernoulli, Fourier, 
Riemann,. . . ). 



11.19. Remark. Many distinctive properties of some differen- 
tial operators A(y, d Vl , . . . d Vm ), y £ R are defined by the proper- 
ties of the corresponding characteristic polynomials A(y, 771 , ... , r) m ) 
of the variable 77 = (771 , . . . ,77 m ). Thus, the hyperbolic polynomial 
t 2 — £ 2 (t 2 -|£| 2 ) is associated with the string equation u tt — u xx = 0 
(or, more generally, with the wave equation u t t — A u = 0); the ellip- 
tic polynomial |£| 2 = Yl£k associated with the Laplace equation 
k 

A u = 0; the parabolic polynomial t— |£| 2 is associated with the heat 
equation ut — A u = 0. According to the type of the characteristic 
polynomial, partial differential equations can be classified into hy- 
perbolic equations, elliptic equations, parabolic equation, .... (See 
exact definitions, for instance, in [48j .) 



11.20. Example. Consider problem ( 11.16 1 (11.18) with / = 
g = 0, h = 1. This means that at the initial instant t = 0 the 
velocity u and the pressure p are equal to zero, and on the boundary 
x = 0 the pressure p = 1 is maintained. Formulae (11.5 1, (11.20) 
(11.21 1 give the following result: 



u = 0, p = 0, 
u=l/pc, p= 1, 



if t < x/c 1 
if t > x/c J 



(11.23) 



The functions u and p are 
discontinuous that is not sur- 
prising, because condition ( 11.9 1 
does not hold. However, on the 
other hand, formulae ( 11.23 ) are 
in a good accordance with the 
physical processes. This yields 



p = 1 




u=0,p=0 



11.21.P. Find the appropriate definitions of the generalized so- 
lutions for the following problems: 
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(1) u t + u x = F(x, t) in Q = {(x, t) £ 



> 0 ,t> 0}; 



lt=0 



= f(x), x > 0; u\ x=0 = h(t ), t > 0; 



(2) u t + (1 / p)p x = F(x, t), p t + p- c 2 u x = G(x, t ) in Q; 

= f(x), p\ f ~ = g{x), x > 0; u\ = h(t), t > 0. 



it=o 



(3) ptt — c 2 p xx = F(x,t) in the half-strip f l = {( x,t ) £ Mi | 
0 < x < 1}; 

p| t=0 = f(x), Pt | t=0 = g(x), 0 < x < 1; p | x=Q = h 0 (x), 
p\ x=1 = fti(x), t > 0. 

Revise the requirements to the functions f, g , h, F, G, under which 
the solutions of these problems belong, say, to the space C 1 , PC 1 
or Lj oc . Prove the theorems of existence, uniqueness and continuous 
dependence (compare with P 11.11). 



We conclude this section by consideration of the non-linear equa- 



tion 



u t + (u 2 / 2) x = eu xx + b x (x, t), u = u(x, t), 



(11.24) 



where e > 0, and b is a given function. This equation is called the 
Burgers equation and is considered in hydrodynamics as a model 
equation for e > 0 of the Navier-Stokes system, and for e = 0 of the 
Euler system (see 1521 i. 

First, consider 




equation (11.241 for 



e = 0 and b = 0. In 
this case, the regu- 
lar (of the class C 1 ) 
solution of this equa- 
tion satisfies the sys- 
tem dx/dt = u, 
du/dt = 0. Thus, 
the solution u(x, t) 

is constant along the characteristic , i.e., along the curve defined 
by the equation dx/dt = u(x,t ); hence, this curve is, in fact, the 
straight line x = a + f(a)t that depends only on the parameter 
atl and a function /. The function / is determined by the rela- 
tions /(a) = dx/dt , dx/dt = u(x,t), i.e., f(x) = u( x,0). If / is a 
decreasing function, for instance, f(x) = ^th(x), then the charac- 
teristics intersect at some t > 0, and at the point of the intersection 
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we have 

u(x — 0, t) > u(x + 0, t). (11.25) 

The continuous solution cease to exist. Thus, the Cauchy problem 
u t + {u 2 /2) x = 0 in R+, u| n = f(x), (11.26) 



has, in general, no continuous solution even for analytic initial data. 
This effect is well known in hydrodynamics. It is connected with 
arising of the so-called shock waves which are characterized by a 
jump-like change of the density, velocity, etc. Thus, physics suggest 



that the solution of problem ( 11.26 1 should be sought as a generalized 
solution of the class PC 1 . 

^ See Addendum. 



Suppose that u is a generalized solution of problem (11.261, and 
u has a jump along the curve 

7 = {(x,t) G R 2 | x = A (t), A G C 1 [cv, /3]}; 



more exactly, suppose that (x,t) G 7 satisfies condition (11.251. 
11.22. P. Prove (compare with P\12.(fy that the Hugonidt condi- 



tion 



(11.27) 



d\(t) /dt = [u(A(t) + 0, t) + u(X(t) — 0, t)\/ 2. 
holds along this line 7 called the break line. 

One can show (see, for instance, 52j ) that relations (11.251, 



(11.271 replace the differential equation u t + (it 2 / 2) x = 0 on the 



break line. 



One of the approaches to the study of problem (111. 261) is based 



on consideration of the Cauchy problem for equation ( |11.24[ ) for e > 0 
(and 6 = 0) with the passage to the limit as e — > 0 (see, for instance, 



]). The point is that (for 6 = 0) equation (11.241 can be reduced, 
however surprising it is, to the well studied heat equation. Actually, 
the following theorem holds 6 ^. 

This theorem was proven in 1948 by V.A. Florin and in 1950 was 
rediscovered by E. Hopf. 



11.23. Theorem. The solution of equation (11.241 can be rep- 



resented in the form u = — 2e(ln G) x , where G is the solution of the 
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linear parabolic equation 



G* — eGn — 



b(x,t ) 
2e 



G. 



(11.28) 



PROOF. Let u be the solution. Setting 
P(x,t) = u(x,t), Q(x,t) = —u 2 (x,t )/ 2 + eu x (x,t) + b(x,t ), 

we have P t = u t , Q x = -u-u x + eu wx + b x (x,t). Therefore, P t = Q x . 
Thus, the function is defined 

(x,t) 

F(x,t) = / Pdx + Qdt. 

(0,0) 

We have F x = P, F t = Q. Hence, F t + (F x ) 2 /2 — eF xx = b. Introduc- 
ing the function G = exp[— F/2e], we obtain that G is the solution 
of equation (11.281 and u = — 2e(ln G) x , because u = F x . □ 



CHAPTER 2 



The spaces T>\ V # and V . Elements of 
the distribution theory (generalized 
function in the sense of L. Schwartz) 



12. The space V 7 of the Sobolev derivatives 

The definition of the generalized solution u £ Lj oc to one or an- 
other problem of mathematical physics given by Sobolev [61] and, 



in particular, Definition 11.7 



is based on Theorem 10.1 and the 



Ostrogradsky-Gauss formula (7.2 1 . Let us recall that Theorem 10.1 



asserts the equivalence of the two following representations of an 
element u £ L] oc : 



1)S19 



l{x) and 2) C£°(fl) 3 ip 



u(x)ip(x)dx, 



and formula (7.2 1 implies that, for the differential operator d a = 
/ dx . . . dx^™ and any function u £ the following iden- 

tity is valid: 

f (d a u(x))ip(x)dx = (— l)l“l f u(x)d a ip(x)dx Vip£C^°(fl). 



Thus, the functional 



d a u : CS° (Si) 3 ^ (d a u,<p) = (-1) H J u(x)d a <p(x)dx 

n 

Vp £ Cg°( Q) (12.1) 



determines the function d a u(x ), if u £ Cl a qfl). Since functional 
) is also defined for u £ Lj oc ( fl), one can, tracing S.L. Sobolev’s 
approach, give the following definition. 
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2. THE SPACES T)\ T>* AND T>' 



12.1. Definition. Let a = (a±, . . . ,a n ) be a multiindex. The 
weak derivative of order a of the function u £ Lj oc ( ft) is defined as 
the functional d a u given by formula (12.1). 



Using Theorem 10.1 formula (7.2 1, and the identity 



a(x)(d a u(x))ip(x) dx = (— 1)^ / u(x)d a (a(x)ip(x)) dx, 



«ec |a| (n), ip g c™(n), 



which is valid for any function a £ we introduce the opera- 

tion of multiplication of the functional d a u , where u £ Lj oc (Cl), by 
a function a £ C°°(C2) with the help of the formula 



ad a u : C£°(ft) 9 ip i — ■» (— 1)^ J u{x)d a (a(x)ip(x))dx £ C. (12.2) 

n 

12.2. Definition. The space of Sobolev derivatives is the space 

of functionals of the form d a u a , where a is a multiindex and 

|q:| <oo 

u ^ iioc(n), equipped with the operation of multiplication by for- 
mula This space is denoted by X> b (fl). 

12.3. Example. Let the function x+ £ L] oc { R) be defined in 
the following way: x+ = x for x > 0, x+ = 0 for x < 0. Let us find 
its derivatives. We have 



(x' + ,ip) = —(x + ,ip r ) = — J x + tp' (x)dx = — J xip'(x)dx 

R R + 

oo 

= ir + / = / ew* ><*)* =<«.*>>. 

0 R 

i.e., x’ + = 0 is the Heaviside function. Now find x'f , i.e. , 0'. We 
have 

OO 

(0', ip) = -(0, ip') = - J ip'(x)dx = -<p{x) |“ = ip(0) = (S, ip), 

0 

(12.3) 
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i.e., 0' = <5(x) is the Dirac (5-function. In the same way one can find 
any derivative of the 5-function of order k. We have 



(5 {k \p) = -(^V) = ■■ • = (-l) fc <<5,^) = (-l)V fc) (0). 

(12.4) 



12.4. P. Let O e £ C°°(K), 0 < O e (x) < 1, and O e (x) = 1 for 
x > e and O c (x) = 0 for x < — e. Let us set 5 e ( x) = 0((x). Show 
that lim e ^o $ k) ,<p) = (-l)V (fe) (0) Vfc >0,V^6 Cg°(R). 



12.5. Remark. Formulae ( |12.4[ ) allow us to extend the func- 
tional 5^ from the functional space C“(K) to the space of func- 
tions k - times continuously differentiable at the point x = 0 (see 
Definition |2)2]). On the other hand, formula (12.3) is not defined on 
the space C(R), because the functional 0 is not defined on 



Define the function O-t : R ra 9 x 
formula 

0±(z) = Iq ± (x), 

Q± = {x = (xi, . . . ,x n ) £ R" | ±Xk > 0 Vfc}. 



0 ±(x) by the following 



(12.5) 



If n = 1, then 0 + = 0 is the Heaviside function, and 0_ = 1 — 0 + 

( in L}oc). 



12.6. P. (cf. P 11.22). 
C 1 (R 2 ). Let, for (x ,t) e!lc 



Let F £ C' 1 (R), A e C^R), u± £ 
R 2 , u[x,t) = u + (x,t)@ + (x — A (t)) + 



U- (x,t)Q- (x — A(t)). Findut and ( F(u)) x , noting that F{u(x,t)) = 
F(u+(x, t))0+(x — A(t)) + F(u-{x , t))0_(x — A(t)). Show that ut + 
(F(u)) x = 0 almost everywhere in Ll if and only if, first, ut + 
( f(u)) x = 0 in Cl \ 7 , where 7 = {(x,f) £ R 2 | x = A(f)}, and 
second, the Hugoniot condition = F holds 

on the line 7 . 



12. 7. P. Check (see (12.51,) that 



d n 

dxi ...dx^ 



- 0 + = 5(x), x £ 



12.8. P. Show that the function E(x,t) = 0(f — |cc|) /2 is the 
fundamental solution of the string operator, i.e., 



(d 2 /dt 2 — d 2 /dx 2 )E(x, t) = 5(x, t). 



Here, S(x,t) is that 5-function in R x R, i.e., ( 6(x,t),ip } = 73 ( 0 , 0 ) 
V 73 g C'q°(R x R). 
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12. 9. P. Noting that, for ip £ C'^°(K), 



lim J In \x\ ■ ip'(x)dx = lim In e(y>(— e) — <p(e)) — J 

|a;| >e |x|>e 



ifi(x) 



dx 



prove that ^ln|x| = v.p.\, i.e., (^ln|x|,^) = v.p. 

\/ip £ C'g°(]R), where v.p. x~ 1 (p(x) dx is so-called the principal 
value = valeur principal (French) of the integral f)™ x~ 1 ip(x)dx 
defined by the formula 



OO 

v.p. J x~ 1 (p(x)dx = lim J x~ 1 (p(x)dx. 



M>e 



(12.6) 



12.10. P. Taking into account that ln(x ± ie) = ln|a; ± ie\ + 
iaig(x ± ie) — > In |x| ± *7r0(— x) as e — > +0, prove the simplest ver- 
sion the Sokhotsky formulae (very widespread in the mathematical 
physics (see, for instance, Eg;; 

= v.p.— ± inS(x), (12.7) 

x =F iO x 

i.e., prove that = V -P- f - a 0 m^dx ± imp^) \hp £ 

cg°(R)- 



12.11. Remark. Formulae (12.71 imply that 



S(x) = f(x — iO) — f(x + 10), where f(x + iy ) = (x + iy) 1 , 

2m 

( 12 . 8 ) 

i.e., the 5-function being an element of P^R) admits the represen- 
tation in the form of the difference of boundary values on the real 
axis of two functions analytic in C+ and in C_, respectively, where 
C± = {z = x + iy £ C j ±y > 0}. This simple observation has deep 
generalizations in the theory of hyperfunctions (see, for instance, 

mm)- 



12.12. Remark. Any continuous function F £ C(R) has, as 
an element of the space Lj oc , the Sobolev derivative F' £ 2T(M). 
If this derivative is a locally integrable function, in other words, if 
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F 0) = Sa f(y) d v + F ( a )> where / e L l 

implies that 



then Theorem 



1.27 



F' (x) = lim a 1 (F(x + a) — F(x)) for almost all ieK. (12.9) 

cr — >0 



In this case, formula ( 12.9 ) totally determines the Sobolev derivative 
F' . Emphasize that the last assertion does not hold (even under the 
assumption that formula (12.9) is valid), if F' ^ L i 1 oc (K). Thus, for 
instance, the Cantor ladder (see [36] or [56]) corresponding to the 



Cantor set of zero measure (see hint to P8.23), i.e. , a continuous 
monotone function K £ C[0, 1] with the value (2k — 1) • 2~ n i n fcth 
(k = 1, . . . , 2™ _1 ) interval I n =]a k , b k [ of rank n (see hint to P 8.23 ) 
has, for almost all x £ [0, 1], zero derivative but its Sobolev derivative 
K' is non-zero. Namely, 



oo 2 71 



K' = ^^(2k-l)-2- 

n—1 k—1 



\6(x-b k )-6(x-a k )). (12.10) 



12.13.P. Prove formula (12.101. 

13. The space T)^ 1 of generalized functions 

The elements of the space V'’ were defined as finite linear com- 
binations of the functionals d a u a (12.1), i.e., of the derivatives of 
the functions u a £ Lj oc . If we neglect the concrete form of the 
functionals, i.e., consider an arbitrary linear functional 

/:C 0 °°(fi) 9 *>.—></,¥>>€ C, (13.1) 

then we obtain an element of the space 2?#(fi), which will be called 
a generalized function (in the domain fl) . Let us give the exact 

13.1. Definition. T>#(VL) is the space of all linear functionals 
(13.1 1 in which the operations of differentiation d a and multipli- 
cation by a function a £ C°°(Q) are introduced by the following 
formulae: 

(d a f,p) = (-1 )'“'</, (af,<p) = (f,a<p) \/tp £ C£°(fi). 

(13.2) 



13.2. Example, f = Y ^ k \x - k), 
k = o 

oo 

£(-i)V fc) (fc) € cs°c 

k= 0 



X £ 



Obviously, / £ 



ie., (f,<p) = 
, and f i 
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2? b (]R). Thus, 2? b (fi) C D#(Q). By the way, the following lemma is 
valid. 

13.3. Lemma (P. du Bois Reimond). If f G 2?#(R) and f = 0, 
then f = const. (Thus, f G T ,b (K). / ) 

Proof. We have (f,p) = (/, p') = 0 \/p G Cq°(K). Let us 
take a function po G C“(M) such that J po = 1. Any function 
p G Co°(lR.) can be represented in the form p = p\ + (f p) po, where 
pi = p—(f p) Po- Note that f p i = 0. Setting ip(x) = pi(£)d£, 

we have ip G C“(M) and ip' = p\. Therefore, (f,p) = ( f,ip ') + 
(/, (f p) p 0 ). Since ( f,ip ') = 0, it follows that (f,p) = C f p, where 
C=(f,P o). ' □ 

Generalizing the notion of a 5-sequence, we introduce 

13.4. Definition. A sequence of functionals f v G T># is said 
weakly converges to / G T> # on the space 4> D Cg°, if /„ — > / in 
V # on the space 4>, i.e., lim^^ 00 (/ !/ , y>) = (/, p) \/p G 4>. If 4> = Gg”, 
then the words “on the space Gg°” are usually omitted. 

13.5. Definition. We say that a subspace X of the space V # is 
complete with respect to the weak convergence, if, for any sequence 

°f functionals /„Gl satisfying the condition 

(U - fu> P) — » 0 VpGC“ as yg — > oo, 
there exists / G X such that f v — > / in T>#. 

13.6. P. Show that V 7 is not complete with respect to the weak 
convergence. 

13. 7. P. Show that T> # is complete with respect to the weak con- 
vergence. 

13.8. Lemma. If /„ — > / in T> # on the space 4) D Cq°, then 
d a f v — > d a f in T># on the space 4> for any a. 

Proof. {d°f v ,p) = (-1)'“' (/„d» — ► (-l)H (f,d a p) = 
{d a f, p). □ 

13.9. Example. Let f v = i.e., (f„,p) = f K ^f^p(x)dx. 

Then f' v = cosvx, f" = -u-smux,... We have ( f v ,p ) — > 0 
W p G Cq° as v — » oo. Thus, cosvx —> 0 in T > #, v sini/x — + 0 in 
V#,... 
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13.10. Lemma. Let a = (a±, . . . , a n ) £ f 2 c R ra . Suppose that 
a sequence of functions /„ £ Lj oc (fl) and a point b = 

(&i, . . . , bn) 1 ') £ II, where 

II = {t = (*i, . . . ,x n ) £ K n | \x k - a k | < cr k , cr fc > 0 Vfc} C n, 

are such that, for F v (x) = f^ 1 . . . f^ n f v (y)dyi . . . dy n , the following 
two properties hold: 

P As b = (b\, . . . , b n ) one can take any point of II such that b k < a k 
Vk. 



(1) \F v (x)\ < G(x), x £ ft, where G £ Lj oc (Q), 

(2) F v (x) — > B+(a; — a) almost everywhere in fl, where 0+ is 



defined in (12.5). 



Then /„ weakly converges to S(x — a) on the space 

$ = {ip £ C(Sl) | ip £ d n p/d Xl . . . dx n £ L 1 ^)}. (13.3) 



PROOF. Using Theorems |8.20[ [8724] and |8.27| we obtain, for any 

<P € 



(fv,<P) = 



d n F„ 



dx\ . . . dx r 



<p) = (-l) n (F v , 



d n p 



dxi . . . dxr, 



OO OO 






d n ip{x)dx 
dx n 



OO OO 



= -(-ir 



d n 1 ip(x)dx 
dx 2 ■ ■ ■ dx n 



dx 2 . • . dx n = (p(a). 



□ 



13.11.P. Using Lemma 13.10 solve problems p\TI I and P\4-4\ 



Let us generalize the notion of the support of a function (see 
Section [3]), assuming, in particular, an exact meaning to the phrase 
usual for physicists: “6(x) = 0 for x 7^ 0”. 



13.12. Definition. Let f £ 2?#(f2), and u be an open set in 
O. We say that / is zero (vanishes) on u> (and write / 1 =0 or 

f(x) = 0 for x £ f2), if (/, ip) = 0 Wip £ Gq°(uj). 
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13.13. Definition. The annihilating set of a functional f £ 
D#(Ll) is the maximal open set Ll 0 = Ll 0 (f) C D on which / is zero, 
i.e., f\ n = 0, and the condition / = 0 implies u) C Do. 

It is clear that Do(/) is the union of w C Ll such that / 1 =0. 

13.14. Definition. Let / £ V#(Ll). The support of the func- 
tional /, denoted by supp /, is the completion to the annihilating set 
D 0 (/), i.e., the set Ll \ Ll 0 {f). 



13.15.P. Let f £ V#(Ll). Check that x £ supp / if and only if, 
for any neighbourhood uj C Ll of the point x there exists a function 
tp £ Cq°(uj) such that (f,tp) ^ 0. Verify also that Definition 13. 1J 
is equivalent to Definition 3. m if f£C{Ll). 



13.16.P. Find supp 8^ a \x) and supp[(:ri 



a; n )(5 (a) (a;)]. 



13.17.P. Let f £ V#{Ll), a £ C°° , and a{x) = 1 for x £ supp/. 
Is it true that a ■ f = f? 



13.18. P. Let uj be an open set in Ll such that oj D supp/, / € 
D#(Ll). Show that af = f, if a{x) = 1 for x £ to. 

13.19. P. Let f £ D#(Ll) be a generalized function with a com- 
pact support. Show that the formula (F,ip) = (f,ifip) \/<p £ C°°(Ll), 
if £ Cff°(Ll), where if = 1 on an open set co D supp / defines the ex- 
tension of the functional f onto the space C°°(Ll), i.e., F is a linear 
functional on C°°(Ll) such that (F,ip) = (/, ip) \/<p £ Cq°(LI). 



14. The problem of regularization 



The idea of representability of a function / : Ll — ► C with 
the help of its “averaging” functional (10.1 1 concerned only locally 
integrable functions. However, in many problems of analysis, an 
important role is played by functions which are not locally integrable. 
This is the reason of arising of the so-called problem of regularization: 
let g : LI 9 x <— > g{x) be a function locally integrable everywhere in 
Ll except a subset N C Ll. It is required to find functionals / £ V# 
such that 

(/,¥>) = J g{x)tp(x)dx \/ip £ C™(Ll\N). 



(14.1) 
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In this case one says that the functional / regularizes the (di- 
vergent) integral f n g(x)dx. 

It is clear that the functionals / satisfying (14.1 1 can be repre- 
sented in the form 



/ — fo + fi, /o G Fq j 

where /1 is a particular solution of the problem of regularization (i.e., 
fi satisfies ( |14.1 1), and Fo is the linear subspace of the functionals 
fo £ such that 



</o,¥>)=0 Vip £ C£°(n\N). 



(14.2) 



The question of description of the subspace F 0 is connected only 
with the set N D supp/o- In the case when N = Xq £ ft, this 
question, i.e., the problem concerning the general form of function- 
als with a point support, is considered in Section |15| As for the 
particular solution of the problem regularization, we conclude this 
section by consideration of the regularization for 1/P, where P is a 
polynomial in the variable x £ R. 



14.1. Example. Consider the regularization of the function 1/x. 
In other words, find the functional / £ P#(R) which satisfies the 
condition: x ■ f = 1. Note (see (12.6)) that 



OO 

up.-,<p\= [ -ip(x)dx \/y £ Cff 
x / J x 



(R\0). 



Thus, the functional v.p.(l/a;) regularizes the function 1/x. Since 
(5, tp) = 0\/ip £ Cq°(R\ 0), it follows that v.p.(l/x) +C ■ S(x), where 
C £ C; therefore, (see ( 12.7 1) functional l/(x±*0) also regularize 
the function 1/x. 



14. 2. P. Check that 



v-P-~,V ) = 
x 



ipfx) - y{-x) 
2x 



dx \/ip £ 



68 



2. THE SPACES T)\ T>* AND T>' 



14.3. P. Let m > 1, and a £ C^°(K). Define, for k > 2, the 
functional v.p.(l/x k ) £ 2?#(R) by the formulae: 



1 



OO 



¥>( 0 )- 



for k = 2m and 



y.k — 2 



(fc-2)! 



v? (fe " 2) (0) 



da; 



V{x) 



+ tp(—x) — 2 



x<//(0) 



„ fc -2 



(fc-2)! 



^ (fe - 2) (0) 



dx 



for k = 2 to + 1 . 

57iow that the functional v.p.(l/x k ) regularizes the function l/x k . 



14.4. P. (Compare with P 16.25). Find the solution f £ £> # (R) 



of the equation P{x)f = 1. In other words, regularize the integral 
P~ 1 (x)ip(x)dx, where P is a polynomial. 

15. Generalized functions with a point support. The Borel 

theorem 

It has been shown in Section [14] that the mean value problem for 
a function locally integrable everywhere in fl C R” except a point 
£ £ fl leads to the question on the general form of the functional 
f £ U#(fl) concentrated at the point £, i.e., satisfying the condition: 
supp / = f. It is clear (see P |13.16 ), that a finite sum of the 5- 
function and its derivatives concentrated at the point £, i.e., the 
sum 

c a 5 {a) {x-f), c a £C,N£N (15.1) 



E 

|a|<A 



is an example of such a functional. 



However, is the sum (15.1 1 the general form of a functional / £ 



D# , whose support is concentrated at the point £? One can show 
that the answer to this question is negative, however, the following 
theorem holds. 
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15.1. Theorem. If f £ T> # and f = Y c a^ a \x — £), then 

ot 

c a = 0 for |a| > Nf for some Nf. 

PROOF. According to the Borel theorem below, there exists a 
function <p G C™ (fi) such that, for any a 

d a v(x)\ x=l _ = (-l) |a| /ca, if c a ± 0 

and 

d a ip(x) 1^ = 0, if c a = 0. 

For such a function p, we have ( Y c a S^°‘\x — £), ip ) = Y 1, where 

\ ot / ot 

the sum is taken over a for which c a yf 0. □ 

15.2. Theorem (E. Borel). For any set of numbers a a € C, 

parametrized by the multiindices a = (ati, . . . , a n ), and for any point 
£ G hi C K n , there exists a function ip G such that d a tp\ x _^ = 

fl/Q \/ Ot . 

PROOF. Without loss of generality, we can assume that £ = 0 G 
0. If the coefficients a a grow not very fast as \a\ — + oo, more exactly, 
if there exist M > 0 and p > 0 such that Y a « < Mp~ k Mk G N, 

|q;|=/c 

then the existence of the function required is obvious. Actually, since 
in the case considered the series Y a a x a /a\, where a! = aq! a n \, 

a 

converges in the ball B p = {x G K" | \x\ < p}, we can take as the 
required function the following one 

p(x) = if(x/p)J2a a x a /a\ G C?{B P ) c Cg°(fi) % 

ot 

where 

^GC 0 °°(R n ), ^> = 0 for \x\ >l,ip=l for |*| < 1/2. 

However, in the general case, the series Y o-aX a /o\ can diverge 

Ot 

in B p . What is the reason of the divergence? Obviously, because it 
is impossible to guarantee the sufficiently fast decrease of a a x a /a! 
as | a | — > oo for all x belonging to a fixed ball B p . One can try to 
improve the situation, by considering the series 

TMx/Pa ) ' a a x a /a \ , 



(15.2) 
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where p a converges sufficiently fast to zero as |a| — ; 
that series (15.2 1 converges to a function ip £ C c 
Cq°(CI) and d a p\ 



oo. If it occurs 
then, as one 

can easily see, <p G Cq°(£1) and d a ip\ x 0 = a a . Indeed, setting 
7 = ( 71 , ■ ■ • , 7 n) < P = {Pi, ■ ■ ■ , Pn) by definition, if < f3 k Vfc, and 
/3 - 7 = (/3i - 71 , - 7 „), we have 



»>L.o = ( £ ^Y)'- (s-^) u (9T-’) |,. 0 j 

= Ylwm 
0 

= J2( a 0/P-) ( da x : ; | x=0 ) + o« = a«. 



It remains to show that series (15.21 converges to tp £ C°°(Q). Note 
that since Y = Y + ^ , it is sufficient to verify that there exist 

a |a|<fc |a|>fc 

numbers p a < 1 such that 



EE ip(x/p a )a a x a /a\ £ C k ( fi) Vfc. 

j>k \a\=j 

Let us try to find p a = pj depending only on j = |a|. If we can 
establish that V/3 such that |/3| < k, the following inequality holds 



\d^ {ij){x / p\ a \)a a x a / a\\ < C a p a , (15.3) 



where C a 
obtain 



C a (ip) < 00 j then, taking pj = 2 



—3 



Y c a 

J“l=j 



we 



E E |5 /3 (^(^/p| tt |)aaX“/a!| 

j>k \a\=j j>k 



pi E 

i«i=j 



< 1. 
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Thus, it remains to prove ( |15.3| ). For |a| > k > |/3|, we have 

/?! 






x \ a„x 



P\a\ 



a\ 



a n 



<- 

a\ 



E 



1<0 
|a„|x' v 



^ a! 



7 ! (/3 — t) ! 
(5\ 



dli> 



X 

P\a\ 



|9 /3_7 x“| 



E 



1 



M 



^ 7 !(/3-a)! \ P|q| , 
x|W)| t=x/pM -x“-^ 7 |.a 

5 t 7 ^W| t=x/P|o 



<K| E 



/?! 



T</3 



7 ! (/3- «)! 



P|a|- 



□ 



Now return to the question on the general form of the generalized 
function f £ V#(Q) with the support at the point £ = 0 £ Q. First 
of all, note (see P |13.19 ) that, for any function a £ Cq°(CI) such 
that a = 1 in a neighbourhood of the point £ = 0, the formula 
(f,<p) = (f,aip) V</3 £ C 00 (fl) is valid. In particular, the functional 
/ is defined on the polynomials. Setting c a = (— l)l“l (/, x Q /a!), we 
obtain 

(/,¥>)= E c <x{s {a \v) + (iinv) ViV, 

|a|<JV 

where 

y>(x) — ^ , 0 < ejv < 1- 

|a|<IV 

(15 - 4) 

It is rather tempting to assume that, for an appropriate sequence 
{ e -/v}? 7 =i> 0 < ejv < 1, the following condition holds: 

(/, r N ) — >0 as N — > 00 , (15.5) 

because in this case Theorem 1 1 5 . 1 1 implies the obvious 



rjv(x) = a — 



15.3. Proposition. If f £ supp / = 0 £ tt and ( |15.5[ ) 

is valid, then 3 N £ N such that f = c a 5^ a \ 

\a\<N 



However, in general, condition (15.5) does not hold, if / £ V #. 
An appropriate example can be constructed with the help of so-called 
Hamel basis (see, for instance, [136]). 
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16. The space V' of generalized functions (distributions by 

L. Schwartz) 

The wish seems natural to have a theory of generalized functions 



in which condition (15.51 is satisfied, hence, Proposition 15.3 holds. 



This modest wish (leading, as one can see below, to the theory of 
the Schwartz distributions) suggests the following: 

1) introduce a convergence in the space Cff (Ll) such that for 
this convergence 

lim rjv = 0 £ (16.1) 

N — »oo 



where r n is defined in (15.4) for some suitable e n £]0,1[; 

2) consider below only the functionals / £ D#(fl), which are 
continuous with respect to the convergence introduced. 

It is clear that one can introduce different convergences accord- 
ing to which rjv— > 0 as Af - >oo. Which one should be chosen? 
Considering this question, one should take into account that the 
choice of one or another convergence also determines the subspace 
of linear functionals on Cfi° that are continuous with respect to this 
convergence. Therefore, it seems advisable to add to items 1) and 
2) above the following requirement: 

3) the space of functionals continuous with respect to the 
convergence introduced must include the space V a of the 
Sobolev derivatives (since this space, as has been shown, 
plays very important role in the problems of mathematical 
physics). 



According to Theorem 16.1 below, requirement 3) uniquely de- 

moreover, (see Proposi- 



termines the convergence in the space Cq°; 
tion 16.101 condition (16. 1[) is also satisfied. 



O 

°0 



16.1. Theorem. Let {<fj} be a sequence of functions ipj £ 
Then the following two conditions are equivalent: 

(/, Tj) -> 0 as j -> oo V/ G V b (Ll); 
a) there exists a compact K C LI such that supp tpj C K 

Vj; 

0 as j — > oo Va. 



(fi)- 

l c 

2 C 



b) max \d a tpj{x)\ 

x£il 



PROOF. The implication 2° ==> 1° is obvious, 
assertion follows from Lemmas 116.21 116.51 below. 



The converse 
□ 
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16.2. Lemma. Va 3 C a such that max \ip^\x)\ < C a Vj. 

Proof. For any a, consider the sequence of functionals 

<^ q) : I^fi) 9 / i — >■ J f(x)d a ipj(x)dx, j > 1, 
n 

defined on the space L 1 ( fl). The functionals are, obviously, 
linear and continuous, i.e. , (by the Riesz theorem 9.141 G L°° . 
According to condition 1°, /) — > 0 as j — > oo V/ € L 1 . There- 

fore, by virtue of the Banach-Steinhaus theorem 1 ^ there exists a 
constant C a such that 11^4°^ ||oo < C a Vj. □ 



11 The Banach-Steinhaus theorems (1927) asserts the following 
(see, for instance, m or )56| ). Let X be a Banach space and {pj} be a 
family of linear continuous functionals on X. If for any x G X there exists 
C x < oo such that x)| < Cj Vj, then there exists a constant C < oo 
such that \{ipj,x)\ < C for ||a:|| < 1 and Vj. 

Proof. Suppose the contrary be true and note that if a sequence of 
functionals (pj is not bounded for ||a:|| < 1, then it is not also bounded 
in the ball B r (a) = {x G X | ||a; — a|| < r}. Let us take a point xi G 
B i(0), a functional ifik! and a number n < 1 such that |(y>*, 1 ,a;)| > 1 for 
x G B ri (xi) C -Bi(O). Then we take a point X 2 G B ri {x i), a functional 
ifik 2 and a number r 2 < n such that | {ip k2 , x) \ > 2 for a: G B r2 (x) C 
B ri (x i). Continuing this construction, we obtain a sequence of closed 
balls B rk (xk) embedded in each other, whose radii tend to zero. In this 
case, \(pkj ,*o)| > j for xo G C\B rk (the intersection C\B rk is non-empty 
by virtue of completeness of A). □ 



16.3. Lemma. Va Vxq G d a ipj(x o) — > 0 as j — > oo. 

Proof. d a tpj( xo) = (S^ a \x — ieo), (— 1 )^<£ j (®)) — > 0 , since 

6( a \x- x 0 )£V\n). □ 

16.4. Lemma. There exists a compact K c such that for all 
j we have supp^- C K. 



74 



2. THE SPACES D\ D* AND T>' 



PROOF. Suppose the contrary is true. Let Kj = (J supptpfc. 

k<j 

We can assume that the intersection (supp^j) D (O \ K) is non- 
empty, i.e. , 3 Xj £ Q\Kj such that <Pj(xj) yA 0. For any j, we choose 
A j > 0 such that 

> \ Vx£Vj = {|® - Xj\ < A,} c Mj = slippy \ Kj. 

(16.2) 

Note that Vj ft Vk is empty for j yA k and consider the function 
/ £ L\ oc { fi) which is equal to zero outside U 3 >i Vj and such that 

f(x) = aj\ipj(xj )\~ 1 exp [— i arg ^Pj(x)\ for x £ Vj, j> 1, (16.3) 

where aj > 0 are constants which will be chosen such that we obtain 
an inequality contradicting 2 ^ inequality 



fiPjdx 



>3- 



(16.4) 



Note that supp f<pj C (Vj U Kj), because supp^j C (Mj U Kj). 
Therefore, the last integral in the equality 



/ 



ftfij dx 



ipj dx + 



f<Pj dx 



n Vj (supp f<Pj)\Vj 



can be estimated by f K f^fjdx < max|</5j| f K \f\dx < Aj , where 

Aj = C J2 l a fc| ' A t(Vfc) . Taking a 7 - = 2 (Aj + j), we obtain (|16.4|) , 
k<j 

since, by virtue of (16. — ( |16.3| ) , j v f^Pjdx > a.j/2. □ 



Inequality ( 16.4 1 contradicts the initial condition 1° 



16.5. Lemma. For any a , e > 0, a,’o £ ft there exist A, v > 1 
such that |<^.°^(a;)| < e for \x — cco| < A and j > v . 

PROOF. Suppose the contrary is true. Then 3cr 3e 0 > 0 3xo £ ft 
such that, for any j, 3 Xj £ {x £ ft | \x — Xo| < 1 /j} such that the 
inequality \tp^\xj)\ > Cq holds. However, on the other hand, 

\^ a) (xj)\ < I <Pj a) (xj) - (Pj a) (x o)| + \v { 3 a) (x 0 )\ — > 0, 
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because 

\<Pj a \xj) - ip ( j a \x 0 )\ < C\xj - a" 0 | -*• 0, and d a tpj(x 0 ) -»■ 0, 
according to Lemmas |16.2| and |16.3| □ 

16.6. Remark. Actually, we have proved a bit more than has 
been stated in Theorem |16.1| Namely, condition 2° follows from the 
proposition: (f,(pj) — > 0 as j — > oo for any f £ L\ oc (Cl) and for any 
Sobolev derivative / = d a g, where g £ L 1 ( fi). 

Now we can define the spaces V and V introduced by L. Schwartz 

m- 



16.7. Definition. The space T>(Cl), which is sometimes called 
the space of test functions (compare with Section [TJ is the space 
Co°(f2) in which the following convergence of sequence of functions 
ipj £ Cfi°( ft) to a function p £ Cq°(CI) is introduced: 

a) there exists a compact K such that supp^- C K Vj; 

b) \/p = (Pi, , p n ) Vcr > 0 3 N = N(P, a) £ N such that 

| d /3 ipj(x) — d / V(:r)| < a Vx £ Cl for j > N. 

In this case we write ipj — > in T> as j — ► oo (or limj^oo ipj = ip in 

V). 



16.8. Remark. It is clear that V(Cl) = f] T> s (Cl), where T> s (Cl) 

s>0 

is the space of functions Cq(D) equipped with a convergence which 
differs from the one introduced in Definition 16. 7| only by the fact 
that the muliindex in condit ion b ) satisfies the condition \P\ < s. 
One can show (see P 16.23 1 that T> b (Cl) = (J T>' S (CY) (i.e., / £ 

s>0 

V b 3s > 0 such that / £ V s ), where V s (Cl) is the space 

of linear functional on T> S (CV) continuous with respect to this con- 
vergence in T> s (Cl). The spaces T> s and T>' s have been introduced by 
Sobolev 



16.9. Definition. The space T>'(CT) of the Schwartz distribu- 
tions (called also the space of generalized functions) is the space of 
linear continuous functional on T>(Cl), i.e., of linear functionals on 
D(Cl) which are continuous in the convergence introduced in T>(Cl). 



16.10. Proposition. There exists a sequence {eiv}jv = i> 0 < 
ejv < 1, such that lim tat = 0 in V, where vn is defined in (15.4 1. 

N—>oo 
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PROOF. By the Teylor formula, 

1 

rj\ r(x) = a 'j /(l — t)^ (p( a \tx)dt. 

' £N ' l«|=JV+i a ' { 

This and the Leibniz formula imply that {d^r^lx)] < CnUn) N ~^^ < 
(1/2) N / 2 for N>N 0 = 2|/3|, if £jv < iC” 2/JV . □ 

Propositions |15.3| and |16.10| imply 

16.11. Theorem (L. Schwartz). If f £ T>'(n), supp / = 0 £ SI, 
then there exist N £ N and c a £ C such that f = ^ CaS^ . 

\a\<N 



16.12.P. Let fk £ 2T(M), where k = 0 or k = 1, and x- fk(x) = k. 
Show (compare with Example 1^.1) that fo(x) = CS(x), fi(x) = 
v.p .L + C6(x), where C £ C. 



The following series of exercises P |16.13| -P |16.25| concerns the 
question on the structure (general form) of distributions. Some hints 
are given at the end of the section. 



16.13. P. Verify that the following assertions are equivalent: 

a) f is a distribution with a compact support, i.e., f £ V(£i) 
and supp f is a compact in S2 : 

b) f £ £'{VL), i.e., f is a linear continuous functional on £ (S2), 
i.e., on the space C°°(£l) with the following convergence: 
lim ipj = ip in £ •<=>■ lim atpj = atp in V Va £ C£°(fl). 

j — KX> j — >00 

16.14. P. Prove that f £ T>' (Vl) if and only if f £ and 

for any compact K C S2 there exist constants C = C{K , /) > 0 and 
N = N(K , /) £ N such that 



I (/,¥>) I < C-p KiN (ip) 

Vip £ C™ (K, Q) = {if£ Cg°(ft) | supp C K}, (16.5) 

where 

Pk,n{<p) = E sup \d a <p(x)\. (16.6) 



|a|<JV 



xGK 



16.15. P. ( Compare with P 16.14). Let Un>iKn = SI, where 
Kpj are compacts in R n . Show that f £ £'{{ 2) (see P 16.13) if and 
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only if f £ and there exist constants C = C(f) > 0 and 

N = N(f) > 1 such that \(f,p)\ < C ■ Pn{p) W> £ C£°(fl) ; where 

Pn(p) = E sup \d a p(x)\. (16.7) 

\u\<N x&Kn 



16.16.P. (Continuation). Let f £ £'(fi), supp / CwgOc 
Using ( 16.7 1 and noting that |V’( a; )l < Jq d 



-4>{x) 



dx Uip £ 



dx± ...dx n 

C^°(il), show that there exist numbers C > 0 and k > 1 such that 

Qnm 



\(f,P)\<C 



n 



dx r { 1 ...dx\ 



;P(x) 



dx Up £ Cg°(uj). (16.8) 



16.17.P. ( Continuation ). Checking that the function p £ 
can be uniquely recovered from its derivative ip = ax ,f dx „ 2 p, show 
that the linear functional l : ip > (/, <p) defined on the subspace 
Y = {ip £ C 0 (u>) | ip = dx £"™ dx rn P, P £ Cg 0 )} of the space L^(uj) is 
continuous. 



16.18. P. (Continuation). Applying the Hahn-Banach theorem 
on continuation of linear continuous functionals (see, for instance, 
|36j ). show that there exists a function g £ L°°(uj) such that 

/ PlTlTTl 

9 ^ dx?...dxz ,p ^ dx = e c °°^' 

c o 

16.19.P. ( Continuation). Show that the following theorems hold. 

16.20. Theorem, (on the general form of distributions with 
a compact support). Let f £ £'(f l). Then there exist a function 
F £ C(Q) and a number M > 0 such that f = d a F, where a = 
(M, ■ ■ ■ ,M), i.e., 

(/,¥>> = (-i) w J F(x)d a <p( x )dx up £ qj°(n). 

a 

16.21. Theorem (on the general form of distributions). Let f £ 
Then there exists a sequence of functions F a £ C(fl), para- 
metrized by multiindices a £ Z™, such that f = y) d a F a . More 

Ot 

oo 

exactly, F a = Y) F aj , F aj £ C(fl), and 
j = i 
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(1) suppFo,. C f lj, where {0?}j>i is a locally finite cover of 

( ->' 

(2) Vjf > 1 3Mj > 1 such that F a . = 0 for |a| > Mj. 



16.22.P. (Peetre |47j ). Lei A : £(fl) — > £(fl) be a linear continu- 
ous operator with the localization property, i.e., supp(Au) C supp(it) 
Mu £ £(f 1). Then A is a differential operator, more exactly: there 
exists a family {a Q } ae z™ of functions a a £ C°°(fi) such that, for 
u £ (. Au)(x ) = a a (x)d a u(x), where m(x) < N(K ) < 

|a|<m(x) 

cx) for any compact K C fl. 



16.23. P. 



(See Remark 



16.8). 



Check that = |jL>'(fi). 

S 



16.24. Remark. We say that a functional f £ V has a finite 
order of singularity, if there exist k > 1 and functions f a £ L\ oc , 
where |a| < k, such that / = ^ d a f a . The least k , for which such 

|q:|=/c 

a representation of / is possible, is called its order of singularity. 
In these terms, the space of Sobolev derivatives V b is, according 
to Definition |12.2| the space of all distributions which have a finite 
order of singularity. 



16.25. P. Resolve the following paradox. On one hand, the dis- 
continuous function 

f3?exp(— 1/z 4 ) for z / 0, z = x + iy£ C, 

|^0 for z = 0 

(being the real part of a function analytic in C \ 0 with zero second 
derivative with respect to x and y at the origin) is a solution of 
the Laplace equation on the plane. On the other hand, by virtue of 



f{x,y) = 



(16.9) 



Theorem 16.20 and a priori estimate (21.6) (see also Corollary 22.24 
and [22l no 2, § 3, item 6]/' if f £ V {A) and A / = 0 in LI, then 
f £ C°°(A). 



16.26.P. Show that S is metrizable and T> is non-metrizable. 



16.27. Remark. One can introduce in T> (respectively, in E) a 
structure of so-called (see |36L 51 i) linear locally convex topological 
space (LLCPS 3 ) and such that the convergence in this space coincides 
with the one introduced above. For instance, the neighbourhood of 
zero in T> can be defined with the help of any finite set of everywhere 
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positive functions 7 m £ C(fi) (m = 0, 1, . . . , M; M £ Z+) as the set 
of all the functions p £ Cg°(fi) such that \d a p\ < 7 | a |, if |a| < M. 
The topology in £ can be defined by introducing the distance by the 
formulae given in the hint to P |16.26| Thus, £ is the Frechet space , 
i.e. , a complete metric LLCTS. One can extend to the Frechet spaces 
(see, for instance, ( 511 ) the Banach-Steinhaus theorem: a space of 
linear continuous functionals on a Frechet space (in particular, the 
space £') is complete relatively the weak convergence. Although V 
is not a Frechet space (by virtue of P 16.261, V is also complete 
relatively the weak convergence (the direct proof see, for instance, 



3) A linear space X is called a linear locally convex topological space, if 
this space is topological [361 , the operations of addition and multiplication 
by a number are continuous and, moreover, any neighbourhood of zero in 
X contains a CONVEX neighbourhood of zero. 

Hints to P fl6T3h P ll6^6l 



16.13 If we suppose that b) does not imply a), then there exists 

dfl, and f ^ 0 in the vicinity 



a sequence of points Xu such that Xk 
of Xk- 



16.14 



If f £ V but estimate (16.51 is not valid, then 3 1\ = 
K c O \/N > 1 3p £ suppp C K n , and |(/, <^iv)| > 

N sup|y>j(^|. We have = Pn ■ \(f, — > 0 in V but 

|a|<AT K 

\(f,}M\= i- 



16.15 Since (/, tp) = {f,pip), where p £ Cfi°, p = 1 on supp /, 



it follows that K = supp p. Warning: in general case, K ^ supp /. 
Indeed, following t54 i vol. 1, p. 94], consider the functional / £ £'(R) 
defined by the formula 



(/, ip) = lim 



“ mp(0) - (lnm)^'(O) 



l v<.m 



Obviously, supp / is the set of the points of the form 1/u, v > 1, 
together with their limit point x = 0. Consider the sequence of 
functions pj £ Cq°(R) such that Pj(x ) = 0 for x < ttj, Pj{x) = 



for 1/j < x < 1. Taking K = supp / in (16.6 1, we have 
Pk,n{Po) -> 0 as j -> oc VA r > 1, while (/, p) = j/y/j -> 00 . 
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16.16 



16.17 



16.18 



sup|9“(^(x)| < Cj( y K)s\vp\-^d a (p{x)\ 



Apply (16.8). 



By the Riesz theorem (see Theorem 9.141, (L 1 )' = L°° . 
Complete the definition of g outside ix by zero (see 



16.20 

16.18| and take F(x) = (— 1) 



J y<x g(y) dy. 



16.21[ Let — 1 be a partition of unity. We have 



(/,¥>) = TMfM = 






3 M <m 3 



(d a F aj , ip) — ^(<9“ ^2 Fa i > ¥’)• 



16.22 



Verify that the functional A a : £ 9 u i— > Au\ x a belongs 
to £' and suppA a = a. Thus, Au | = Y2 (a a {x)d a u) | 

|a|<m(a) 

Using the Banach-Steinhaus theorem (see note 1 in Section 16 1, 



prove that sup m(a ) < oo for any K C K C Cl. Applying A to 
aeif 

(y — x) a /a \ , show that a £ C' 00 (U). 

Apply Theorem 16.21| 

Function ( 16.9 ) does not belong to D’ (i.e., does not admit 



16.23 



16.25 



the regularization in V) as well as any other function f £ C° 
which for any m £ N and C > 0 does not satisfies the estimate 
\f( x )\ < C|x| -m for 0 < |x| < e, where 1/e 1. The last fact 

can be proved, by constructing a sequence of numbers Cj > 0 such 
that , for the function c Pj(x ) = ejip(jx), where ip £ Co(K), <p = 0 
outside the domain {l<|x|<4}, J ip = 1, the following conditions 
are satisfied: J R „ f(x)ipj(x)dx — > oo as j — »■ oo but <pj — > 0 in V as 
j -> oo- 



16.26 The distance in £ can be given by the formula p^ip^) = 

OO 

d(tp — ip), where d(ip) = Y2 min(pjv (<p), 1), and Pn is defined 



in ( |16.7 1. D is non-merizable since for the sequence tpk,m{x ) = 
(p(x/m)/k , where ip £ 75(151), the following property, which is valid 
in any metric space, does not hold: if p>k,m — > 0 as k — > oo, then Vm 



3k(m) such that <pu r 



0 as to 



oo. 



CHAPTER 3 



The spaces H s . Pseudodifferential 
operators 



17. The Fourier series and the Fourier transform. The 

spaces S and S' 



In 1807, Jean Fourier said his word in the famous discussion (go- 
ing from the beginning of 18th century) on the sounding string |50j . 
Luzin wrote [42] that he accomplished a discovery which “made a 
great perplexity and confusion among all the mathematicians. It 
turned over all the notions” and became a source of new deep ideas 
for development of the concepts of a function, an integral, a trigono- 
metric series and so on. Fourier’s discovery (however strange it seems 
at first sight) consists in the formal rule of calculations of the coef- 
ficients 



«fc 



p/2 




u(y)e-° l(k/p)y dy, 



\ = 2ni , i = \/-T 



(17.1) 



(which are called the Fourier coefficients) in the “expansion” 

OO 

u(x) ~ ^2 a, k e ' l(k/p)x , \x\<p/2 (17.2) 

k — — oo 



of an “arbitrary” function u : fi 9 x i — > u{x) £ C, where fi = 
] — p/2,p/2[, by the harmonics 



e i(k/ P ) x _ cog 2n(k/p)x + isin2n(k/p)x, k £ Z. (17.3) 



The trigonometric series ( 17.2 1 is called the Fourier series of the 
function u (more exactly, the Fourier series of the function u in the 
system of functions ( 17. 3| ) ) . The first result concerning the conver- 
gence of the Fourier series was obtained by 24 years old L. Dirichlet 
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in 1829 (see, for instance, [ 72 ]): if u is piecewise continuous on 
[—p/2, p/2] and the number of its intervals of monotonicity is finite, 
then the Fourier series of the function u converges to u at every 
point of continuity of it; moreover, if the function / is continuous 



and u(—p/2) = u{p/2), then series (17.21 converges to u uniformly. 



Although the Fourier coefficients (17.1) are defined for any function 



u £ L 1 , the Fourier series can diverge at some points even for contin- 
uous functions (see, for instance, 



also compare with P 17.9 1 . 



As for integrable functions, in 1922, 19 years old A.N. Kolmogorov 
constructed m the famous example of a function u £ L 1 , whose 
Fourier series diverges almost everywhere (and later an example of 
the Fourier series, everywhere diverging, of an integrable function). 



See in this connection formulae (|17.17|H 17. 



The following theorem (see, for instance, |56| 1 on convergence 
of the Fourier series in the space L 2 is of great importance: “for any 
u £ where S2 =] — p/2,p/2[, series (17.21 converges to it in 

L 2 (Siy\ i.e. 



\U- X! a k e k\\L 2 



0 as N 



oo, 



where 



|fc|<JV 



efc : 9 x 



e k (x) = exp ( i—x 

p 



(17.4) 



(17.5) 



This theorem shows the transparent geometric meaning of the Fourier 
coefficients. Actually, consider in L 2 (Q) x L 2 (Q) the complex- valued 
functional 

p/2 



(it, v) = j u(x)v(x)dx, 



-p/2 

where v is the complex-conjugate function to v. Obviously, this 
functional defines a scalar product ^ in the space L 2 (fl) with respect 
to which (one can readily see) functions (17.5) are orthogonal, i.e., 

(efc, e m ) = 0 for k ^ m. (17-6) 



Therefore, by choosing N > |m| and multiplying scalarly the func- 
tion (u - a k e k ) by e m , we obtain |(u, e m ) - a m ■ (e m ,e m ) | = 
|fc|<iV 
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|(u- E a k e k ,em)\ < \\u- E a, k e k \\ L 2 ||e m || L 2 . This and (|17.4|) 

|fc|<_/V |fc|<JV 

imply that 

dm ~ (^b ^m) / (f-m: ^ m ); (17.7) 



Thus, the coefficient a k is the algebraic value of the orthogonal pro- 
jection of the vector u £ L 2 (Q) onto the direction of the vector e k . 



2 ^ This means that the functional ( u , v ) is linear in the first argument 
and (u, u) > 0, if u yl 0, and (u, v) = ( v,u ), where the bar denotes the 
complex conjugation. Note that the function u i — *• ||u|| = y/ ( u , u ) is a 
norm (see note 5 in Section[8|, and |(n, u)| < ||w|| • ||v|| (compare with ( |9.3 1 
for p — 2). Recall also that the Banach space X (see note 5 in Section [8 1 
with a norm || ■ || is called a Hilbert space , if in X there exists a scalar 
product (•,•) such that ( x,x ) = ||a;|| 2 Va; £ X. Thus, L 2 ( fl) is a Hilbert 
space. 



When the geometric meaning of the Fourier coefficients became 
clear, it might seem surprising that, as Luzin wrote, “neither sub- 
tle analytical intellect of d’Alembert nor creative efforts of Euler, 
D. Bernoulli and Lagrange can solve this most difficult problem 3 ^, 
i.e., the problem concerning the formulae for the coefficients a k in 
(17.21. However, one should not forget that the geometric trans- 
parency of formulae (17.71 given above became possible only thanks 
to the fact that the Fourier formulae (17.1 1 put on the agenda the 
problems whose solution allowed to give an exact meaning to the 
words such as “function” , “representation of a function by a trigono- 
metric series” and many, many others. 



3) The reason of arising this question is historically connected with 
the problem of a sounding string (see 1421 l50l j — the first system with 
an infinite number of degrees of freedom which was mathematically in- 
vestigated. As far back as in 1753, D. Bernoulli came to the conclusion 
that the most general motion of a string can be obtained by summing the 
principal oscillations. In other words, the general solution u = u(x, t) of 
the differential equation of a string 



Utt - u xx = 0, |x| < p/2, t > 0 (17-8) 



which satisfies, for instance, the periodicity condition 

u(-p/2, t) - u(p/ 2, t) = 0, u x (-p/ 2, t ) - u x {p/ 2, t) = 0 (17.9) 
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can be represented in the form of a sum of harmonics propagating to the 
right and to the left (along the characteristics x ± t = 0, compare with 
Section 111, more exactly: 






a + e iX k(x +t ) 



+ a k e 






(17.10) 



fee z 



where a k £ C, and Afc = 2n k/p. Indeed, equation ( 17.8 1 and the boundary 



conditions (17.91 are linear and homogeneous. Therefore, a linear com- 



bination of functions satisfying ( 17.8 1 (17.91 satisfies these conditions as 
well. This fact suggests an idea to find the general solution of problem 



(17.81 (17.91, by summing (with indeterminate coefficients) the particu- 
lar solutions of equation ( 17.81) , which satisfy the periodicity conditions 



(17.91. Equation (17.81 belongs to those which have an infinite series of 
particular solutions with separated variables, i.e. , non-zero solutions of the 



form p(x)ip(t). Actually, substituting this function into (17.81, we obtain 
(x)i/j(t) = p(x)i/j t t(t). Hence, 






Pxx(x)/p{x) = = const . 



(17.11) 



The periodicity condition (17.91 implies that ip £ X, where 

I = {^C 2 (fl)nC 1 (Q) I *>(-§,*)«*¥>(§ ,*), ¥>'(-$,*) =¥»'( §,*)}■ 

(17.12) 

Thus, the function p necessarily (see |l7 . 1 1| ) ) must be an eigenfunction of 
the operator 



- d 2 /dx 2 : X 



L 2 (Q), n=]-p/2,p/2[. 



(17.13) 



This means that p is a non-zero function of the space A', satisfying the 
condition 

— d 2 p/dx 2 = p ■ ip, (17.14) 

for a constant p £ C, that is called the eigenvalue of operator ( |17. 13[ > . 
Since p belongs to the space A', it follows that the number p can be only 



positive (because otherwise ip = 0). Let us denote p by A 2 . Then (17.141 
implies that p(x) = ae lXx , A £ R, a £ C \ {0}. Obviously, this formula is 
consistent with the condition p £ X if and only if A = Afc = 2kTv/p, k £ Z. 
Thus, taking into account ( 17.11 1, we obtain 



/ \ / /.\ + i\ux iXut , — iXux 

Pk{x)ipk{t) = a k e fc e +a k e e 



— iXf^x — iXfct 



at e C, 



and, hence, D. Bernoulli’s formula (17.101 is established. 



The Bernoulli formula brought into use the principle of composition 
of oscillations as well as put many serious mathematical problems. One of 



them is connected with the finding of the coefficients a k in formula (17.101 



17. THE FOURIER SERIES AND THE FOURIER TRANSFORM 



85 



for any specific oscillation that (compare with P 11.21 1 is determined by 
the initial conditions 



u(x,0) = f(x), u t (x, 0) = g{x), 



(17.15) 



i.e., by the initial deviation of the string from the equilibrium position 
and by the initial velocity of the motion of its points. In other words, the 
Bernoulli formula posed the question of finding the coefficients a * from 
the condition 



— \ i\v.x 

i . \o K — 



g{x)- 



+ a k y Xkx = /(*)’ - 1 

fcEZ 

It is interesting that in 1759, i.e., in six years after the work by 



D. Bernoulli, formulae (17.1 1 answering the formulated question were al- 
most found by 23-year Lagrange. It remained to him only to make in 
his investigations the rearrangement of the limits in order to obtain these 
formulae. However, as Luzin writes, “Lagrange’s thought was directed in 
another way and he, almost touching the discovery, so little realized it 
that he flung about D. Bernoulli the remark “It is disappointing that such 
a witty theory is inconsistent.” 

As has been said, half a century later the answer to this question was 



given by Fourier who wrote formulae ( 17.1 1. This is the reason why the 



method, whose scheme was presented on the example of solution of prob- 



lem ( 17.8 H 17.9 1, (17.151, is called the Fourier method (see, for instance, 
The term is also used (by obvious reasons) the method of sep- 
aration of variables. This method is rather widespread in mathematical 
physics. 

The reader can easily find by the Fourier method the solution of the 
Dirichlet problem for the Laplace equation in a rectangle, by considering 
preliminary the special case: 

(x,y) G]0, 1[ 2 ; u\ x=0 = u\ x=1 = 0; 

Ui = 9 ^)- 

Equally easy one can obtain by the Fourier method the solution 



A u(x,y) = 0, 



u\ y=0 = f(x), 



u(x,t ) = UN(x,t) + ^2 



2 sin A k 



Afc[l + a sin 2 A*,] 



' cos A kX, N > 0, «o = 0 
(17.16) 



of problem ( 6.11 1 for the heat equation. In formula ( 17.16 1, is the fcth 



{k € N) eigenvalue of the operator 

—d 2 /dx 2 : Y — >L 2 (Q), ST =]0, i[, 
defined on the space 

y = W g c 2 (fi) n c 1 (si) | (*> ± <V)L =±1 = o}, 
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where a > 0 is the parameter of problem (6.111. Note that A k G](fc — 
l)n,(k — l/2)-7r] is the fcth root of the equation cotan A = a A and the 
eigenfunctions ipk(x) = cos A kX of operator (17.161 satisfy (compare with 
(17.61) the orthogonality condition: 



1 

/ 



(<Pk,<Pm) = / ifik{x)ifim{x)dx = 0 for k^m. 



(17.17) 



Indeed, integrating by parts (or applying the Ostrogradsky-Gauss for- 
mula in the multidimensional case) and taking into account the boundary 
conditions (ip ± ap')\ x _ ±1 = 0 and the fact that — tp'f. = A we have 



C 

(A k A m )(</5fe, <p m ) — I [}PkPm PmPk) dx 



-1 



<Pk<P' 



'll f ' ' ' I 1 f ' ' 

‘m |_j — / VkVm ~ + / l Pk l fim 



= 0 . 



One can show (see, for instance, |68j ) that the eigenfunctions ifk, k G N, 
form (compare with ( 17.4 1) a complete system in L 2 = Y, i.e. , Vu G 
L 2 Ve > 0 there exist N > 1 and numbers oi, . . . , oj v such that ||u — 

N 



< e - Therefore, (compare with (17.21 ( 17.7 1) the formal 

k= 1 



E 



{u,lpk) 
\ {Vk,Vk) 



V>k, 



(17.18) 



which is called by the Fourier series of the function u in the orthogonal 
(by virtue of (17.171) system of functions ipk, converges to u in L 2 . The 



reader can easily verify that (1, <fik)/(<Pk, <fik) = 2 sin A*,/(Afe[l + crsin 2 A*]) 
as well as that series |I7T6| ) converges uniformly together with all its 
derivatives for t > e for any e > 0 and determines a smooth (except the 
angle points ( x,t ) = (±1,0)) and unique (see, for instance, [ 20 ] ) solution 
of problem (6.11 1. 



Let us note one more circumstance. Series (17.161 converges quickly 
for large t. One can show that, for any k > 1, 



\u(x,t) — Mjv(x,t)| < 10 k /N fort>k/(4:.3N 2 ). 



(17.19) 



However, for small t, series ( 17.16| ) converges very slowly. Therefore, for 
small t, it is advisable to use another representation of the solution of 
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problem (6.11 1 that is obtained below in Section 18 with the help of con- 



cept of dimensionality (see Section[6| and the so-called Laplace transform. 



Substituting formally (17.1 1 into (17.21, we obtain 

p/2 

oo i 

u( x)= V -e^ k/p)x / e-° l( ~ k/p)v u{y)dy. (17.20) 

„ P J 



k—— 



-p/2 



Tending p to infinity and passing formally in ( 17.20 ) to the limit, 
we obtain, as a result, for the “arbitrary” function u : R — > C the 
(formal) expression 

OO OO 

u(x ) = J e lx t( J e~ lv ^u(y)dy^jdf. (17.21) 

— OO — OO 

From these formal calculations we may give the exact 

n 

17.1. Definition. Let £ € R n , x £ 1", x£ = x k£k, i.e., 

fe = i 

a:£ = ( x , £) is the scalar product of x and £. The function 

m(£) = J e~ zx ^u(x)dx, i = 2m, i = \f — T (17.22) 

R" 

is called the Fourier transform of the function u £ L 1 (K rl ), and the 
mapping F : L 1 (R") 9 u i — > u = F u £ C is called the Fourier 
transformation (in L^R”)). 

17.2. Lemma. If u £ L 1 (R"), then F u £ C(R ra ) and ||Fn||c < 

hiu,. 



PROOF. It follows from the (Lebesgue) theorem 8.20 that u = 
F u £ C(R); moreover, |zt(£)| < f R „ |u(x)| dx. □ 

17.3. Example. Let u±( x) = 0±(a:)e :Fax , where x £ R, a > 0, 



and 0± is defined in (12.51. Then u ± (£ ) = Let us note that 

^±2^ 

u± L 1 although u± £ L 1 . Note as well that the function u± can 
be analytically continued into the complex half-plane C=p 



Below, in Theorem 17.6 we give some conditions under which 



expression (17.211 acquires the exact meaning of one of the most 



important formulae in analysis. Preliminary, we give 
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17.4. Definition. Let p > 1 and k £ Z. We say that a function 
u £ L p (fl ) belongs to the Sobolev space W p,k (Ll), if all its Sobolev 
derivatives d a u, where |a| < k, belong to L P (H). The convergence 
in the space W p,k is characterized by the norm 

Mwv.> = E ( 17 - 23 ) 

\oi\ <k 

i.e., Uj — > u in W p,k as j — > oo, if ||u — uj\\ WP ,k — > 0 as j — > oo. 

One can readily verify that IT p,fc is a Banach space. 

17.5. Lemma. 4 ) W 1, "(K") c C'(K ra ), i.e., for any element 
{u} € W 1,n there exists a unique function u £ C which coincides 
almost everywhere with any representative 5 ' of the element juj, and 

IMIc<Nki.». 



Lemma 17.5 is a simple special case of the Sobolev embedding 
theorem (see, for instance, IS USES! HOI)- Note that the embedding 
W p,fc (M") C (7(1®™), valid for n/p < k, does not hold if p > 1 and n/p = k 



(see, in particular, Section |20[ where the special case p = 2 is considered). 



5 * See note 1 in Section [h] 



PROOF. It follows from the (Fubini) theorem 8.24 and Theo- 



rem 



8.27 that the function u admits the representation in the form 

dyi, 



Xl 

f 


X2 

f 


1 

*3 

CD 

3 

ST 


■ ' y dy 


...dy 2 


J 

-oo 


J ■ 

— OO 


J dy 1 dy 2 ■ 

— OO 


• ^Un 



x= (xi, . . . x n ), 

that implies its continuity and the estimate ||tt||c < / I I- D 
17.6. Theorem. Let u £ W 1,Tl (M"). Then, for any x £ K", 



N N 



u(x) = lim un(x), where ujv(x) = / ... / e. lx ^u{f)d^i . . . df n , 
N->oo J J 

-N -N 

(17.24) 

and u = Fu is the Fourier transform of the ( continuous 6 ) ) function 



u. 
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By virtue of Lemma 17.5 



PROOF. It follows from the Fubini theorem that 



u N {x) = J ... J j u(y) 

— OO L L-OO L“00 

dd N {y 2 - x 2 ) 



dO N (y 1 - Xi) 



dyi 



X — 



dy -2 



dy 2 



dy i 

dON(y n - Xn) 
dy n 

N 



dy n , 



where fov(cr) = JXi $N(s)ds and Sn{s) = J_ N e' ls ^d^ = 



sin 27 t Ns 



Note (compare with P4.3 and P 13.11| ) that fov(c x) — > 0(a), a £ R, 
and V7V |0jv(ct)| < const. Actually, setting A*, = 5 n (a) da 

for k £ Z+, we have that Afc does not depend on TV, | A* | J, 0 as k — * 

OO 

oo, moreover, X 2 k > —X 2 k+i and = - / 0 °° sin x dx = 

k—0 

Therefore, 9n ( a) — > 0(c) and |0 jv(ct) | < 2Aq. Then (as in the proof 



of Lemma 13.101 we should integrate by parts, apply the Lebesgue 
theorem and obtain that un(x ) — > u(x). 



□ 



17.7. Remark. The proof of the theorem containing, in par- 
ticular, the solution of Exercise P - 
into account the proof of Lemma 



4.3 (and P 13.11 1 shows (taking 



17.5 1 that the assertion of The- 



orem m is valid under more broad hypotheses: it is sufficient to 
require that the function u and n its derivatives d k u/dx±dx 2 . . . dxk , 
k = 1, . . . , n, be integrable in R". 



17.8. Remark. Of course, the assertion of Theorem 17.6 makes 



the sense only for u £ L 1 tl C. However, this necessary condition is 



not sufficient for validity of (17.241, as it follows from 



17. 9. P. Constructing (compare with 
<Pn £ L 1 (R) D C(R), such that J^^ y 
N — > oo) and ||<£iv||ii + llvdvllc < 1 
Steinhaus theorem (see note 1 in Section 16), show that there exists 
a function ip £ L 1 ( R) n 0(11), for which equality (17.24) does not 
hold at least at one point x £ R. 



) a sequence of functions 
sin NyipN(y)dy — *■ oo (as 
and applying the Banach - 



Formal expression ( 17.21 1 and Theorem 17.6 suggest the idea on 



expediency of introducing the transformation 



i 



£>u> — >F _ 1 «eC, 
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defined by the formula 



(F 1 m)(t) = / e lx ^u(f)df, i = 2ni, i = \f—i, x £ 



(17.25) 



K" 



This formula differs from formula ( 17.22 ) by the sign of the expo- 
nent. The transformation F _1 is called the inverse Fourier trans- 
form, since u = F _1 Fu, if u £ VF 1,7l (]R") and Fu £ T 1 (R). 

Define (following L. Schwartz |54j ) the space of rapidly decreas- 
ing functi ons S = 5(K n ) C W 1,n (R"). In the space S (see Theo- 
i below), the transformations F _1 and F are automorphisms 



rem 



17.16 



(i.e., linear invertible mappings from S onto itself). 



17.10. Definition. The elements of the space 6>(R n ) are the 
functions u £ C 100 ^ 11 ) that satisfy the following condition: for any 
multiindices a = (or, . . . , a n ) and (3 = (/3i , . . . , /3 n ), there exists a 
number C a p < oo such that \/x = (x ±, . . . , x n ) £ R n 

„ „ r)\B\ 

I x a d%u{x)\ < C a p , where x a = xf 1 . df = ^ . 

In this case, we say that the sequence of functions Uj £ S converges 
in S to u ( Uj — > u in S) as j —> oo, if Ve > 0 Vm £ N 3 j Q £ N 
Vj > jo the following inequality holds: p m (uj — u) < e, where 



Pm(v) = SUp 
aiGR 77 



£ (i+Nri^ 



, |a|<r 



Obviously, e “ 2 £ 5(R) but e x2 sin(e x2 ) ^ <S(M). 



17.11. P. Check that the space S is a Frechet space (see Re- 
mark 16.21) in which the distance p can be defined by the formula: 

OO 

p(u,v) = d{u — v), where d(<p) = 2~ m inf(l,p m (y>)). 

771=1 



17.12. P. Show (see P \T6J3\ ) that V(R n ) C S(R n ) C £(K n ). 
In particular, show that the convergence in T> (in S) implies the 
convergence in S (in £). Verify that T> is dense in S and S is dense 
in £. 



17.13. P. Integrating by parts, verify that the following lemma 
holds. 
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17.14. Lemma. For any multiindices a, (3 and any u £ S 

(-?)^F[d“(^u(z))](£) = 5 ( 0 , u = F u. (17.26) 

17.15. Corollary. F<S c S, i.e., F u £ S, if u £ S. 



PROOF. Since u £ S, it follows that for any fixed N £ N and 
any a, (3 £ Z” there exists d a p > 0 such that \df{x l3 u{x))\ < d a p{ 1 + 
|a;|)~ Jv . Therefore, by virtue of Lemma 

|r^n(0l < ||F[9“(x^)]||c < D a p sup | df (x^u) \ . (17.27) 



17.14 



Thus, u £ S. □ 

17.16. Theorem. The mappings F : S — » S and F _1 : S — > S 
are reciprocal continuous automorphisms of the space S. 



PROOF. F is linear and, by Theorem |17.6| monomorphic. Let 
us check that Vu £ S 3u £ S such t hat F u = u. Let u 0 = F u. 
Since uq £ S, according to Theorem 



17.6 



= F -1 Fu = 



F _1 u 0 . 



Consider the function u(x) = u o(—x). We have u = F _i uo = Fu. 
Definition 17.10 immediately implies that F Uj — > 0 in S, if Uj — > 0 
in S. The same arguments are valid for F -1 . □ 



17.17. Lemma (the Parseval equality). Let f and g £ 6>(R n ). 
Then 

(F/,F 5 ) i 2 = (/, 5 ) i 2 , i.e., J f(0g(0d£ = j f(x)g(x)dx. 

M” K" 

(17.28) 

Moreover, 

(Ff,g) = (f,Fg), i.e., j f{f)g{f)df= j f{x)g{x)dx. (17.29) 

R™ R n 



Proof. The Fubini theorem implies ( |17.29 ) 



since 



J f{x)g(x)dx = J J f(x)e lxi g{£) dx df = J f(Qg(£) df 

R 71 R n R n R n 



Let h = Fg. Then g = Fh, since 
g{f) = (F _1 /i)(£) = J e^h{x)dx = J e~° l ^h{x)dx = (F 
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Substituting g(£) = h(£) and g(x) = h{x) into (17.291, we obtain 
(Ff,Fh ) L 2 = (f,h) L 2 V/i. G S, i.e., (up to the notation) (17.281. □ 

Note that both sides of equality ( |17.29 ) define linear continuous 
functionals on S: 



/ : 5 9 5 i — > J f(x)g(x) dx, f :S 3 g i — > J f(Qg(Q df. 

In this connection, we give (following L. Schwartz [54]) two defini- 
tions. 



17.18. Definition. 6> , (M n ) is the space of tempered distribu- 
tions , i.e., the space of linear continuous functionals / : <S(K") — > C 
equipped with the operation of differentiation 

(5“/,^) = (-l)l“l(/,9», 

where a G Z™, and with the operation of multiplication ( af,ip ) = 
(/, atp) by any tempered function a, i.e., a function a G C'°°(K 11 ) sat- 
isfying the condition Va 3 C a < oo 3 N a < oo, such that |<9 Q a(:r)| < 

C , a (l + |a|) J V 

17.19. Definition. Let / G S', g G S' . Then the formulae 

(F f, tp) = (/, Fip) WipGS and (F _1 gr, = ( g , F" V) Vi/>g5 

(17.30) 

specify the generalized functions f = Ff G S' and F^ 1 ^ G S' , which 
are called the Fourier transform of the distribution f G S' and the 
inverse Fourier transform of the distribution g G S ' . 



17.20. Example. It is clear that d G S', 1 G S'. Find F<5 and 
FI. We have 

(Fd,tp) = {6, Ftp) = £>(0) = lim J e~ lx ^ip(x)dx = J tp(x)dx = (l,<p), 

i.e., F5 = 1. Similarly, F _1 i5 = 1. Furthermore, (Fl,y>) = (l,Fy>) = 
(F” 1 ^, F(^) = (<y,F _1 F V ?>, he., FI = 6. Similarly, F _1 l = S. 

l ) c S'(R n ) c 



V'l 



17.21.P. Verify (compare with P 17.12) that £'( 
R n ). 



17.22.P. Prove (compare with P \16.1E\ and [57 ] ) that f G < S , (M”) 
if and only if there exists a finite sequence {/a}| a |<iv of functions 
f a G C(R") satisfying the condition |/ a (ai)| < C( 1 + |a:|) m , where 
C < oo, m < oo and such that f = (F d a f a - Thus, S' C D' . 

\a\<N 
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17.23. P. Verify that the mappings F : S' — * S' and F -1 : S' — > 
S' are reciprocal automorphisms of the space S' which are continuous 
relatively the weak convergence in S' , i.e, if v — > oo, then (F f v , ip) — ► 
(F /, p)VpeS^ (/„, ip) - (/, Vp G 5. 



17.24. P. Calculate F5 and FI (compare with Example 17.20), 
by considering the sequence of functions 

5„(x) = v ■ l[_i/„ ) i/i,](a;) and l v (x) = 1 (a;), 

where l\ a ,b\ — 0{x — a) — 6{x — b). 

17.25. P. Considering the sequence of the functions 

f„(x) = 0±{x)e Tx / v , ieK, 

show ( see (12.7 1 ) that 0±(£) = ±- . 

± 0 

17.26. Remark. The space S' is complete with respect to the 
weak convergence, since S is a Frechet space (see P 17.11 and Re- 
mark 16.271. 



17.27. P. Show that formula (17.261 is valid for any u G S'. 

17.28. Lemma. LetC f G Then f = F f is a tempered 

function ( see Definition 17.lS\ ) and 

7(0 = </(*), e" Sa *>- 



(17.31) 



7) Recall that the space E' is defined in P 
PROOF. By virtue of Theorem 
f a G Co(R"). Therefore 



16.13 



16.20 



/ = E da fa, where 

|o:|<m 



<mp@) = £;<dzfa(x),( F v’)(*)) 

a 

a 

= £(*)'“' (/.M.FIfWOK*)) 



= E(<) 



M 



fa (x) 






dx. 
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Since f a (x)e x R?)j it follows (by the Fubini 

theorem) that 






J2Ci) H MxV x k a dx 



<p(Z)<% 



= J (J2Ux),(d x ) a e-°^M0d( 



= I (f(x),e~^MOdC 



Thus (by virtue of P 13.19), /(£) = (f(x),e zx ^}. Similarly, 

d p f{0 = {f{x),{-ix) f, e-°'*). (17.32) 



Since / £ S' C S' , it follows that (by virtue of Definition 17.10) 
37V > 1 such that 

\(f(x),ip{x)}\ < N sup ^2 {l + \x\) N ■\d a i/j(x)\Vilj £S. 

XeR " |q;|<JV 

Therefore, |c^/(£)| = \(f(x),a(x)(— ix)^e~ lx ^)\ < (7(1 + \£,\) N ■ □ 



18. The Fourier Laplace transform. The Paley Wiener 

theorem 



Formula (17.261 proven in Section 17 (which is valid for u £ S', 
see P 17.27 1 contains an important property of the Fourier transfor- 
mation, which is often expressed in the following words: “after ap- 
plying the Fourier transformation the derivation operator becomes 
the multiplication by the independent variable”. More exactly, the 
following formula holds: 



f(d“«(®)) = rs(o, 



(18.1) 



where = ( i)~^ d and u = Fit, u £ S’ . 

Property ( |18.1 1 allows us to reduce, in a sense, problems involv- 
ing linear differential equations to algebraic ones. Thus, applying 
the Fourier transformation to the differential equation 

A(D x )u(x) = ^2 a a D“u(x) = f(x), a a £ C, f £ S', (18.2) 

|a|<m 
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we obtain an equivalent “algebraic” equation 



■ 5(0 = £ a ^° I • “(0 = /( 0, / 6 5'. (18.3) 

\ \oi\<m 

As Hormander [29] and Lojasiewicz IS] have proved, equation (|18.3[ ) 
has always a solution u £ S' (see in this connection Remark 19.2 1. 
Therefore, the formula u = F _1 u determines a solution of differential 



equation (18.2). Indeed, by virtue of (18.1) and Theorem 17.16 

/(®) = f- 1 (7(0) = f- 1 (a(0-5(0) 



= A{D X ) F- 1 ^)) = A(D x )u(x). 



This approach to construction of a solution of a linear differential 
equation with the help of the Fourier transformation is rather similar 
to the idea of the operational calculus 1 ) (see, for instance, [38]) using 
the so-called Laplace transformation (introduced for the first time 
(see j34j ) by Norwegian Niels Abel who was a mathematical genius, 
little-known while alive, and who died from consumption when he 
was 26 years old). The Laplace transformation transfers a function 
f oft £ R + integrable with the “weight” e st for any s > 0, into the 
function 

OO 

L[f](s) = J e~ st f(t)dt, s> 0. (18.4) 

o 



Let us illustrate the idea of the calculus of variations on the example 
of problem (6.11 1, restricting ourselves by the case a = 0. In other words, 
we consider the problem 



Ut = u x 



t > 0, |x| < 1; u 



— ±i 0; 



= l. 



(18.5) 



As has been noted in Section 17 series ( 17.16 1 constructed by the Fourier 
method, which gives the solution of this problem, converges very slowly 
for small t. By the way, this can be foreseen, since the Fourier series 
converges slowly for discontinuous functions, and the function u(x,t) is 
discontinuous at the angle points of the half-strip {|x| < 1, t > 0}. In this 
connection, consider preliminarily the problem 



dT 

Iff 



d 2 T 



0 Q£2 1 £ > 0 ’ 



> 0 ; T 



e=o 



= Ti; T\ 



= To, 



(18.6) 



which simulates the distribution of the temperature in the vicinity of an 
angle point. 
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Passing (see Section [6| to the dimensionless parameters in the stan- 
dard way 

r = Z/V^r, u = (T — Ti)/(T 0 - Ti), 



from (18.61 we obtain that u (t, £) = f(r), where the function / satisfies 
the following conditions: 

f”(r) + 7 ,/V) = /(0) = 0, /( oo ) = 1. 

Hence, w(t, £) = erf(£/(2 y/ar)) = 1 — erfc(£/(2 y/ar)), where 



y 

erf(y) = ^ J e ~ V dr h erfc (y) = 1 - erf (y). 
o 

These preliminary arguments suggest that the solution u(x,t) of problem 
(18.5 (for small f, seemingly, must be well approximated by the following 



1 — [erfc((l — x)/2 Vt) + erfc((l + x)/2 Vi)]. (18-7) 

This point allows us to obtain the representation of the solution of problem 
( 18.5 1 in the form of a series rapidly converging for small t with the help of 
the Laplace transformation. Denoting the function L[u { ■, *)](s) by v(s, x), 
where u is the solution of problem ( 18.5 1, we rewrite, taking into account 
the two following obvious properties of the Laplace transformations: 

L[l](s) = 1/a, L[f'](s) = s ■ L[f](s) - /( 0), (18.8) 

problem ( 18.5|l in the form (“algebraic” in the variable s) 



( s ■ v(s,x) - 1) - v xx (s, x) = 0, v(s, x)\ a 



=±i 



= 0, s > 0. 



This problem can be solved explicitly. Obviously, its solution is the func- 
tion 

v(s,x) = • ch (y/sx)/ ch(v / s). 



Thus, the solution u of problem ( 18. 5|) satisfies the relation 



L[u(-, a:)](s) = - — - • ch (yfsx)/ ch(v / s). 
s s 



(18.9) 



Formulae ( 18.7 1 and ( |18.9| ) suggest that in order to obtain the represen- 
tation of the solution of problem ( |18.5| ) in the form of a series rapidly 
converging for small t we should 

• first, find the Laplace transform of the function /*, (t) = erfc(£:/2\/t) 

• second, represent the right-hand side of formula ( 18. 9| ) in the 
form of a series whose members are function of the form L[fk]. 
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Below, it will be shown that 

L[f k ](s) = ^ exp(-fc%/i). (18.10) 

On the other hand, expressing ch via exp and representing (l + g) -1 , where 
q = exp(— 2y/s) < 1, as the series 1 — q + q 2 — q 3 + ■ • ■, we obtain that 



ch (y/sx) y - y 

s-ch^i ^ ’ 

v n = 0 



,+i exp[--/s(2n + 1 - a;)] + exp[—/s(2n + 1 + a;)] 



(18.11) 



Basing on ( 18.8 1 ( 18.11 1, one can show that the solution of problem ( 18.5 ) 
can be represented as a series 



u(x, t) = 1 + 



' N 

E(- 1 ) n+1 ^ 



_n = 0 






(18.12) 



where a n = erfc((2n + 1 — x)/2 yt) + erfc((2n + 1 + x)/2y/i) and rjv = 

E ( — l) n + 1 an- 

n>N 



18.1.P. Show that in (18.121 



M - ^\ln eX ^~ N2/t) - 



(18.13) 



18.2. P. By comparing estimate (18.13) with estimate dl7l9| , 
show that , for t < 1/4, it is more convenient to use the representation 
of the solution of problem (18.5) inform (18.121, but for t > 1/4, it 
is better to use that form (17.161. 



Let us establish formula (18.101. It follows from formula 

L [f'k\{s) = exp (-ky/s) (18.14) 



and the second formula in (|18.8|, since 

/fc(0)=0, 



fk(t) = 1/2 t 3/2 exp(— fc 2 /4t). 
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As for formula (18.14), it can be proved, taking into account P 18.3 



in the following way 

OC 

k r 



L[f' k ](s) = 



2^ , 



t -3/2. e -fc a /4t dt = 



oo 

J exp -[rj 2 



4 rj 



^}dri 



= -e~ k ^~ s 



OO 

J g-in-i/n ) 2 dr] = e - fc ^. 



k k 

(The change of variables: r] = — y=, a = — >/s.) 

18.3. P. Let F(a) = / Q °° exp [— (rj — ^) 2 ]dr/, where a > 0. Then 
p = Tjl 

Hint. F'(a) = 0. 

There is a close connection between the Fourier and Laplace 
transformations. It can be found, by analyzing the equality 

d p m = </(*), f e £'(R n ). P £ Z» , 

which has been proved in Lemma 17.28| The right-hand side of this 



equality is meaningful for any complex £ £ C n and is a continuous 
function in C n . Thus, as is known from the theory of functions of a 
complex variable, the analytic function is defined 

/ : C" 3 £ i — ► /(£) = (f(x),e~* x £) £ C, 

which can be treated as the Fourier transform in the complex do- 
main. This function is sometimes called the Fourier- Laplace trans- 
form. This name can be justified by the fact that, for instance, for 
the function / = 6 + f £ L 1 (K) (compare with Example 17.31, the 
function 

OO OO 

C_ 3 £ i — * J e~ lx ^ f(x)dx = J e~ lx ^ f(x)dx £ C, 

— oo 0 

being analytic in the lower half-plane C_, is, for real £ (respectively, 
for imaginary £ = —is/ 27t, where s > 0), the Fourier transform 
(respectively, the Laplace transform) of the function /. 

The important role of the Fourier-Laplace transformation con- 
sists in the fact that, due to so-called Paley-Wiener theorems, certain 
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properties of this analytic function allow to determine whether this 
function is the Fourier-Laplace transform of a function / as well as 
to characterize the properties of this function /. In Section [22] we 



shall use (compare with Example 17.3 I 



18.4. Theorem (Paley-Wiener) . Let f be analytic in C_ and 



sup / |/(£ - iv)\ 2 d£ < oo. 

i)>0 J 



Then f is the Fourier transform in C_ of the function f = 6+f £ 
L 2 (R). 

The proof of this theorem and of the inverse one see in HD- 

19. Fundamental solutions. Convolution 

At the beginning of Section [18] it was noted that the differential 
equation 



A(D x )u(x) = ^2 



a a dfu(x) = f(x), a a £ C, / £ £' , 



(19.1) 



\ot\<r, 



has a solution u £ S'. In contrast to equation (18.2), the function / 



in ( 19.1 ) belongs to £' C S'. This fact allows us to give an “explicit” 



formula for the solution of equation (19.1 1, in which the role of the 



function / is emphasized. In this connection note that the formula 
>exp(— q\x - y\) 



< x ) = 4 ^ f f(v)- 



\ x — y\ 



-dy, q>0, fe£'nPC b , (19.2) 



giving (compare with (7.8 1 (7.9 1) a solution of the equation — A u- 



q 2 u = /, is meaningless for q = 0 , if supp / is not compact (for 
instance, f = 1). 

In order to deduce the desired “explicit” formula for the solution 



u of equation (19.1 1, first, we represent the function u = F u in the 
form u(£) = /(£)e(£), where e £ S 1 is a solution of the equation 



the function u = F 1 (/-e) via / = F i / and the function e = F e 
satisfying (by virtue of the relation A(£)u(£) = 1) the equation 

A(D x )e(x) = S(x). (19.3) 



A(£) • e(£) = 1 (see Remark 19.2 below). Then, it remains to express 
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19.1. Definition. A function E £ V' is called a fundamental 
solution of the operator A(D X ), if A(D x )E(x) = S(x). 

19.2. Remark. Any differential operator with constant coeffi- 
cients has (as has been proved in [29i . 41]) a fundamental solution 
in the class S'. However, the presence of the space V in Defini- 
tion [T9TT] is justified by the fact that, for some differential operators, 
it is possible (as has been shown by Hormander) to construct in V 
a fundamental solution locally more smooth than the fundamental 
solution in S' . (Note that the two fundamental solutions Ei and E 2 
of the operator A(D X ) differ by a function v = Ei — E 2 satisfying 
the equation A(D x )v = 0.) 

If A(£) 7^ 0 for any £ £ R n , then the formula E(x) = F _1 (l/A(£)), 
obviously, determines a fundamental solution of the operator A(D X ). 
In this case, E £ S' , because 1 /A(£) £ S'. In the general case, the 
fundamental solution can be constructed, for instance, with the help 
of regularization of the integral f <p(£)d£/A(£) (compare with P 14.4) 
that can be most simply made for ip £ D, since in this case the regu- 
larization (by virtue of analyticity of the function <p = J?ip) is possible 
by passing of £ into the complex domain, where A(£) 7^ 0 (see, for 
instance, (57l ). 



19.3. Examples. It follows from P 7.1 that function (7.9) is the 



fundamental solution of the Laplace operator. According to P 12.8 



the function 9(t — |x|)/2 is the fundamental solution of the string 
operator. For the heat operator dt — d xx , the fundamental solution 



is E(x,t) = 9(t)P(x,t) 1 where the function P is defined in (6.15l. 



Indeed, by virtue of the properties of the function P proven in Sec- 
tion [gJ for any ip £ 7?(K 2 ), we have 

OO 

(E t - E xx ,ip) = -{E, (p t + (p xx ) = - £ li j \<Pt + p xx )Edxdt 



lim 

e -++0 



OO OO 

J J {P t - P xx )p dx dt + J P(x,t)tp(x,t)dx 






x—t—0 ’ 



Before representing the solution u = F 1 (/ • e) of equation (19.1) 



via f £ S' and the fundamental solution e £ S' of the operator 
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A(D X ), we express the function F _1 (/ • g) via / and g in t he as- 
sumption that g = F^ 1 ]/ £ £' . According to Theorem 16.20 there 
exist multiindices a and (3 and continuous functions f a and gp with 
compact supports such that / = Dff a , g = D%gp. Therefore, 

F ~\m ■ g(0) = F _1 r +/3 / a (?) • MO) = D^F~\f a ■ gp). 

Let us calculate F -1 (/ Q • gp). We have 



f a - 9f3 = 



a -i(y+i t 



K fa(y)g/3(cr) dy da. 



Setting y + a = x and using the compactness of supp/ Q and Slippy, 
we obtain f a -gp = F [f a * gp\, where tp * ip denotes (see note 4 in 
Section [ 5 ]) the convolution of the two functions tp and ip, i.e., 

(<P * ip)(x) = J <p(y)fp(x - y)dy. 

19.4. P. Verify that if tp £ . ip £ then tp * ip = ip * tp, 

D a (tp * ip) = (D a tp) * ip = tp * (D a tp). 



19.5. Definition. (Compare with P 19.4 1 . The convolution f*g 
of two generalized functions / = D a f a £ £' and g = O fj gp £ £' , 
where f a , gp £ Co(R ra ), is the generalized function D a+I3 (f a * gp). 

Note that for f £ £' and g £ £ ' the following formula holds: 

f ~ 1 (7-g) = f*g- (19.4) 

Indeed, F[D° +/3 (/ Q * gp)] = £° +/, F (/„ * gp) = (Ofa)(f%)- 
Let us find F _1 (/ • g) in the case when f £ £' and g £ S' . 

19.6. P. Let f £ £', g £ S' . Setting g v (x) = tp(x/v)g(x), where 
tp £ C^°(K"), tp = 1 for |x| < 1, verify that g„ £ £' ; g v g in S' ; 
7 -9v -> f-g in S', and (f*g„,<p) -> (F _1 (/ • g),tp) \/tp £ S. 



On the base of P 19.6 generalizing formula (19.4), we give 

19.7. Definition. The convolution f*g of two generalized func- 
tions / £ £ ' and g £ S' is the function from S' , given by the formula 
f*g = F _1 (/ • g). (Note that / £ C°°, g £ S'.) 

19. 8. P. (Compare with P\19.JA). f £ £' , g £ S' =>■ f * g = g* f ; 
D a (f * g) = (D a f) *g = f* (D^g); 6*g = g. 



It follows from abovesaid that the following theorem holds. 
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19.9. Theorem (compare with (19.21 and note 4 in Sec t. |5| . 
The desired “explicit” formula for the solution of equation (19.1) 
has the form u = f * e, where e £ S' is the fundamental solution of 
the operator A(D X ). 



19.10.P. Prove the Weierstrass theorem on uniform approxima- 
tion of a continuous function f £ C{K) by polynomials on a compact 
I\ C R ra : Ve > 0 3 a polynomial p such that \\f(x) — p{x)\\c(K) < e - 

Hint. Let 9 be a neighbourhood of K. Following the scheme 
of the proof of Lemma |13.10| and taking into account P |4.4| set 

P ( x ) = J f{y)5u{x - y)dy, 

n 



where 



sa x ) = n 



m—1 



2 \v' s 



n T V 



X = (xi, . . -,X n ). 



19.11. Lemma. Let u £ L 1 f)C, v £ L 2 . Then u* v £ L 2 and 

||m * v\\ L 2 < ||m|| l i • ||i>||l 2 - (19.5) 

Proof. |/u(£ - rj)v(r])dy\ 2 < (f \u(£ - ri)\dij) ■ A(£) = || u|| £ i • 
H(£), where T(£) = / |u(£— rf)\-\v{if)\ 2 dip However, (see Lemma 8.26 1 
f A(£)d£ = \\v\\ 2 L 2 ■ ||w||ii. □ 



20. On spaces H s 

The study of generalized solutions of equations of mathematical 
physics leads in a natural way to the family of Banach spaces W p,m 
introduced by Sobolev. Recall Definition 17.4 of the space W p,m (TL). 
For p > 1 and m £ Z+, the space W p ’ m (VL) is the Banach space of 
the functions u £ L p (Ll) whose norm 



IMIiTp.m(n) — 




( 20 . 1 ) 



20. ON SPACES H a 



103 



is finite. Here, d a u is the generalized derivative of the function u, 
i.e., 



d a u = v £ L] oc (VL) and 

J v ■ ipdx = (— l)l“l J ud a ipdx (20.2) 



S.L. Sobolev called the function v satisfying conditions (20.2) the 
weak derivative of order a of the function u. Maybe, this is the 
reason, why letter W appeared in the designation of the Sobolev 
spaces. 

For p = 2, the spaces W p,m are Hilbert spaces (see note 4 in 



Section 17 1. They are denoted (apparently, in honour of Hilbert) by 
H m . These spaces play a greatly important role in modern analysis. 
Their role in the theory of elliptic equations is outlined in Section [22] 
A detailed presentation of the theory of these spaces can be found, 
for instance, in P HO]; a brief one is given below. 



20.1.P. Using formula (18.1 1, verify that the space H m { R") in- 
troduced above for m £ Z + is the space of u £ iS^R”) such that 
(1+ \f\) m (Fu)(0 £ L 2 (R n ). 



20.2. Definition. Let s £ R. The space H s = H s ( R") consists 
of u £ S' = S'(R n ) for which the norm 

IMU = ll(£) S • «(0IU a (R»)’ where (£) = 1 + |£|, u = Fu. 
is finite. 



20. 3. P. Verify that if a > f3, then S C H a C H@ C S' and the 
operators of embedding are continuous and their images are dense. 

20.4. Theorem (Sobolev embedding theorem). If s > n/2 + m, 
then the embedding H s ( R”) C C™( R") holds and there exists C < oo 
such that 

NI( m) <C'|M| s Vue7r, (20.3) 

where 

IMI(m) = X! SUP \ d ° U ( X )\- 

i | xGM 71 

Proof. We have to prove that Vu £ H s { R") 3v £ C™{ R n ) such 
that v = u almost everywhere and ||f||( m ) < C||m|| s . It is sufficient 
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to prove it for m = 0. The inequality 



HO 1 = 



5 ( 0 ( 0“ -(0 



~ s e zx ^df 



< Hull 



(0- Js d£ 



1/2 



where it £ S (M) , implies estimate (|20.3| ). If it € H a , u n £ S and 
||it n — it|| s — ► 0, then by virtue of (20.3) 3v £ C° such that ||u n — 
u||( 0 ) -*• 0 and ||i>||( 0 ) < C'||u|| s . Since ||u- v\\ L 2 ( Q) < Cfi(||u-ii„|| s + 
||it n — u||( 0 )) — ► 0, we have it = v almost everywhere. □ 

20.5. P. Let u(x) = ip( 2x) In | In |a:||, where x £ R 2 and, ip is the 



function from Example 3.6 Show that u £ i(f 1 (R 2 ). Thus, iJ"/ 2 (R") 



cannot be embedded into C(R"). 



20. 6. P. Verify that 6 £ H s (M . n ) for s < —n/2. 



20.7. Theorem (of Sobolev on traces). Let s > 1/2. Then, 
for any (in general, discontinuous) function u £ H s (S. n ) the trace 
qit £ jJ s_1 / 2 (]R" _1 ) is defined, which coincides, for a continuous 
function u, with the restiction it | _ 0 of the function u onto the 

hyperplane x n = 0. Moreover, 3 C < oo such that 

Ml'_i/ a < qi«||. Vu £ H s (R n ), (20.4) 

where || ■ ||/. is the norm in the space H a (W l ~ 1 ). 



Proof. Let x = ( x',x n ) £ R" 1 x R. For it £ 5(11"), we have 
u(x', 0) = /,£„_! [/ R u(£',£„)d£„] df; therefore, {hu\\' s _ 1/2 ) 2 = 

J^-AC) 28 - 1 \J R d(CAn)df n \ 2 dC- Then 



5(0,6,)d& 



<A(0 J (0 2S \um 2 dfn 



where A(f) = J(0~ 2s df n < C S (A)~ S+1/2 , and C s = Cf( 1 + 
z 2 )~ s dz < oo for s > 1/2; (z = £ n (l + |£'| 2 ) _1 / 2 ). Therefore, 
IIt u II«_i ./2 — C|MU f° r u £ S. If it £ H s (R n ) and ||u„ — u\\ s — > 
0 as n — > oo and u n £ S, then 3w £ TP _1 / 2 (R" _1 ) such that 
||7it„ — w \\' s -i /2 w does not depend on the choice of the se- 

quence {un\. By definition, yit = w; and estimate (20.4) is also 
valid. □ 



20.8. Definition. The operator P = Pq : X>'(R" ) 9 / 1 — * 
Pf £ V(Li), where is a domain in R" and (Pf,ip) = ( f,ip ) \/ip £ 
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V(Q), is called the restriction operator of generalized functions given 
in R n onto the domain Cl. 



20.9. Definition. Let H s (Cl) denote the space PnH s (M. n ) with 
the norm 

||/||s,n = inf ||L/|| S , / £ H S (C), (20.5) 

where the infimum is taken over all continuations Lf £ H s (M. n ) of 
the function / £ H s (Cl) (i.e., PnLf = /)■ If it is clear from the 
context that a function / £ H s (Cl) is considered, then we may omit 
index Cl write ||/|| s and instead of ||/|| s ,n. 



20.10. Definition. The space H S (T ), where T = dCl is the 
smooth boundary of a domain Cl <e R ra , is the completion of the 
space C 00 (r) in the norm 

K 

lKr = £M' a - (20.6) 

k = 1 



K 

Here, || • ||' is the norm of the space FP(R" _1 ), J2 Pk = 1 is the 

k= 1 

partition of unity (see Section |3| subordinate to the finite cover 

K 

(J Tfc = T, where T* = Cl k n T, and flfc is a n-dimensional do- 
fe= l 

main, in which the normal to T do not intersect. Furthermore, the 
function ip k p £ (^(R 71-1 ) is defined by the formula (< PkP)(y ') = 
Vk((?k {y')) ' p(&h w l iere a k is a diffeomorphism of R", (affine 

outside a ball and) “unbending” r^. This means that, for x £ Cl k , 
the ?rth coordinate y n = y n {x) of the point y = ( y',y n ) = cr k (x) is 
equal to the coordinate of this point on the inward normal to T. If 
it is clear from the context that we deal with a function p £ H S (T), 
then equally with ||p||' sr we also write ||p||'. 



20.11. Remark. Definition 20.10 of the space H S (T) is correct, 
i.e., it does not depend on the choice of the cover, the partition 
of unity and the diffeomarphism a k - In the book | 59 ] this fact is 
elegantly proven with the help of the technics of pseudodifferential 
operators. 



20.12.P. Let s > 1/2. Show that the operator C(Cl) n H s (Cl) 3 
u i — > u| r £ C(T) can be continued to a continuous operator 7 : 
H s (Cl) -> FP- 1 / 2 ^). 
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20.13. Remark. The function 7 u £ fP _ 1 / 2 (r), where s > 1/2, 
is called the boundary value of the function u £ H S (Q). One can 
rather easily show (see, for instance, [70]) that £P _1 / 2 ( T), where 
s > 1/2, is the space of boundary values of functions from H S (Q). 
The condition s > 1/2 is essential, as the example of the function 
u £ R 1 / 2 (R_|_) given (compare with P) 
ip( 2x) In | In |cc| | shows. 



20.5) by the formula u(x) = 



20.14. Remark. The known Arzela theorem (see [Mirfiji as- 
serts that if a family {/„} of functions f n £ C(Cl) defined in Q (= R n 
is uniformly bounded (i.e., sup ||/„|| < 00 ) and equicontinuous (Ve > 

n 

0 35 > 0 such that \f n (x) — f n (y)\ < eVn, if \x— y\ < 5), then one can 
choose from this sequence a subsequence converging in C(fl). With 
the help of this theorem one can prove (see, for instance, j8] (701 ) 
that the following theorem is valid. 



20.15. Theorem (on compactness of the embedding). Let Ll <s 
R", and a sequence {u n } of functions u n G H s ( R) ( respectively , 
ti„ £ H s (dil)) be such that ||wn|| s < 1 (respectively, ||Rn||s < 1/- 
Then one can choose from this sequence a subsequence converging in 
R*(R) (respectively, in R*(9R)/, ift<s. 



21. On pseudodifferential operators (PDO) 

The class of PDO is more large than the class of differential 
operators. It includes the operators of the form 

Au(x) = J K(x,x- y)u(y)dy, u £ 
n 

Here, K £ X>'(R x R n ) and K £C°°(VLx (M n \ 0)). If K(x, x-y) = 
E a a (x) ■ S^^x — y), then Au(x) = E a a (x)d a u(x). Another 

|a| <ra |ct| <m 

example of PDO is the singular integral operators |45| . However, 
the exclusive role of the theory of PDO in modern mathematical 
physics (which took the shape in the middle sixties |30L 1351 1661 ) 
is determined even not by the specific important examples. The 
matter is that PDO is a powerful and convenient tool for studying 
linear differential operators (first of all, elliptic one). 
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Before formulating the corresponding definitions and results, we 
would like to describe briefly the fundamental idea of the applica- 
tions of PDO. Consider the elliptic differential equation in R™ with 
constant coefficients 

a(D)u = ^ a a D a u = /. (21.1) 

\ot\<.m 

The ellipticity means that 

a m(0 = ^2 a “£“ ^ 0 for l£l ^ 

| Q! | =772 

This is equivalent to the condition 



1(01 = 



| q;| <m 



> C |C| m for | > M > 1. 



(21.2) 



Now prove the following result on the smoothness of the solutions 
of equation (21.1). If u £ H s ~ l and a(D)u £ FP _m f or some s, then 
u £ H s . Of course, this fact can be established by constructing the 
fundamental solution of the operator and a(D) and investigating its 
properties (see, for instance, |31] 1. However, instead of solving the 
difficult problem on mean value of the function l/a(f), where f £ R ra , 
(this problem arises after application the Fourier transformation to 
equation (21.1 1 written in the form F -1 a(f)Frt = /) it is sufficient 
to note “oniy” two facts. The first one is that, taking into account 
(21.21, we can “cut out” the singuiarity of the function 1/a by a 
modifier p £ C°° such that p = 1 for |f | > M + 1, p = 0 for |£| < M. 
The second is that 

(F -1 (p/a)F)(F -1 aF)it = u + (F _1 rF) 

Therefore, by virtue of obvious inequafities 

|p(0M0l<C(l + |e|)- m , \rm<C N (l- 

which impfy the inequafities 

||(F- 1 (p/a)F)/|| s < C\\f\\ s -m, IKF-VF)' 



r = p-1. 


(21.3) 


|) _Ar VJV> 1, 


(21.4) 


< C |u|| s _jv, 
e 


(21.5) 


, ueH s , 


(21.6) 



ll«ll- < C (ll/IU-rn + IMU-i) , / = a(D)u, 
where C does not depend on u. The above-mentioned resuft on the 
smoothness of the soiution of efliptic equation (21.1) follows from 
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(21.6 1 . Here the word “a priori” for estimate (21.61 of the solution 



of equation (21.1) means that the solution was obtained before the 



investigation of solvability of equation (21.1 1, i.e., a priori. 



The simplicity of the deduction of the a priori estimate (21.6) 



characterizes brightly enough the role of the operators of the form 
F _1 aF. Such operators are called pseudodifferential operators con- 
structed by the symbol a = a(£). We shall also denote them by 
Op( a (Q) or a(D). Depending on the class of symbols, one or another 
class of PDO is obtained. Ifa(x,£) = ^2 a a (a;)£“, then a(x, D)u(x) = 
Op(a(x, £))u(x) = a a (x)D%u(x). If a(x, £) is a function positively 
homogeneous of zero order in £, i.e., a(x, t£) = a(x, £) for t > 0, then 
a(x,D) — Op(a(x, £)) is the singular integral operator [45], namely, 



Op(a(x, £))u(x) = b(x)u(x) 



lim 

e — >0 



\x-y\>e 



c(x,x- y) 

-y \ n 



[ {y)dy- 



Here, c(x, tz) = c(x, z ) for t > 0 and /|q =1 c(x , z ) dz = 0. In partic- 
ular, in the one-dimensional case, when a(£) = a + 0 + (£) + a_6L(£), 
where 9 + is the Heaviside function and 6L = 1 — 0 + , we have 



( ( /-\\ a + " 1 “ a ~ / \ i 

Op(a(x, €))u = — — — u ( x ) + 2 tt V ' p ’ 



a+ — a_ 

'« 

x-y 



{y)dy , 



that follows from P 17.25 and (12.7 1. 



Let us introduce a class of symbols important in PDO. 



21.1. Definition. Let m £ R. We denote by S m = 5 m (K") the 
class of functions a £ C 00 ^" x K") such that a(x, £) = a,o(x,£) + 
ai(£) and Va V/3 3C a @ £ <S(R") 3Cp £ K such that 



\d^a 0 (x,0\<C a p(x)-(O m ~ m , 

|5?o 1 (0|<C f /3 <0 Tn " l/J| ,<0 = 1 + 10- (21-7) 

If a £ S m , then the operator a{x , D ) given by the formula 

Op(a(x, £))u(x) = J e lx ^a(x,t;)u(t;)dl;, u = Fu ( 21 - 8 ) 

and defined, obviously, on Cg°, can be continued to a continuous 
mapping from H s ( R n ) into H s ~ m (R. n ), as the following lemma shows. 
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21.2. Lemma (on continuity). Let a £ S m . Then\/s £ M 3 C > 0 
such that 



\a{x, D)u\\ s - m < C||«|| B Wu £ C 0 °°(M”). 



(21.9) 



Proof. If a(x,£) = ai(£), then estimate (21.9 1 is obvious. There- 
fore, it is sufficient to establish this estimate for a = ao (see Defini- 
tion 21.1 1. Setting A 0 v = a 0 {x,D) v, we note that 



0M(0 = 



e l ( x >S ^ao(x,r])dx) v(rj)dr]. 



By virtue of (21.1 1 and Lemma 17.14 we have 



I0M(0I < c a H Hv)\dv, M » 1 . 

The triangle inequality |£| < \t)\ + |£ — rj\ implies the Peetre inequality 

( 0 s < (v) s (ti-v) lsl - (21.10) 

Therefore, (C) s_m |(^ov)(C)l < C a /(?7) s (£ — ry) l s— ^1 — 1“! \v(ri)\dr}. It 
remains to apply inequality (19.51. □ 



21.3. Example. Let a(£) = e + l/(|£| 2 + q 2 ), e > 0, q > 0. Then 
a £ S° for e > 0 and a £ S~ 2 for e = 0. Moreover, (compare with 
( 19.2 1) for n = 3, a(D)u(x) = eu(x) + 7r f R3 \x — y I” 1 exp(— 27rg|a; — 
y\)u(y)dy, u £ C£°. Indeed, set / = Op( 1/(|£| 2 + Q 2 )) u that is 
equivalent to it = {\D\ 2 +q 2 )f, i.e. , —Af+(2nq) 2 f = 47 t 2 u. By virtue 
of the estimate ||/|| s < C||ti|| s + 2 , the solution of the last equation 
is unique in H s . It ca n be represented in the form / = 47 t 2 G * u, 
where (compare with P7.ll G(x) = exp(— 27rg|a;|)/47r|a;| £ H° is the 
fundamental solution of the operator —A + (27 iq) 2 . 

21.4. Definition. Let a £ S m . An operator Op(a(x ,£)) is 
called elliptic if there exist M > 0 and C > 0 such that 

Kar.OI >C|C| m Viel" and |£| >M 



(compare with (21.2 1). 



21.5. P. Following the above proof of estimate (21.6 1 and using 



Lemma 21.6 prove that, for the elliptic operator Op(a(x, £)) with 
the symbol a £ S m the prior estimate 



IMU < C(||a(a;,D)|| s _ m + ||u|| 



,- N ) \/u £H s , C= C(s, N), N>1, 

( 21 . 11 ) 



110 



3. THE SPACES H‘ . PSEUDODIFFERENTIAL OPERATORS 



holds. 



Hint. Setting R = Op(p(£)/a(x,£)), where p £ C°°, p = 1 for 
|£| > M + 1, p = 0 for |£| < M, show that R ■ Op(a(x,f))u = 
u + Op(r(x, f))u, where r £ S™ -1 . 

21.6. Lemma (on composition). Let a £ S k , b £ S m . Then 
VN > 1 

a(x,D) ■ Op(b(x,£)) = ^ Op £)(D*b{x, £)] /a! + T)v, 

|a|<IV 

where ||Tjvu|| s+iV -(>+m) < C\\v\\ s Mv £ H s . 

The proof see, for instance, in [551 . 

21.7. Definition. An operator T : Cq° — > S' is called smooth- 
ing , if ViV > 1 Vs £ R 30 > 0 such that ||Tu|| s+ jv < C||u|| 5 Vu £ H s . 

21.8. P. Let there be a sequence of functions aj £ S mj , where 

mj J. — oo as j | +oo. Then there exists a function a £ S mi such 
that {a— a j ) e S mN \/N > 1. 

j<N 

Hint. Following the idea of the proof of the Borel theorem 1 15. 2 [ 
one can set 

OO 

a{x,0 ='^2tp(£/t j )a j (x,t), 

3 = 1 

where ip £ 0 00 (K rl ), tp(£) = 0 for |£| < 1/2, tp(£) = 1 for |£| > 1, 
and choose tj tending to +oo as j — > oo so rapidly that, for \x\ < 1 
and M + |/3| + l<j, the following inequality holds: 

\d^D^/t J )a J (x,0\<2-H0 rn3 ~ 1 

Solution see, for instance, in [ 59] , 

21.9. Definition. All operator A : Oq° — > S 1 is called pseudo- 
differential of class L, if A = Op(a(x,f)) + T, where a £ S m for 
some to £ R, and T is a smoothing operator. Any function a a £ S m 
such that (a a — a) £ S~ N \/N is called the symbol of the operator 
A £ L. 



21.10. P. Applying Lemma 21.6\ show that the operator A £ 
L has (compare with P 16.22 ) the pseudolocality property, in other 
words, if ip and if belong to C(f° and if = 1 on suppip, then ipAfl—if) 
is a smoothing operator. 
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21.11. Remark. The class L is invariant with respect to the 
composition operation (see Lemma 21.6 and P 21.8 1 as well as with 
respect to the change of variables. The following lemma is valid (see, 
for instance, mm- 



21.12. Lemma (on change of variables). Let a £ S m . Then the 
operator a(x, D x ) = Op(a(x, £)) in any coordinate system given by a 
diffeomorphism (affine outside a ball) a : x i — » y = cr(x), for any 
N > 1, can be represented in the form 



Y2 Op ip a (y,r})(dga(x,£) | ?= 



\a\<N 



^=ta'(x)ri; x=a ~ 1 (y) 



■ t n , 



(21.12) 



where t a'(x) is the matrix transpose to cr'(x) = da/dx, and p a (y,p) 
is a polynomial in p of degree < \a\/2 given by the formula 



Va(y, V) = —, D z ex P - °(. x ) - - x), rj) 

a\ 

Moreover, ||TAri;|| s+ [jv/2+i]- m < C||f|| s Vi> £ H s . 



I z=x,x=a~ 1 (y) * 



22. On elliptic problems 

In Section |5] we have considered (for some domains O) the sim- 
plest elliptic problem, namely, the Dirichlet problem for the Laplace 
equation. One can reduce to this problem the investigation of an- 
other important elliptic problem - the problem with the directional 
(the term skew is also used) derivative in a disk Q, = {(x,y) £ R 2 | 
x 2 + y 2 < 1} for the Laplace equation 

Am = 0 in fi, du/d\ = f on T = dfl, f£C°°{T). (22.1) 



Here, d/d\ = ( ad/dx — bd/dy) is the differentiation along the direc- 
tion A (possibly, “sloping” with respect to the normal to the bound- 
ary T). This direction depends on the smooth vector field 

er : T 9 s i — > cr(s) = (a(s),b(s)) £ M 2 , a 2 (s) + b 2 (s) / 0 Vs £ T. 



Let us identify a point s £ T with its polar angle ip £ [0, 27r] . If 
a (jp) = (cos p, — sin <p), then A = v is the (outward) normal to T; 
if cr(p) = — (sin ip, cos ip), then A = r is the tangent to T. All these 
cases as well as others are important in applications. However, our 
interest in problem (22.1 1 is caused, first of all, by the fact that it 
well illustrates the set of general elliptic problems. 
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It turns out that the solvability of problem (22.1 I depends (see 



P 22.1 -P 22.4 below) on the so-called degree of the mapping a rela- 



tively the origin, namely, on the integer 

N = {arg[a(27r) + ib{ 2n)\ — arg[a(0) 



■ib( 0)]}/27 



It is clear that N is the number of revolutions with the sign, made 
by the point cr(tp) around the origin, when it moves along the closed 
curve a : [0, 27r] 3 ip i — » cr(ip) = (a(<p), b(<p)). 



For N > 0, problem (22.1 1 is always solvable but non-uniquely: 



the dimension a of the space of solutions of the homogeneous problem 
is equal to 2 N + 2. 



If N < 0, then for solvability of problem (22.1 1, it is necessary 



and sufficient that the right-hand side / is “orthogonal” to a certain 
subspace of dimension /? = 2|iV| — 1. More exactly, there exist (2|iV| — 
1) linear independent functions <l> 7 £ L 2 (r) such that problem (22.1) 
is solvable if and only if 



I f&jdT = 0 Vj = 1, . . . ,/3 = 2|7V| — 1. 

r 

In this case the dimension a of the space of solutions of the homo- 
geneous problem is equal to 1. A special case of problem ( |22.1 1, 
when A = v is the normal to T, is called the Neuman problem for 
the Laplace equation. In this case, N = —1, since a(<p) + ib(ip) = 
exp(— iip). The Neuman problem is solvable if and only if f r fdT = 0; 
and the solution is defined up to an additive constant. Indeed, if 
f r fdT = 0, then a continuous function 

S 

g{s ) = ff(so) + J f(<p)d<p. 

So 

is defined on T. Using the function g, we construct a solution v of the 
Dirichlet problem Av = 0 in U, v = g on T. Then the real part of the 
analytic function u+iv, i.e., the function u (defined up to an additive 
constant) is a solution of the considered Neuman problem, because 
du/dv = dv/d t = /, where d/dr is the differentiation along the 
tangent to T. Conversely, if u is a solution of the Neuman problem, 
then the “orthogonality” condition f r fdT = 0 holds. This follows 1 - 1 
immediately from the Gauss formula (|7.5|). Finally, the first Green 



